FREQUENCY  CONVERSION  OF  SINGLE  PHOTONS: 
PHYSICS,  DEVICES,  AND  APPLICATIONS 


A  DISSERTATION 

SUBMITTED  TO  THE  DEPARTMENT  OF  APPLIED  PHYSICS 
AND  THE  COMMITTEE  ON  GRADUATE  STUDIES 
OF  STANFORD  UNIVERSITY 
IN  PARTIAL  FULFILLMENT  OF  THE  REQUIREMENTS 

FOR  THE  DEGREE  OF 
DOCTOR  OF  PHILOSOPHY 


Jason  Scott  Pelc 
July  2012 


Report  Documentation  Page 


Form  Approved 
OMB  No.  0704-0188 


Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 
VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 


1.  REPORT  DATE 

JUL  2012 


2.  REPORT  TYPE 


3.  DATES  COVERED 

00-00-2012  to  00-00-2012 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


10.  SPONSOR/MONITOR'S  ACRONYM(S) 

11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 


4.  TITLE  AND  SUBTITLE 

Frequency  Conversion  Of  Single  Photons:  Physics,  Devices,  and 
Applications 

6.  AUTHOR(S) 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Stanford  University, Stanford, CA, 94305 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 


12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 


14.  ABSTRACT 

The  ability  to  manipulate  the  carrier  frequency  of  quantum  states  of  light,  through  a  process  called 
quantum  frequency  conversion  (QFC),  has  numerous  applications  for  both  technology  and  basic  science. 
For  example,  one  can  upconvert  a  single-photonlevel  signal  in  the  1.5-  m  telecommunications  band  (where 
single-photon  detection  has  been  challenging)  to  a  visible  wavelength  to  take  advantage  of  well-developed 
single-photon  detectors  based  on  silicon  avalanche  photodiodes.  On  the  more  fundamental  side,  the 
manipulation  of  a  single  photon’s  frequency  may  enable  the  construction  of  networks  of  dissimilar 
quantum  systems,  whereby  one  can  imagine  generating  many-body  entangled  quantum  states  over  vast 
distances. 


15.  SUBJECT  TERMS 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 

18.  NUMBER 

19a.  NAME  OF 

ABSTRACT 

OF  PAGES 

RESPONSIBLE  PERSON 

a.  REPORT 

unclassified 

b.  ABSTRACT 

unclassified 

c.  THIS  PAGE 

unclassified 

Same  as 
Report  (SAR) 

174 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


©  2012  by  Jason  Scott  Pelc.  All  Rights  Reserved. 
Re-distributed  by  Stanford  University  under  license  with  the  author. 


This  work  is  licensed  under  a  Creative  Commons  Attribution- 
Noncommercial  3.0  United  States  License. 
http://creativecommons.Org/licenses/by-nc/3.0/us/ 


This  dissertation  is  online  at:  http://purl.stanford.edu/zs827ns8173 


I  certify  that  I  have  read  this  dissertation  and  that,  in  my  opinion,  it  is  fully  adequate 
in  scope  and  quality  as  a  dissertation  for  the  degree  of  Doctor  of  Philosophy. 

Martin  Fejer,  Primary  Adviser 


I  certify  that  I  have  read  this  dissertation  and  that,  in  my  opinion,  it  is  fully  adequate 
in  scope  and  quality  as  a  dissertation  for  the  degree  of  Doctor  of  Philosophy. 

Robert  Byer 


I  certify  that  I  have  read  this  dissertation  and  that,  in  my  opinion,  it  is  fully  adequate 
in  scope  and  quality  as  a  dissertation  for  the  degree  of  Doctor  of  Philosophy. 

Yoshihisa  Yamamoto 


Approved  for  the  Stanford  University  Committee  on  Graduate  Studies. 

Patricia  J.  Gumport,  Vice  Provost  Graduate  Education 


This  signature  page  was  generated  electronically  upon  submission  of  this  dissertation  in 
electronic  format.  An  original  signed  hard  copy  of  the  signature  page  is  on  file  in 
University  Archives. 


Abstract 


The  ability  to  manipulate  the  carrier  frequency  of  quantum  states  of  light,  through  a 
process  called  quantum  frequency  conversion  (QFC),  has  numerous  applications  for 
both  technology  and  basic  science.  For  example,  one  can  upconvert  a  single-photon- 
level  signal  in  the  1.5-/irn  telecommunications  band  (where  single-photon  detection 
has  been  challenging)  to  a  visible  wavelength  to  take  advantage  of  well-developed 
single-photon  detectors  based  on  silicon  avalanche  photodiodes.  On  the  more  funda¬ 
mental  side,  the  manipulation  of  a  single  photon’s  frequency  may  enable  the  construc¬ 
tion  of  networks  of  dissimilar  quantum  systems,  whereby  one  can  imagine  generating 
many-body  entangled  quantum  states  over  vast  distances. 

Quantum  frequency  conversion  will  only  be  useful  if  it  can  be  done  both  effi¬ 
ciently  and  with  little  added  noise.  We  demonstrated  a  conversion  efficiency  exceed¬ 
ing  99.99%  using  reverse-proton-exchange  waveguides  in  periodically  poled  lithium 
niobate  with  approximately  150  mW  of  pump  power.  Noise  has  been  a  more  serious 
issue:  the  generation  of  noise  photons,  due  to  inelastic  scattering  of  light  from  the 
strong  pump  laser  used  to  drive  the  frequency  conversion,  has  limited  the  utility  of 
QFC  devices  in  many  applications.  We  present  an  analysis  of  the  two  primary  noise 
processes  in  QFC  devices  (spontaneous  Raman  scattering  and  spontaneous  paramet¬ 
ric  fluorescence),  and  offer  solutions  on  how  they  may  be  either  mitigated  or  avoided 
completely. 

We  then  discuss  applications  of  QFC  devices  for  up-  and  downconversion  of  single- 
photon  signals.  We  used  a  long-wavelength  pump  to  enable  high-efficiency  and  low- 
noise  single-photon  detection  for  1.5-/jm  telecom  band  signals,  and  demonstrated  a 
cascaded  frequency  conversion  approach  that  enabled  low  timing  jitter  as  well.  We 


also  demonstrated  a  downconversion  quantum  interface,  in  which  the  emission  from 
a  single  semiconductor  quantum  dot  at  a  wavelength  of  910  nm  was  downconverted 
to  1.56  n m  while  maintaining  the  single-photon  character  of  the  light.  The  results 
presented  in  this  dissertation  indicate  a  promising  future  for  QFC  devices  as  the  field 
of  quantum  communications  matures. 
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Chapter  1 
Introduction 


Computer  networks  form  the  basis  of  the  internet  and  today  represent  our  primary 
means  of  communication.  They  are  based  on  binary  switches  using  digital  logic 
in  semiconductor  devices,  and,  even  with  the  rapid  rate  of  miniaturization,  are  at 
their  heart  classical  devices.  That  is,  the  state  of  an  individual  logic  device  is  either 
a  “1”  or  a  “0”.  Quantum  information  processing  enables  superpositions  of  these 
binary  states  to  evolve  at  the  same  time,  allowing  a  degree  of  parallelism  not  possible 
in  classical  information  processing  [1],  It  is  known  that  a  quantum  computer  could 
solve  certain  classes  of  problems  exponentially  faster  (in  terms  of  computational  steps) 
than  their  classical  counterparts.  Although  a  large-scale  quantum  computer  remains  a 
distant  goal,  those  features  which  enable  quantum  computation  (namely  entanglement 
and  coherent  superposition)  are  the  very  ones  that  have  given  quantum  mechanics 
its  reputation  for  being  counterintuitive  and  bizarre.  As  such,  the  exploration  of 
quantum  computer  primitives — the  controllable  quantum  two-level  systems  and  their 
interactions — makes  a  very  worthwhile  topic  for  research. 

Quantum  networks  are  a  tantalizing  concept.  At  its  heart,  a  quantum  network 
consists  of  a  number  of  controllable  quantum  systems  which  can  interact  with  each 
other  while  separated  by  some  distance  [2],  It  therefore  requires  both  a  stationary 
quantum  node  (where  information  is  stored  and  processed)  and  an  appropriate  in¬ 
terface  by  which  distant  nodes  can  communicate  over  a  quantum  channel.  Most 
research  on  quantum  networks  has  focused  on  optically  controlled  quantum  nodes 
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that  interact  via  the  emission  and  absorption  of  photons.  Recently,  a  German  team 
demonstrated  an  elementary  two-node  network  whereby  distributed  entanglement  was 
created  across  a  distance  60  m  as  a  single  8'Rb  atom  emitted  a  single  photon  which 
was  transmitted  down  an  optical  fiber  and  captured  by  a  second  trapped  87Rb  atom 
[3], 

There  exist  several  promising  candidate  systems  for  quantum  network  nodes: 
namely  the  trapped  neutral  atoms  mentioned  already,  nitrogen-vacancy  color  cen¬ 
ters  in  diamond  [4],  and  III-V  semiconductor  quantum  dots  [5].  However,  all  these 
candidate  systems  have  optical  transitions  in  the  visible  or  near  visible  spectral  re¬ 
gion.  If  one  is  interested  in  extending  the  distance  over  which  distributed  entangled 
states  can  be  created,  it  is  evidently  necessary  to  pursue  a  system  whereby  the  pho¬ 
tonic  channel  consists  of  transmission  of  single  photons  in  the  1 . 5- /irn  spectral  region, 
where  the  transmission  loss  in  optical  fiber  is  at  its  lowest  value,  below  0.2  dB/km. 
Such  links  require  either  development  of  a  quantum  node  with  1.5-/im  optical  tran¬ 
sitions,  or  use  of  a  frequency  conversion  technique  in  which  the  carrier  frequency  of 
light  can  be  translated  while  maintaining  the  quantum  characteristics  of  the  signal. 
This  dissertation  is  concerned  with  the  latter,  which  is  accomplished  through  a  tech¬ 
nique  called  quantum  frequency  conversion,  first  proposed  by  Kumar  in  1990  [6],  and 
demonstrated  shortly  afterwards  [7]. 

While  we  have  here  framed  the  discussion  of  quantum  frequency  conversion  in 
the  context  of  quantum  networks,  it  has  other  applications  in  various  scientific  and 
technological  fields.  Frequency  conversion  in  x2'1  media  has  been  a  subject  of  active 
research  since  first  observation  of  x ^  processes  by  Franken  just  over  50  years  ago 
[8].  It  has  been  used  primarily  for  the  generation  of  radiation  at  optical  frequencies 
inaccessible  with  current  laser  technology,  but  in  recent  years  the  development  of 
highly  efficient  x(2>  frequency  converters  based  on  waveguides  has  enabled  its  use 
in  nonlinear  all-optical  signal  processing  [9].  The  results  of  this  dissertation  can  be 
viewed  as  an  extension  of  classical  signal  processing  applications  to  the  quantum 
realm.  As  quantum  communications  systems  are  in  their  infancy,  we  believe  that  the 
potential  for  future  research  along  this  direction  is  promising. 
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1.1  Overview  of  this  dissertation 

Chapter  2  presents  the  theory  of  frequency  conversion  in  materials  using  both 
classical  and  quantum  mechanical  descriptions.  We  give  an  overview  of  experimental 
achievements  in  this  field,  and  point  out  the  primary  impediments  to  progress  to¬ 
wards  a  long-distance  quantum  network.  Chapter  3  gives  a  description  of  the  device 
technology  we  use  for  QFC:  reverse-proton-exchange  waveguides  in  periodically  poled 
lithium  niobate.  The  principles  of  quasi-phasematched  frequency  conversion  devices 
are  presented,  along  with  device  fabrication  by  periodic  poling  of  ferroelectrics.  We 
then  describe  the  waveguide  technology,  giving  some  special  attention  to  the  fabrica¬ 
tion  and  characterization  procedures  used  in  our  lab  to  produce  the  highest  quality 
devices  possible. 

A  primary  contribution  of  this  dissertation  is  a  detailed  discussion  of  the  non¬ 
idealities  of  the  QFC  process  via  the  generation  of  noise  photons,  which  will  be 
presented  in  Chapter  4.  Apart  from  technical  noise  (e.g.  stray  light  or  improperly 
filtered  pump  sources),  the  sources  of  noise  in  QFC  devices  can  be  traced  to  two 
physical  effects.  The  first  is  parametric  fluorescence,  a  process  endemic  to  non¬ 
linear  materials.  Naive  calculations  of  parametric  fluorescence  show  that  in  the  highly 
phase-mismatched  regime  it  should  not  be  a  significant  noise  source,  but  imperfec¬ 
tions  in  the  fabrication  of  periodically  poled  gratings  enhance  the  generation  rate 
substantially,  by  up  to  a  factor  of  about  104.  Tolerances  for  these  fabrication  errors 
are  derived.  The  second  noise  source  is  spontaneous  Raman  scattering.  By  a  careful 
investigation  of  the  Raman  spectrum  of  LiNbOs  we  find  certain  system  configurations 
by  which  Raman  scattering  can  be  avoided,  leading  to  design  rules  by  which  we  can 
achieve  low-noise  quantum  frequency  conversion  for  diverse  applications  in  science 
and  technology. 

Chapters  5  and  6  present  applications  of  QFC  technology  to  the  up-  and  down- 
conversion  of  single  photons.  Upconversion  has  been  used  as  a  technique  to  improve 
detection  of  infrared  radiation  since  the  early  days  of  nonlinear  optics  [10].  It  has  also 
been  the  primary  technological  driver  of  quantum  frequency  conversion:  quantum  key 
distribution  and  other  quantum  communications  technologies  require  high-efficiency, 


CHAPTER  1.  INTRODUCTION 


4 


fast,  and  low-noise  single-photon  detectors  for  light  in  the  spectral  range  in  which 
optical  fibers  (the  primary  quantum  channels)  have  the  lowest  attenuation:  namely 
the  region  around  A  =  1.5  fim.  While  some  emerging  1 . 5- /irn  single-photon  detection 
technologies  are  improving,  the  best  single-photon  detectors  are  those  based  on  Si, 
which,  due  to  the  1.1-eV  bandgap  of  Si,  are  unable  to  detect  photons  with  wavelengths 
longer  than  approximately  1  /irn.  Chapter  5  describes  our  efforts  at  developing  up- 
conversion  single-photon  detectors:  we  achieved,  to  the  best  of  our  knowledge,  the 
highest  efficiency  long-wavelength-pumped  (low-noise)  upconversion  detector  for  the 
1.5-/im  band,  and  demonstrated  a  l-/im  upconversion  detector,  with  approximately 
30  noise  counts/s  at  peak  conversion  efficiency. 

Chapter  6  presents  a  downconversion  interface  suitable  for  linking  an  InAs/GaAs 
quantum  dot  with  an  optical  transition  at  910-nm,  a  promising  quantum  network 
node  candidate,  with  the  1.5-/im  spectral  window  suitable  for  long-distance  transmis¬ 
sion  down  optical  fiber.  We  demonstrate  time-resolved  downconversion  of  the  single 
photons  emitted  from  the  quantum  dot,  and  show  that  the  photon  statistics  of  the 
quantum  dot  emission  are  preserved.  Additionally,  through  the  observation  of  Rabi 
oscillations  in  the  converted  photons,  we  prove  that  downconversion  system  is  fully 
compatible  with  the  types  of  coherent  control  operations  on  the  quantum  dot  spin 
needed  for  its  use  as  a  quantum  network  node. 

In  Chapter  7,  we  move  past  simple  up-  and  downconversion  devices  and  explore 
additional  design  space  to  bring  added  functionality  to  single-photon  frequency  con¬ 
version  systems.  We  present  two  main  concepts:  first,  we  describe  the  demonstration 
of  cascaded  frequency  conversion.  As  it  is  possible  to  achieve  near-unity  conversion 
efficiency  of  a  single  process  in  RPE  PPLN  waveguides,  using  approximately  twice 
as  much  pump  power  it  is  possible  to  convert  an  input  signal  twice.  As  a  proof-of- 
concept,  we  demonstrate  two-stage  frequency  conversion  of  a  1.5-/im  single-photon 
signal  to  the  green  spectral  region  for  high-speed  single  photon  detection.  The  second 
device  concept  is  based  on  a  multi-channel  QFC  device,  whereby  one  can  switch  on 
and  off  multiple  pumps  to  demultiplex  a  single-photon  signal  to  enable  higher  clock- 
rate  quantum  communications  systems.  In  Chapter  8,  we  conclude  the  dissertation, 
and  describe  some  avenues  for  further  research  and  development. 


Chapter  2 

Theory  and  review  of  quantum 
frequency  conversion 


In  the  following  sections  we  provide  a  theoretical  basis  for  both  classical  and  quantum 
mechanical  descriptions  of  three-wave  interactions  in  y ^  nonlinear  waveguides,  and 
describe  how  these  interactions  can  be  used  to  translate  the  carrier  frequency  of  a 
quantum  state  of  light.  We  will  also  briefly  review  experimental  work  in  this  field. 

2.1  Three- wave  mixing  in  waveguides 

Classical  electromagnetic  field  theory  can  be  used  to  describe  nonlinear  optical  in¬ 
teractions  [11].  Classical  electromagnetism,  encapsulated  in  Maxwell’s  equations  and 
Newton’s  laws,  describes  the  motion  of  free  and  bound  charges  in  response  to  elec¬ 
tric  and  magnetic  fields,  and  then  how  the  motion  of  charges  acts  as  a  source  of 
electromagnetic  radiation.  As  a  prototypical  three-wave  interaction,  we  will  study 
sum- frequency  generation  (SFG),  in  which  a  signal  at  uq  is  combined  with  a  pump 
at  ojp  to  produce  an  output  at  u>2  =  oq  +  c op.  The  goal  of  this  treatment  will  be  to 
derive  coupled  wave  equations  for  field  amplitudes  at  each  of  the  frequencies  ujj  with 
j  e  {1,  p,  2}.  This  subject  is  treated  in  all  standard  texts  on  nonlinear  optics  [12], 
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so  we  offer  only  an  abbreviated  discussion.  Maxwell’s  equations  are  given  by 


„  -  <9B 

VxE  =  “9V 

(2.1a) 

-  dD 

VXH=J+-, 

(2.1b) 

V  •  D  =  p, 

(2.1c) 

V-B  =  0, 

(2.  Id) 

where  E  and  H  are  the  electric  and  magnetic  fields,  D  is  the  electric  displacement 
held,  B  is  the  magnetic  flux  density,  and  J  and  p  are  the  free  current  and  charge 
densities  [13].  We  use  tildes  to  denote  quantities  that  are  rapidly  varying  in  time  at  an 
assumed  optical  carrier  frequency,  and  boldface  to  denote  vectors.  The  relationships 
between  E  and  D  and  between  H  and  B  are  given  by  the  constitutive  relations, 
where  we  assume  a  nonmagnetic  material  and  separate  the  polarization  response  of 
the  medium  into  linear  and  nonlinear  components  Pl  and  Pnl-  We  incorporate  the 
linear  component  into  the  electric  held  via  the  use  of  a  relative  permittivity  er: 


B  =  /i0H,  (2.2a) 

D  =  e0erE  +  Pnl;  (2.2b) 


where  er  =  1  +  yd)  =  n2  where  yd)  js  the  linear  susceptibility  of  the  medium  an  n  is 
the  refractive  index.  Taking  the  curl  of  Eq.  (2.1a)  and  using  Eq.  (2.1b)  we  obtain 


<9E 


<92E 


V  E  +  +  c2  gt2 


1  <92P 


NL 


£0c 


dt 2 


(2.3) 


where  we  have  assumed  a  source- free  medium  (i.e.  p  =  0),  and  used  the  vector 
identity 

V  x  (V  x  e)  =  V  (V  •  E)  -  V2E.  (2.4) 

Neglecting  the  the  V(V  •  E)  component  [12],  Eq.  2.3  is  evidently  a  driven,  damped 
vector  wave  equation.  Note  that  we  have  introduced  an  effective  conductivity  a  such 
that  J  =  crE.  which  incorporates  the  effect  of  losses. 
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We  will  be  concerned  primarily  with  nonlinear  optical  interactions  in  waveguides, 
which  consist  of  discrete  bound  modes  and  a  continuum  of  radiation  modes  (which 
can  usually  be  neglected).  We  first  give  a  useful  factorization  of  the  full  electric  field 
of  a  particular  frequency  component  at  ujj,  where  the  normalizations  serve  to  provide 
a  bridge  to  the  quantum  mechanical  description  in  Sec.  2.2.  Since  the  waveguides 
used  in  this  dissertation  are  weakly  confining  and  support  only  a  single  polarization 
of  light  [14],  we  now  drop  the  vector  nature  of  the  electric  fields: 

Ej(x,y,z,t )  =  ijuj{x)y)aj{z)U{ujt~kjZ\  (2.5) 

Here,  a,j(z)  is  a  slowly  varying  envelope  with  units  [|aj(z)|2]  =  photons/s,  kj  = 
oJjUj /c  is  the  propagation  constant  (where  rij  is  taken  to  be  the  effective  index  of  the 
waveguide  mode),  and  Uj(x,  y)  is  the  spatial  mode  of  the  waveguide  and  is  normalized 
such  that 

J j  \uj(x,v)\2  dx  d y  =  1.  (2.6) 

The  scale  factor  =  (2 Zjhujj)1^2,  where  Z3  =  Z0/nj  where  Z0  =  \] /r0 / e0  ~  376.73  O 
is  the  impedance  of  free  space,  and  h  =  h/2n  is  Planck’s  constant. 

To  arrive  at  coupled  wave  equations,  it  is  necessary  to  specify  the  form  of  the 
frequency  components  of  the  nonlinear  polarization  Pnl-  For  sum-frequency  genera¬ 
tion,  we  are  interested  in  the  upconversion  of  a  signal  at  uj\  via  the  interaction  with 
a  pump  at  ujp,  producing  a  target  field  at  oj2  =  cui  +cup.  The  scalar  magnitudes  of  the 
nonlinear  polarization  at  these  three  frequencies  are 


PnlO^)  =  2eodE(uji)E(ujp) 

(2.7a) 

Pnl(wi)  =  2eodE{oj2)E*{ojp) 

(2.7b) 

PnlO^p)  =  2eodE(u}2)E*(uji) 

(2.7c) 

where  d  =  y^/2  is  the  nonlinear  coefficient  and  is  a  function  of  position:  it  has 
lateral  dependence  due  to  the  waveguide  structure  and  axial  dependence  due  to  quasi¬ 
phasematching  (QPM,  see  Sec.  3.1). 

Coupled  wave  equations  can  be  derived  by  inserting  the  nonlinear  polarization  of 
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Eqs.  (2.7)  into  Eq.  (2.3).  Inserting  the  complete  electric  field  (consisting  of  a  sum 
of  the  waves  at  the  three  frequencies  involved  in  the  problem)  and  using  the  slowly 
varying  envelope  approximation  [12],  one  obtains  the  following  coupled  first-order 
nonlinear  differential  equations: 


(cF  +  ir)  ai  =  ~iKapa‘2e~lAk'z 

(cF  +  ~2^)  ap  =  — *Kaia2  e~lAk'z 

(f  +  y)  0,2  =  ~iKaiaPe%Ak'zi 


(2.8a) 

(2.8b) 

(2.8c) 


where  A k'  =  k2  —  kp  —  k\  is  the  phase  mismatch,  oij  is  the  power  loss  coefficient 
(j  G  {1,  p,  2}),  and  the  coupling  coefficient 

K  =  tod(2N^A)1/2e.  (2.9) 

v  nrnpn2  J 

Here,  rij  is  the  effective  index  of  the  waveguide  mode  at  frequency  coj,  and  0  is  the 
mode-overlap  integral  in  the  waveguide,  given  by 


0 


d(x,y)u,(x,y)up(x,y)u2(x,y)  di  dy, 


(2.10) 


where  we  have  incorporated  the  spatial  form  of  the  nonlinear  coefficient  d  =  dod(x,  y), 
where  max  { |  d(x:  y)  |  }  =  1. 

As  an  illustrative  example,  we  shall  proceed  to  solve  Eqs.  (2.8)  for  the  case 
a\(z  =  0)  =  aio  and  a2(z  =  0)  =  0,  corresponding  to  the  upconversion  of  a  weak 
signal.  We  make  several  simplifying  asumptions.  First,  the  input  signal  aio  is  much 
weaker  than  the  pump  ap,  such  that  we  assert  ap(z)  =  ap( 0),  i.e.  an  undepleted 
pump,  and  thus  can  set  the  right-hand  side  of  Eq.  (2.8b)  equal  to  0.  We  also  consider 
perfect  phasematching,  that  is,  that  A k'  ^  0  has  been  corrected  by  some  phase- 
matching  technique  (e.g.  QPM,  see  Sec.  3.1),  and  also  that  the  propagation  losses 
of  the  waveguide  are  negligible.  This  last  assertion  is  not  strictly  true  for  the  wave¬ 
guide  technology  described  in  Chapter  3,  but  the  low  losses  of  our  waveguides  do  not 
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Figure  2.1:  Photon  flux  evolutions  ai(z)  and  CI2O)  i11  a  classical  SFG  process  where 
ai(0)  =  aio,  and  02(0)  =  0.  Complete  conversion  occurs  after  a  distance  z  =  71/27. 

substantially  change  the  conversion  dynamics  but  rather  only  reduce  the  maximum 
conversion  efficiency.  Rewriting  Eqs.  (2.8)  with  these  assumptions,  we  have 


dcq 

d^  =  “7fl2 

(2.11a) 

da2 

d^=7ai 

(2.11b) 

where  7  =  nAp  and  we  write  ap  =  Apel^p  (with  Ap  real  and  positive)  and  have  chosen 
the  arbitrary  overall  phase  (j)p  =  7t/2.  The  solution  to  Eqs.  (2.11)  proceeds  straight¬ 
forwardly  by  differentiating  Eq.  (2.11b)  and  inserting  the  equality  of  Eq.  (2.11a), 
whereby  we  are  led  to  the  solutions 


a±(z)  =  aiocosyz,  (2.12a) 

a-2(z)  =  aio  sin  72:.  (2.12b) 


These  solutions  are  plotted  in  Fig.  2.1.  We  see  that  complete  conversion  of  a  signal 
at  (jj\  to  u>2  occurs  after  a  propagation  distance  of  z  —  73/27. 
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2.2  Quantum  frequency  conversion 

Eqs.  (2.12)  give  the  classical  solution  to  the  upconversion  of  a  signal  at  by  in¬ 
teraction  with  a  pump  at  c op  via  SFG.  One  attains  a  proper  quantum  mechanical 
treatment  of  three-wave  interactions  via  the  Hamiltonian  of  the  quantized  electric 
holds  in  a  nonlinear  dielectric  [15,  16].  This  is  a  subtle  process,  and  is  beyond  the 
scope  of  the  dissertation.  However,  it  can  be  shown,  and  all  experimental  observations 
have  verified,  that  Eqs.  (2.11)  give  the  proper  quantum  mechanical  description  if  the 
held  envelopes  cq  and  a2  are  replaced  with  the  annihilation  operators  di  and  d2  of 
photons  at  frequencies  uq  and  cu2.  The  operator  di  (dj)  annihilates  (creates)  one  pho¬ 
ton  at  frequency  u\.  These  operators  obey  the  commutation  relations  [dj,  dj,]  =  Sjk, 
where  djk  is  the  Kronecker  delta,  and  have  the  properties 

d|n)  =  \/n\n  —  1),  (2.13a) 

d)\n)  =  y/n  +  l|n  +  1),  (2.13b) 

where  |n)  is  the  Fock  (photon-number)  state  containing  n  photons,  and  a|0)  =  0  [17]. 

Here,  we  present  the  input-output  relation  (after  propagation  over  a  distance 
L  corresponding  to  the  length  of  the  device)  for  the  annihilation  operators  d\  and 
d2  as  the  two  helds  exchange  photons  through  the  conversion  process.  A  quant um- 
mechanically  accurate  version,  unlike  the  results  of  Eqs.  (2.12),  must  include  both 
input  operators  such  that  the  effects  of  vacuum  fluctuations  are  included.  Working  in 
the  Heisenberg  picture  (whereby  the  operators  evolve  in  time  but  the  quantum  state 
does  not),  this  input-output  relation  is  given  by 

di(L)  =  d^O)  cos 7-L  +  d2(0)  sin yL,  (2.14a) 

d2(L)  =  —  di(0)  sinyL  +  d2(0)  cosyL.  (2.14b) 

We  see,  evidently,  that  if  yL  =  tt/2,  up  to  an  arbitrary  overall  phase  factor  the 
annihilation  operators  have  switched  places  [6] 


d]_(0)  =>  — d 2(L)  and  d2(0)  di (L). 


(2.15) 
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Figure  2.2:  Schematic  diagram  of  the  quantum  frequency  conversion  process  in  a  x ^ 
material.  A  quantum  state  of  light  at  uq  (012)  is  upconverted  (downconverted)  to  oj2 
(uq)  via  interaction  with  a  pump  field  at  ojp,  where  oq  +  uip  —  u>2. 


Eq.  (2.15)  captures  the  concept  of  the  quantum  frequency  converter,  which  is  shown 
schematically  in  Fig.  2.2.  Through  interaction  with  a  strong  pump,  if  the  device 
length  and  pump  power  are  chosen  such  that  7 L  =  7t/2,  the  annihilation  operators  for 
photons  at  the  input  and  output  have  swapped  frequencies.  Importantly,  this  means 
that  the  quantum  frequency  converter  is  capable  of  simultaneously  upconverting  and 
downconverting  two  input  quantum  states.  I11  most  experimental  implementations, 
one  of  the  input  states  will  be  a  vacuum  state.  However,  based  on  the  similarity  of 
Eqs.  (2.14)  to  a  beamsplitter  input-output  relation  [18],  there  is  a  proposal  to  use  the 
quantum  frequency  converter  as  a  method  by  which  one  could  interfere  two  photons 
of  different  color  [19]. 

We  consider  one  example  of  the  effects  of  quantum  frequency  conversion  on  an 
input  quantum  state,  and  calculate  the  conversion  efficiency  of  a  signal  at  07  to  a 
target  at  uj2.  We  consider  an  input  state  j^)  =  |ni)Wl |0)W2  =  |niO).  We  calculate 
the  conversion  efficiency  as  the  ratio  of  photon  number  at  the  target  frequency  at  the 
output  of  the  device  to  the  photon  number  at  the  signal  frequency  at  the  input  of  the 
device: 

_  (jj\h2{LM)  =  l^\a\{L)a2{L)\^) 

(V>|hi(0)|^)  (^|d|(0)di(0)|^) 

=  sin27L  (2.16b) 


where  the  result  of  Eq.  (2.16b)  is  due  to  applying  Eqs.  (2.14)  and  (2.13)  to  Eq.  (2.16a). 
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Figure  2.3:  Conversion  efficiency  rj  versus  normalized  pump  power  Pp/Pmax  for  a  SFG 
process,  where  Pmax  =  7r2/(4L2?/nor),  as  given  in  Eq.  (2.19).  Note  that  this  result  is 
the  same  as  is  obtained  for  the  classical  case  of  Eqs.  (2.12). 

These  results  are  also  evident  classically  from  Eqs.  (2.12).  Generalizing  further,  we 

note  that  for  a  general  state  | ^12) ,  quantum  properties  of  the  two  modes  are  de¬ 

termined  by  (A)  =  <^i2|/i(ai)hM,  and  (A)  =  <-012 1 A (^2) |'0i2>  where  A  (A)  is  a 
generic  function  of  the  annihilation  operator  Si  (a 2).  Therefore,  for  7 L  —  7t/2,  when 
Eq.  (2.15)  holds,  we  have 

(A  A  =  L ))  =  (lMA[fi2(0)]hM  (2.17a) 

(A  (2  =  £))  =  <^12 1 A  [— ai  (0)]  1-012)  -  (2.17b) 

We  see,  therefore,  that  the  quantum  properties  of  the  mode  Si  (S2)  at  the  device 
output  are  determined  by  the  input  mode  S2  (Si)  for  a  complete  quantum  frequency 
conversion  process  [6]. 

It  is  worthwhile  to  express  the  conversion  efficiency  of  Eq.  (2.16)  and  the  coupling 
coefficient  7  in  terms  of  experimental  parameters  and  physical  constants.  Recalling 
that  7  =  kAp,  where  k  is  given  by  Eq.  (2.9)  and  using  the  held  definitions  of  Eq.  (2.5), 
we  define  a  normalized  efficiency  of  the  device  ?/nor  by  equating  7  =  \JrjnovPp,  where 
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Pp  is  the  pump  power.  Using  this  definition,  we  find  that 


^nor  =  e0^ol®|2 


2u}\U}2Zq 

nin2n3 


(2.18) 


As  the  main  experimental  parameter  once  a  device  is  fabricated  is  the  pump  power, 
we  plot  the  conversion  efficiency  rj  as  a  function  of  Pp, 


Tj  =  Sill 


VnoiPpL )  =  sin2  I  ^  '  Pp 


2  VP 

Zj  \  1  rr 


(2.19) 


in  Fig.  2.3,  where  we  have  defined  Pmax  as  the  pump  power  requried  for  complete 
conversion,  given  by  Pmax  =  7r2/(4L2/7nor).  For  pump  powers  larger  than  Pmax,  back- 
conversion  reduces  the  conversion  efficiency. 


2.3  Historical  overview 

Quantum  frequency  was  first  proposed  by  Prem  Kumar  in  1990,  who  showed  theoret¬ 
ically  that  the  properties  of  a  quantum  state  of  light  are  unchanged  upon  frequency 
conversion  in  a  material  [6].  The  first  experimental  demonstration  of  QFC  fol¬ 
lowed  in  1992,  where  Kumar’s  group  was  able  to  upconvert  a  squeezed  state  from  1064 
nm  to  532  nm,  and  it  was  proposed  that  QFC  can  provide  a  wavelength-tunable  source 
of  nonclassical  light  [7].  In  this  initial  demonstration,  the  nonlinear  medium  was  a 
bulk  KTP  crystal;  due  to  the  low  nonlinearity  it  was  necessary  to  use  a  Q-switched 
mode-locked  Nd:YAG  laser  to  achieve  appreciable  conversion  efficiencies. 

Following  this  demonstration,  it  seems  that  QFC  was  not  actively  pursued  for 
several  years.  A  resurgence  of  interest  followed  the  development  of  entangled  photon 
sources  based  on  spontaneous  parametric  downconversion  (SPDC)  in  y ^  crystals, 
which  provided  a  much  more  convenient  source  of  nonclassical  light  than  the  atomic 
physics  techniques  which  had  been  used  before  that  time  [20].  Additionally,  follow¬ 
ing  the  development  of  experimental  quantum  cryptography  [21],  it  was  realized  that 
the  limiting  factor  for  system  performance  was  the  relatively  poor  single-photon  de¬ 
tectors  for  telecom-band  photons  [22],  Vandevender  and  Kwiat  developed  a  telecom 
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single-photon  detector  based  on  quantum  frequency  upconversion  of  the  telecom-band 
photons,  followed  by  detection  with  a  higher-performance  Si  APD  single-photon  de¬ 
tector  [23].  Again,  short  pump  pulses  (a  few  ns,  at  Xp  =  1064  nm)  were  needed  due 
the  use  of  a  bulk  crystal,  but  here  the  use  of  PPLN  enabled  80%  conversion  efficiency 
although  at  a  peak  power  level  of  12  kW. 

We  note,  however,  that,  dating  to  the  earliest  days  of  nonlinear  optics,  upcon¬ 
version  was  proposed  as  a  technique  to  improve  infrared  detector  technology  [11]. 
In  the  late  1960s,  Midwinter  demonstrated  upconversion  of  both  optical  signals  [10] 
and  images  [24]  at  1.7  /im  to  the  visible  using  a  ruby  laser  as  a  pump.  While  it  was 
realized  that  upconversion  occurs  on  a  photon-by-photon  basis,  the  implications  for 
quantum  optics  were  not  yet  realized:  a  reasonably  comprehensive  theory  of  quan¬ 
tum  optics  had  only  recently  been  developed  by  Glauber  [25].  In  the  first  experiment, 
conversion  efficiencies  of  approximately  1%  were  achieved,  using  a  LiNbCU  crystal  us¬ 
ing  noncritical  birefringent  phasematching  [12],  Upconversion  techniques  were  later 
extended  further  into  the  infrared  for  use  in  astronomy  [26].  In  both  of  the  image 
upconversion  experiments  mentioned  here,  the  conversion  efficiences  were  limited  to 
approximately  10-7.  We  note  that  there  has  been  a  recent  theoretical  proposal  on 
the  frequency  upconversion  of  a  quantum  image,  a  state  of  light  in  which  there  are 
quantum  correlations  between  the  various  spatial  modes  comprising  the  state  [27]. 

The  first  demonstrations  of  continuous-wave  upconversion  for  single-photon  de¬ 
tection  were  based  on  a  RPE  PPLN  waveguide  [28] ,  or  a  bulk  PPLN  crystal  placed  in 
a  pump-resonant  cavity  [29].  In  the  resonant  cavity  of  Albota  and  Wong,  up  to  90% 
conversion  efficiency  was  observed  for  a  circulating  pump  power  of  21.6  W.  Roussev 
and  coworkers  observed  a  signal  depletion  level  of  99%,  with  overall  conversion  ef¬ 
ficiency  limited  to  82%  due  to  waveguide  losses  of  approximately  0.18  dB/cm.  The 
waveguide  used  was  of  the  RPE  type  that  will  be  described  in  the  following  chapter, 
and  had  Pm ax  =  88  mW,  approximately  350  times  lower  than  for  the  bulk  PPLN 
device  of  [29]. 

A  common  issue  in  the  first  generation  of  quantum  frequency  upconverters  was 
a  very  high  level  of  noise  photons.  As  an  extremely  large  flux  of  pump  photons  is 
required  to  obtain  efficient  conversion,  even  very  weak  spontaneous  inelastic  scattering 
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processes  can  significantly  reduce  the  system  performance.  The  two  spontaneous 
scattering  processes  responsible  for  noise  in  QFC  devices  are  SPDC  [23,  29]  and 
spontaneous  Raman  scattering  [30].  In  fact,  even  in  the  early  upconversion  papers  it 
was  found  that  noise  processes  were  a  significant  limiter  of  the  system  performance 
[10].  In  Chapter  4  we  present  a  detailed  discussion  of  these  noise  processes,  which 
had  not  been  studied  in  detail  until  this  work.  However,  it  was  realized  that  the 
use  of  a  pump  frequency  u>p  <  UJ\  is  highly  effective  at  reducing  the  noise,  as  it 
becomes  difficult  for  photons  scattered  from  the  strong  pump  to  overlap  spectrally 
with  the  acceptance  bandwidth  of  the  converter  near  04 .  Two  low-noise  upconversion 
detectors  demonstrated  this  fact,  one  for  the  l-/im  band  [31],  and  a  second  for  the 
1.55-/im  telecom  band  [32], 

To  date,  upconversion  has  been  used  to  bring  significant  added  functionality  to 
quantum  optics  and  quantum  information  systems.  Low-jitter  upconversion  detectors 
were  developed,  which  were  useful  in  QKD  systems  [33,  34].  In  a  pioneering  exper¬ 
iment,  Tanzilli  and  collaborators  showed  that  if  the  signal  photon  of  an  entangled 
photon  pair  was  upconverted,  its  entanglement  (using  the  time-bin  variable)  with  the 
idler  photon  was  preserved  [35].  Recently,  polarization  entanglement  was  also  shown 
to  be  preserved  [36].  Two  additional  concepts  of  technological  importance  were  the 
use  of  the  combination  of  the  narrow  acceptance  bandwidth  of  the  upconverter  with 
a  tunable  pump  frequency  to  enable  a  single-pixel  upconversion  spectrometer  [37], 
and  the  demonstration  of  high-dynamic-range  time-domain  reflectometry  for  charac¬ 
terizing  fiber  optic  links  [38].  Recently,  interfaces  of  quantum  frequency  upconverters 
with  a  1.3-/im  quantum  dot  single-photon  source  were  demonstrated,  showing  the 
preservation  of  photon  antibunching  upon  upconversion  [39],  and  temporal  shaping 
by  modulation  of  the  pump  waveform  [40].  Additionally,  it  has  been  shown  the  phase 
coherence  of  a  quantum  state  is  preserved  [41]  and  that  photon  statistics  (up  to  fourth 
order)  are  preserved  upon  upconversion  [42], 

While  the  bulk  of  experimental  work  on  QFC  has  focused  on  upconversion,  it  was 
realized  early  on  that  downconversion  of  quantum  states  of  light  is  also  possible.  As 
was  discussed  in  Chapter  1,  downconversion  represents  a  critical  functionality  needed 
to  connect  quantum  nodes  with  transitions  in  the  visible  or  near-visible  spectral 
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region  with  the  telecom  band,  allowing  long-distance  quantum  communication  over 
optical  fiber.  This  application  is  in  contrast  to  a  naive  view  of  “downconversion”  as 
a  process  that  involves  the  amplification  of  quantum  noise  [43].  Rather,  in  a  single¬ 
photon  downconverter,  the  single-photon  is  at  the  highest  frequency  and  the  strong 
pump  driving  the  conversion  is  at  a  lower  frequency1,  such  that  the  resulting  OPA 
has  negligible  gain,  and  therefore  negligible  noise  (but  of  course  subject  to  the  same 
constraints  of  parasitic  spontaneous  scattering  processes  as  apply  in  upconversion 
devices).  A  theoretical  discussion  of  single-photon  downconversion  was  presented  by 
Z.  Y.  Ou  in  2008  [44],  where  it  was  shown  that  QFC  enables  frequency  translation 
of  quantum  states  and  the  storage  of  photonic  states  in  stationary  matter-based 
quantum  memories.  A  number  of  single-photon  downconversion  experiments  followed 
in  the  next  three  years  [45,  46,  47,  48,  49].  These  experiments  have  shown  preserved 
coherence  in  interference  experiments,  and  have  also  demonstrated  high  conversion 
efficiencies.  In  a  recent  experiment,  entanglement  properties  have  also  been  shown  to 
be  preserved  upon  downconversion,  as  was  done  by  Tanzilli  in  upconversion  [50].  This 
dissertation  presents  the  first  demonstration  of  a  downconversion  interface  to  a  solid- 
state  qubit  system,  which  is  expected  to  be  of  interest  for  quantum  communication 
based  on  quantum  repeaters  [51]. 

While  this  dissertation  focuses  on  quantum  frequency  conversion  in  media, 
there  have  been  recent  proposals  and  a  few  experimental  demonstrations  using  four- 
wave  mixing  in  yF)  media.  In  one  such  implementation,  known  as  Bragg  scattering, 
two  pumps  at  copi  and  ujp2  are  used  to  convert  a  signal  at  ny  to  targets  at  u>2±  = 
uji  ±  ( ujpi  —  ujP2 ),  although  phasematching  can  be  used  to  preferentially  create  U2+  or 
0J2-  [52],  Since  one  signal  photon  and  one  pump  photon  are  destroyed  as  a  target 
photon  and  a  second  pump  photon  are  created,  this  process,  similar  to  SFG  and 
single-photon  DFG,  should  also  be  noiseless.  The  first  experimental  demonstration  of 
Bragg  scattering  for  a  weak  signal  was  done  in  2006,  and  a  signal  depletion  of  22  dB 
was  observed,  with  a  large  amount  of  noise  (on  single-photon  scales,  but  moderate  on 
classical  telecommunications  scales)  due  to  Raman  scattering  as  all  signals  were  in  the 

Rn  this  dissertation,  we  always  use  the  term  “pump”  to  refer  to  the  field  driving  the  conversion 
process,  in  contrast  with  optical  parametric  amplification  (OPA)  terminology. 
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1.5-/im  telecom  band  [53].  Demonstration  of  QFC  with  x ®  processes  was  achieved 
using  a  photonic  crystal  hber,  with  depletion  efficiencies  of  29%  [54] .  The  use  of  two 
pumps  makes  the  yf3)  QFC  techniques  attractive  for  frequency  translation  between 
nearby  frequencies  where  a  x  2  technique  would  require  extremely  low  ojp.  However, 
it  should  be  noted  that  implementations  of  yy3)  -based  QFC  techniques  have  suffered 
a  great  deal  from  noise  photons  primarily  due  to  spontaneous  Raman  scattering. 


Chapter  3 

Re  ver  se-  prot  on-  exchange 
waveguides  in  periodically  poled 
lithium  niobate 


The  results  presented  in  this  dissertation  rely  heavily  on  the  performing  optical  fre¬ 
quency  conversion  with  high  efficiency.  Optical  nonlinearities  are  famously  extremely 
weak  [11],  and  it  is  only  recently  that  researchers  have  been  able  to  attain  the  con¬ 
version  efficiencies  necessary  for  the  realistic  implementation  of  frequency  conversion 
at  the  single-  or  few-photon  level  [23,  28,  29]. 

The  device  platform  used  in  this  work  is  based  on  annealed  and  reverse  proton 
exchange  (RPE)  waveguides  in  periodically  poled  lithium  niobate  (PPLN)  [55].  The 
development  of  sophisticated  materials  and  process  models  by  previous  generations 
of  Fejer-group  students  has  enabled  a  reliable  technique  where  devices  can  be  de¬ 
signed  for  any  application  of  interest  [56] .  In  the  next  section,  we  describe  the  theory 
and  practice  of  quasi-phasematching  by  periodic  poling.  In  Sec.  3.2,  we  discuss  the 
fabrication  process  used  for  making  RPE  waveguides,  and  discuss  the  “toolbox”  of 
waveguide  components  that  enables  the  applications  discussed  in  following  chapters. 
In  Sec.  3.3  we  outline  the  characterization  procedures  by  which  fabricated  devices  are 
screened  for  performance,  and  describe  the  fiber-pigtailing  technique  by  which  the 
devices  can  be  integrated  into  fiber-based  quantum  optical  systems. 
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3.1  Quasi-phasematching  and  periodic  poling 

For  the  work  in  this  dissertation  we  are  primarily  concerned  with  high-efficiency 
sum-  or  difference-frequency  generation  processes.  Eqs.  (2.8)  gave  the  coupled  wave 
equations  for  an  SFG  process  in  which  a  target  held  at  u 2  was  generated  from  a  signal 
cu i  and  pump  cup,  where  u2  =  u i  +  top-  Neglecting  propagation  losses,  the  growth  of 
the  held  at  the  sum-frequency  (SF)  is  given  by 

^  =  rd0a1ape<Ate,  (3.1) 

az 

where  the  T  =  in/ d0,  with  n  given  by  Eq.  (2.9),  and  the  phase  mismatch  A k'  = 
k2  —  k\  —  kp,  where  kj  =  u> jn{uj j) / c  is  the  wave  vector  at  frequency  j  G  {p,  1,  2}. 
Due  to  dispersion,  the  quantity  Ak  will  in  general  be  nonzero.  In  the  limit  of  low 
conversion  where  a\(z )  =  ai(0)  and  ap(z )  =  ap(0),  the  evolution  of  a2(z)  will  thus  be 
governed  by  the  rotation  of  the  complex  exponential  exp(iAk'z),  causing  oscillations 
in  the  power  P2(z)  =  hco2\a2(z)\2 ,  with  periodic  reversals  in  power  how  occurring  every 
coherence  length  Lc  =  ir/Ak'.  The  earliest  demonstrations  of  nonlinear  frequency 
conversion  operated  in  this  phase-mismatched  regime,  and  as  such  were  unable  to 
attain  appreciable  conversion  efficiencies  [8]. 

Until  quasi-phasematching  (QPM)  was  developed,  achieving  phasematching  of 
nonlinear  optical  interactions  was  a  matter  of  finding  materials  with  suitable  birefrin¬ 
gence  for  the  interaction  of  interest.  Many  non-centrosymmetric  crystals  suitable  for 
yd2)  nonlinear  optics  also  have  linear-optical  anisotropy:  the  refractive  index  varies 
as  the  orientations  of  the  polarization  and  propagation  direction  are  changed  with 
respect  to  the  crystalline  axes.  In  some  circumstances,  it  is  possible  to  choose  the 
polarizations  and  k-vectors  of  interacting  waves  such  that  Ak'  =  0,  enabling  a  phase- 
matched  interaction.  This  technique  is  known  as  birefringent  phasematching,  and 
enables  quadratic  growth  of  power  with  propagation  distance,  as  is  shown  by  the 
green  curve  in  Fig.  3.1. 

In  the  case  of  LiNbCU,  the  largest  nonlinear  coefficient  (by  approximately  a  factor 
of  6)  involves  the  interaction  of  three  waves  polarized  along  the  crystalline  z-axis  ,  and 
is  therefore  not  compatible  with  birefringent  phasematching  [12].  The  technique  of 
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Figure  3.1:  Intensity  evolution  of  generated  sum-frequency  radiation  for  three  differ¬ 
ent  phasematching  conditions. 

quasi- phase-matching,  introduced  by  Armstrong  and  Bloembergen  in  1962,  involves 
a  periodic  modulation  of  the  sign  of  the  x^2'  throughout  the  medium  [11],  Allowing 
the  nonlinear  coefficient  to  include  a  spatial  dependence  d(z),  we  find  that  the  SF 
field  at  the  output  of  the  medium  is  given  by  an  integration  of  Eq.  (3.1): 

a2(L)  =  Ta1{Q)ap(Q)  [  d(z)  e~l  Ak' z  dz,  (3.2) 

Jo 

which  (assuming  the  d(z)  is  only  nonzero  on  the  range  [0,  L\)  is  the  Fourier  transform 
of  d(z).  For  periodic  modulation  with  period  A,  d(z)  can  be  constructed  as  a  Fourier 
series: 

OO 

d(z)  =  d0  GmeiK™z ,  (3.3) 

m=— oo 

where  Km  =  27rm/A.  In  the  case  of  periodic  poling,  d(z)  takes  the  values  ±d0 
alternating  in  a  square  wave  fashion  with  duty  cycle  D.  The  Fourier  coefficients  have 
the  values 

Gm  = - sin  mnD 


rmr 


(3.4) 
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Incorporating  the  form  of  Eqs.  (3.3)  and  (3.4)  into  the  integral  in  Eq.  (3.2),  in  the 
case  where  Krn  is  close  to  Ak',  the  integral  is  dominated  by  the  7nth  term  of  the 
summation  in  Eq.  (3.3),  and  thus 


a2(L )  =  rdmai(0)ap(0)e*Afci/2  sine 


(3.5) 


where  dm  =  d0Gm  =  des  is  the  amplitude  of  the  utilized  harmonic  of  the  square  wave 
modulation,  and  A k  =  Ak'  —  Krn  [57] . 

From  Eq.  (3.5),  we  see  that  the  behavior  of  a  QPM  interaction  is  the  same  as  that 
of  a  birefringently  phasematched  interaction  but  where  the  effective  phase  mismatch 
has  been  shifted  by  an  amount  Km,  and  the  nonlinear  coefficient  d  has  been  replaced 
by  dm.  For  first-order  (m  =  1)  QPM,  d\  is  maximized  for  a  50/50  duty  cycle  D  =  0.5, 
in  which  case  d\  =  (2/n)do.  This  case  corresponds  to  a  reorientation  of  d  every  coher¬ 
ence  length  Lc  =  n /  Ak' .  We  can  calculate  the  power  growth  of  the  SF  field  by  a  direct 
integration  of  Eq.  (3.2).  The  power  growth  for  this  case  where  Ak  =  0  is  plotted  as 
the  solid  red  curve  in  Fig.  3.1,  where  we  see,  on  coarse  spatial  scales,  a  quadratic 
growth  of  the  held  as  the  phase  is  reset  every  coherence  length.  The  normalized, 
or  small-signal,  efficiency  of  the  QPM  process  compared  with  a  birefringently  phase- 
matched  interaction  is  reduced  by  a  factor  ^qpmA/bpm  =  (2 / 7r)2  ~  41%.  For  LiNb03, 
the  use  of  quasi-phase-matching  enables  the  use  of  the  d33  =  27  prn/V  nonlinear 
coefficient,  which  allows  for  frequency  conversion  with  efficiencies  approximately  16 
times  larger  than  for  interactions  using  the  largest  birefringently  phasematcheable 
nonlinear  coefficient  d3i  =  4.1  prn/V  [12]. 

Several  techniques  have  been  utilized  to  achieve  quasi-phase-matching  in  nonlinear 
optical  materials,  including  stack-of-plates  methods  [58] ,  orientation-patterned  semi¬ 
conductors  [59],  and  periodically  poled  ferroelectrics  [60].  The  work  described  in  this 
dissertation  takes  advantage  of  electric- field  poling  techniques  [61].  Lithium  niobate, 
the  material  of  choice  for  the  work  of  this  dissertation,  is  a  ferroelectric  material  and 
thus  has  a  spontaneous  electric  dipole  moment.  Applying  an  electric  field  pulse  above 
the  coercive  field  of  the  material  can  reorient  the  crystal  structure  and  change  the 
sign  of  the  d33  nonlinear  coefficient.  The  dynamics  of  the  electric-field  poling  process 
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Figure  3.2:  Typical  voltage  and  current  waveforms  during  whole- wafer  electric-field 
periodic  poling. 

have  been  studied  in  detail  [62],  By  patterning  the  surface  of  the  material,  electric- 
held  poling  enables  a  structuring  of  the  volume  of  the  medium.  Inverted  domains 
propagate  from  the  +z  to  —  z  surface  of  the  wafer,  and  can  be  poled  with  an  aspect 
ratio  (domain  height): (domain  width)  exceeding  100:1. 

The  poling  procedure  we  use  is  outlined  schematically  in  Fig.  3.3.  We  pattern 
the  +z  surface  of  a  0.5-mm-thick,  72. 6-mm- diameter  congruent-composition  LiNbOs 
wafer.  A  thick  layer  (approximately  3  to  5  /irn)  of  photoresist  is  spun  onto  the  wafer 
surface.  Electrodes  are  patterned  via  UV  illumination  through  a  photolithographic 
mask  followed  by  developing  the  pattern  in  LDD-26W  developer  solution.  As  electri¬ 
cal  contact  to  the  wafer  will  be  made  through  a  liquid  electrolyte  (saturated  solution 
of  lithium  chloride  in  water),  the  photoresist  layer  (used  as  an  insulator)  is  hard-baked 
onto  the  wafer.  Metal  electrodes,  while  used  in  initial  poling  experiments  [63],  are 
not  used  here  as  it  has  been  found  that  it  is  difficult  to  remove  all  traces  of  metal 
from  the  wafer  surface,  which  is  necessary  here  as  waveguides  are  formed  very  close 
to  the  +z  surface.  The  wafer  is  placed  into  an  electrical  circuit  where  the  +z  surface 
receives  a  high-voltage  pulse  created  using  a  function  generator  followed  by  a  Trek 
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x2000  high-voltage  amplifier,  and  the  —  z  surface  is  connected  to  ground  through  a 
1  kh2  resistor.  Poling  requires  electric  fields  exceeding  the  the  LiNbCU  coercive  field  of 
21  kV /mm,  which  for  the  0.5- nun-thick  samples  indicates  a  voltage  of  10.5  kV.  The 
voltage  across  the  1  kfl  resistor  is  monitored,  giving  a  measurement  of  the  poling  cur¬ 
rent.  The  applied  voltage  pulse  (measured  by  the  output  monitor  of  the  amplifier)  is 
shown  in  black  in  Fig.  3.2,  and  consists  of  a  ramp  up  to  a  field  level  slightly  below  the 
coercive  field,  a  voltage  spike  above  the  coercive  field  to  aid  in  nucleation  of  inverted 
domains,  and  a  voltage  hold  at  a  level  slightly  above  the  coercive  field.  As  is  seen  in 
Fig.  3.2,  once  the  domains  start  to  nucleate,  current  begins  to  flow  through  the  wafer 
as  the  ferroelectric  domains  reorient.  Poling  is  done  in  the  self-terminating  regime, 
whereby  the  growth  of  domain  walls  is  stopped  at  a  certain  duty  cycle  D  determined 
by  the  poling  initial  conditions  [64].  We  have  found  that  poling  quality  (uniformity  of 
domain  pattern)  is  highest  under  self-terminating  conditions.  Unlike  previous  work 
[55],  in  which  the  poling  current  I(t)  was  limited  by  the  amplifier’s  current  limit  of 
20  mA,  in  Fig.  3.2  we  see  current  levels  in  the  range  of  a  few  mA  so  the  amplifier 
is  still  in  the  linear  regime.  When  one  designs  a  photolithographic  mask  for  poling, 
one  can  calculate  the  expected  area  Apol  to  be  poled,  which  in  turn  be  related  to  the 
total  charge  Q  transported  through  the  wafer  as 

4*>i  =  TT  =  JT  f  I(t)  dt  (3-6) 

where  Psp  =  78  fiC  cm-2  is  the  spontaneous  polarization  charge  density  associated 
with  the  spontaneous  dipole  moment  in  LiNbCU  [55].  For  various  projects  throughout 
this  dissertation,  we  have  poled  a  range  of  periods  between  7  and  25  /mi  using  these 
techniques,  sometimes  even  within  the  same  waveguide  [65]. 

3.2  Fabrication  of  RPE  waveguides 

Nonlinear  optical  interactions  can  be  performed  most  efficiently  when  the  interacting 
fields  are  confined  to  a  small  area,  hence  high  intensity  and  high  field  strength,  over 
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Figure  3.3:  Process  flow  for  fabrication  of  RPE  PPLN  waveguides. 


long  distances.  Performing  frequency  conversion  in  waveguides  has  two  distinct  ad¬ 
vantages.  (1)  In  bulk  nonlinear  optical  media,  diffraction  sets  limits  on  the  minimum 
spot  size  for  a  given  interaction  length;  an  analysis  of  coupled-wave  interactions  with 
focused  Gaussian  beams  indicates  an  optimum  efficiency  for  near-confocal  focusing 
[66].  Therefore,  for  bulk  interactions,  the  area  of  an  optimally  focused  beam  is  propor¬ 
tional  to  the  length  of  the  crystal,  so  the  conversion  efficiency  7]  increases  only  linearly 
with  sample  length.  (2)  In  waveguides,  the  interactions  involve  discrete  eigenmodes 
of  the  structure.  For  bulk  interactions,  conversion  efficiency  is  a  function  of  intensity 
and  therefore  varies  across  a  beam  profile,  making  it  difficult  to  attain  near  unity 
conversion  efficiencies,  while  for  waveguides,  an  entire  mode  converts  as  a  single  en¬ 
tity.  This  aspect  of  nonlinear  interactions  in  waveguides  has  enabled  extremely  high 
signal  depletions  in  SFG,  at  the  level  of  41  dB  [67],  a  level  only  achievable  with  very 
large  pump  beams  in  bulk  devices  which  therefore  have  a  very  large  value  of  Pmax  for 
complete  conversion. 

There  are  several  techniques  for  forming  waveguides  in  PPLN,  including  Ti  indif¬ 
fusion  [68],  mechanical  ridge  waveguides  [69],  and  annealed  proton  exchange  (APE) 
waveguides  [70].  In  this  work  we  have  focused  on  annealed  and  reverse  proton  ex¬ 
change  (RPE)  waveguides  [71],  which  have  substantially  lower  propagation  losses  and 
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higher  mode-overlap  integrals  than  APE  waveguides,  making  them  the  platform  of 
choice  for  applications  where  high  throughput  and  conversion  efficiency  are  essential. 

Following  preparation  of  PPLN  wafers  via  the  procedure  outlined  in  Sec.  3.1,  we 
form  annealed  and  reverse  proton  exchange  waveguides  via  the  steps  diagrammed 
schematically  in  Fig.  3.3.  RF  sputtering  of  an  80  nm  SiC>2  layer  followed  by  a  second 
photolithography  step  is  used  to  create  channels  through  which  protons  enter  the  sub¬ 
strate.  Following  the  lithography  and  development  of  the  masking  pattern,  channels 
are  etched  in  the  SiC>2  mask  using  a  20:1  buffered-oxide  wet  etch.  We  use  molten 
benzoic  acid,  buffered  with  approximately  5  wt  %  lithium  benzoate,  at  a  temperature 
of  between  170  and  180°  C,  as  a  proton  source.  Following  proton  exchange,  samples 
are  annealed  in  air  at  310°  C  for  a  typical  duration  of  22.5  hrs.  Following  removal 
of  the  SiC>2  mask,  reverse  exchange,  by  which  protons  are  removed  from  the  material 
and  exchanged  with  lithium  ions,  is  done  in  a  eutectic  melt  of  LiNOs/NaNC^/KNOs, 
at  temperatures  of  approximately  300°  C.  The  proton  dose  may  be  measured  in  situ 
via  absorption  spectrophotometry  of  an  O-H  absorption  resonance  at  a  wavelength 
of  2.85  /irn.  A  detailed  model  of  the  concentration-dependent  nonlinear  diffusion  of 
protons  in  LiNbC>3  has  been  carried  out  in  the  dissertation  work  of  Roussev  [56]. 
Combining  the  diffusion  model  with  a  model  for  the  dispersion  of  HxLi1_xNb03  and 
a  2D  waveguide  eigenmode  solver  has  enabled  a  platform  via  which  critical  design 
parameters  may  be  simulated. 

The  existence  of  a  fabrication  model  has  enabled  accurate  prediction  of  linear 
and  nonlinear  optical  properties  of  RPE  PPLN  waveguides,  enabling  the  inclusion  of 
complex  integrated  optical  structures  such  as  directional  couplers,  tight-radius  bends, 
and  mode  multiplexers  and  demultiplexers.  For  most  of  the  experiments  described  in 
this  dissertation,  the  waveguides  used  have  schematic  designs  as  shown  in  Fig.  3.4. 

At  the  input  of  the  device  is  a  mode  filter,  a  single-mode  waveguide  in  which  the 
mode  shape  is  generally  designed  to  closely  match  that  of  an  optical  fiber  used  for 
input  coupling.  As  RPE  waveguides  and  optical  fibers  have  approximately  the  same 
core-cladding  index  contrast,  it  is  possible  to  achieve  spatial  mode  matching  of  ap¬ 
proximately  90%.  For  RPE  PPLN  waveguides,  there  exists  a  noncritical  waveguide 
width  for  which  the  QPM  wavelength  is  first-order  insensitive  to  waveguide  width 
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Figure  3.4:  Schematic  of  an  RPE  PPLN  waveguide. 


variations  [72] .  The  width  of  the  mixing  waveguide  is  chosen  to  be  noncritical  for  the 
interaction  of  interest,  and  typically  occurs  for  channel  widths  between  7  and  9  /irn. 
The  mode  Liter  therefore  serves  a  second  purpose;  the  adiabatic  taper  converts  the 
fundamental  (TM0o)  spatial  mode  of  the  mode  Liter  to  the  TM0o  mode  of  the  mixing 
waveguide,  such  that  all  of  the  power  in  the  input  signal  can  participate  in  the  con¬ 
version  process.  At  the  output  of  the  device,  an  additional  adiabatic  taper  and  mode 
Liter  can  be  included  to  enable  high-efbciency  outcoupling  into  another  optical  Lber. 


3.3  Device  characterization  and  packaging 

In  this  section  we  describe  the  techniques  by  which  RPE  PPLN  waveguides  are  char¬ 
acterized  and  prepared  for  use  in  systems  experiments.  A  typical  waveguide  chip  is 
5.9  mm  wide,  and  contains  approximately  100  candidate  waveguides  spaced  by  50  /irn. 
A  processing  run  consisting  of  devices  from  two  LiNbC>3  wafers  can  therefore  entail 
over  1500  individual  waveguides  to  test.  Therefore,  a  screening  process  is  important 
to  select  the  best-performing  waveguide  for  a  particular  application. 

3.3.1  Loss  Measurements 

The  Lrst  method  by  which  waveguides  are  screened  is  a  measurement  of  waveguide 
propagation  loss.  For  frequency  conversion  of  single  photons,  high  device  throughput 
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is  of  paramount  importance.  We  typically  measure  propagation  loss  via  the  Fabry- 
Perot  technique  [73].  A  single-spatial-mode  waveguide  with  uncoated  facets  with 
transmittivity  T  polished  orthogonal  to  the  waveguide  forms  a  monomode  low-finesse 
Fabry- Perot  resonator.  The  transmitted  power  Pt  through  the  waveguide  is  given  by 

rp2  -aL 

Pt  =  - - —2 - y - PoVc  (3.7) 

(l  —  R)  +  RR sin2  (0/2) 


where  Pq  is  the  input  power,  ?/c  is  the  input  coupling  efficiency,  and  0  =  2 (3L  is  the 
round-trip  phase  accumulation,  where  0  the  propagation  constant  and  L  is  the  total 
device  length,  and  R  =  Rexp(-aL).  If  one  has  knowledge  of  the  facet  reflectivity  R, 
by  measuring  the  contrast  ratio 
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a  function  only  of  R,  one  can  determine  the  propagation  loss  a. 
(3.7)  and  (3.8)  yields 


Combining  Eqns. 


(3.9) 


For  weakly  confining  waveguides,  the  mode  reflectivity  is  found  using  the  effective 
index  of  the  mode  neg  as  R  =  (neg  —  l)2/ (^efr  +  l)2-  For  RPE  waveguides  in  LiNb03, 
R  «  0.1324  for  A  =  1.55  /irn.  The  propagation  loss  a,  measured  in  dB  cm-1  is  then 


found  as 


4.34 


a  = 


In  R  —  In  R 


(3.10) 


We  generally  have  measured  propagation  losses  in  the  1.55-/im  telecom  band.  By 
sweeping  the  wavelength  of  a  tunable  single-frequency  laser,  the  propagation  phase  0 
can  be  varied  continuously.  For  a  typical  device  of  length  L  =  5.2  cm,  the  free  spectral 
range  of  the  Fabry-Perot  cavity  is  A /  =  1.3  GHz,  corresponding  to  a  wavelength 
range  AA  =  11  pm.  A  sample  transmission  curve  versus  wavelength  is  shown  in 
Fig.  3.5,  where  we  observe  interference  fringes  with  average  contrast  ratio  K  =  0.243, 
corresponding  to  a  propagation  loss  of  0.09  dB/cm.  The  best  waveguides  typically 
have  propagation  losses  in  the  range  0.05  to  0.1  dB/cm,  although  on  an  individual 
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Figure  3.5:  Transmission  vs.  wavelength  for  a  5.2-cm  RPE  PPLN  waveguide.  The 
observed  fringe  contrast  of  0.243  corresponds  to  a  propagation  loss  a  =  0.09  dB  cm-1. 


chip  it  is  not  uncommon  to  find  a  large  variation  in  waveguide  quality,  which  can 
often  be  traced  back  to  lithographic  defects,  or  to  scratches  on  the  facet  formed 
during  polishing.  We  note  that  the  expected  error  in  loss  |  Aa|  can  be  calculated  from 
Eq.  (3.7)  by  the  techniques  of  error  propagation  [74]  as 


|Aa| 


4.34  |  A  AT  | 
K~ 


(3.11) 


For  our  cases,  we  have  fractional  uncertainty  in  the  contrast  \AK\/K  «  2%,  which 
implies  a  value  |Aa|  ~  0.02  dB/cm.  For  accurate  loss  measurements,  one  must  be 
careful  that  the  any  detector  dark  voltage  or  signal  due  to  stray  light  is  subtracted 
off  before  calculating  the  contrast  ratio. 
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Figure  3.6:  Phasematching  tuning  curves  for  SFG,  showing  generated  target  power 
at  A2  =  832  nm  versus  wavelength  Ai,  showing  (a)  near-ideal  and  (b)  highly  distorted 
tuning  behavior. 

3.3.2  Phasematching  characterization 

The  measurement  of  a  QPM  tuning  curve  for  a  waveguide  is  a  useful  diagnostic 
of  waveguide  quality.  In  the  low-conversion-efficiency  limit,  the  tuning  curve  shape 
is  given  by  P2(Afc)  =  sine2  (AkL/2).  However,  any  non-ideal  waveguide  will  have 
nonuniformities  in  waveguide  width  and  poling  quality  that  will  cause  the  shape  of 
the  tuning  curve  to  vary  from  this  ideal  shape.  Therefore,  finding  waveguides  whose 
phasematching  behavior  is  near-ideal  is  an  effective  method  of  screening  for  device 
quality. 

As  most  of  the  projects  in  this  dissertation  have  focused  on  interactions  involving 
radiation  in  or  near  the  1.55-/mi  telecommunications  band,  we  can  make  use  of  the 
well-developed  external-cavity  diode  lasers  (ECDLs),  which  can  be  tuned  over  large 
bandwidths  with  no  mode  hops,  enabling  a  single-shot  measurement  of  the  phase- 
matching  tuning  curve.  The  use  of  computer  control  via  GPIB  and  data  digitization 
using  a  data  acquisition  card  enabled  rapid  screening  of  candidate  waveguides. 

Two  example  phasematching  tuning  curves  are  shown  in  Fig.  3.6.  In  this  case, 
single-frequency  radiation  at  Xp  =  1.8  /irn  was  combined  with  the  telecom  band 
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Figure  3.7:  Closeup  image  of  a  fiber  pigtail. 


signal  at  Ai  in  a  fiber-optic  wavelength  combiner.  The  combined  beams  were  then 
launched  using  a  fiber  collimator  and  coupled  into  the  waveguide  using  an  aspheric 
lens  (/  =  8  mm).  The  waveguides  under  test  were  designed  for  sum-frequency  gen¬ 
eration  to  produce  target  radation  at  A2  =  (A71  +  A~  :)  1  =  832  nm.  The  generated 
target  radiation  was  collimated  using  a  second  aspheric  lens  and  detected  using  an 
amplified  Si  photodiode.  Fig.  3.6(a)  shows  a  near-ideal  tuning  curve,  while  the  curve 
in  Fig.  3.6(b)  shows  significant  distortions. 

3.3.3  Pigtailing 

In  the  previous  sections  we  have  described  the  techniques  used  for  fabricating  and 
characterizing  RPE  PPLN  waveguides  for  efficient  nonlinear  optics.  Owing  to  the 
small  size  of  the  optical  mode  of  these  waveguides,  their  use  in  systems  experiments 
requires  some  technique  to  maintain  stable  alignment.  We  have  designed  the  input 
mode  Liters  to  support  optical  modes  which  closely  match  those  of  optical  fiber, 
and  therefore  can  bypass  these  alignment  difficulties  by  attaching  or  “pigtailing”  an 
optical  fiber  ferrule  to  either  the  input  or  output  of  the  waveguide,  or  both. 
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Fig.  3.7  shows  a  closeup  image  of  a  fiber  pigtail.  A  SMF-28  optical  fiber  is  packaged 
inside  of  a  1-mm-diameter  glass  ferrule  and  flat  polished.  The  ferrule  is  carefully 
aligned  to  the  desired  waveguide  while  monitoring  the  transmission  of  a  probe  optical 
signal.  To  ensure  mechanical  stability  during  the  alignment  procedure,  the  ferrule 
is  held  in  a  mechanical  clamp  and  the  waveguide  is  held  to  a  working  surface  using 
a  vacuum  clamp.  Alignment  is  done  using  piezoelectric  actuators.  As  an  adhesive 
we  use  UV-curing  optical-grade  epoxy  (for  example,  Norland’s  NOA-61).  A  tiny 
droplet  of  epoxy  is  applied  between  the  ferrule  and  the  waveguide  using  an  optical 
fiber  tip.  The  ferrule  is  then  brought  into  contact  with  the  waveguide,  and  then  cured 
using  a  handheld  UV  gun.  For  C-band  signals,  we  can  typically  achieve  coupling  loss 
(including  the  4%  Fresnel  reflection  between  fiber  and  LiNbOs)  of  approximately  0.5 
dB.  We  have  also  pigtailed  waveguides  designed  for  other  wavelengths,  but  because 
loss  measurements  are  typically  done  at  1550  nm,  it  can  be  difficult  to  decouple 
propagation  and  coupling  loss. 


Chapter  4 

Noise  processes  in  quantum 
frequency  converters 


The  theory  of  quantum  frequency  conversion  was  described  in  Chapter  2,  and  a  device 
platform  for  achieving  high-efficiency  frequency  conversion  was  described  in  Chapter 
3.  An  important  aspect  of  the  theory  described  in  Chapter  2  is  that  the  conversion 
process  is,  in  principle,  noiseless.  We  recall  here  from  Eq.  (2.19)  that  the  conversion 
efficiency  in,  for  example,  an  upconverter  for  an  input  signal  uj\  pumped  at  ujp ,  is 
given  by 


V 


(A f2(-  -  L)) 

(N,(z  =  0)) 


sm  — 


(4.1) 


where  Nj  =  Pj/(huj )  and  Pp  is  the  pump  power  [30],  and  we  have  assumed  negligible 
propagation  loss  and  pump  depletion.  For  RPE  PPLN  waveguides,  a  typical  value 
of  r/nor  is  approximately  100%/W/cm2,  which,  for  a  device  of  length  L  —  5  cm  gives 
a  value  Pmax  =  100  mW.  This  power  requirement  means  that,  for  a  pump  at,  say, 
1.06  pm,  approximately  5  x  101'  pump  photons  per  second  are  necessary  to  reach 
complete  conversion. 

Since,  even  for  an  extremely  efficient  converter,  such  a  huge  pump  photon  flux 
is  required  for  substantial  conversion,  the  converter  is  highly  susceptible  even  to 
very  weak  scattering  processes.  Any  photon  within  the  acceptance  bandwidth  of  the 
frequency  converter,  whether  present  as  an  input  signal  or  generated  via  an  inelastic 
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Figure  4.1:  Schematic  of  noise  processes  in  (a)  short-wavelength-pumped  and  (b) 
long-wavelcngth-pumped  quantum  frequency  upconverters. 


spontaneous  scattering  event  due  to  the  strong  pump,  will  register  as  a  signal.  For 
instance,  for  a  scattering  process  with  a  scattering  probability  of  ICC12,  one  would 
have  a  noise  floor  (for  signal-to-noise  ratio  of  1)  of  5  x  105  signal  photons/s,  a  large 
signal  on  “single-photon”  scales. 

In  this  chapter,  we  seek  to  describe  the  two  main  spontaneous  scattering  processes 
that  have  been  limiting  factors  in  the  performance  of  QFC  devices.  A  schematic  of 
the  noise  processes  we  will  address  is  shown  in  Fig.  4.1.  A  first  such  process  is  spon¬ 
taneous  parametric  downconversion  (SPDC).  In  Sec.  4.1,  we  show  that  nonidealities 
in  the  QPM  grating  lead  to  a  “white”  phasematching  noise  floor,  such  that  SPDC 
occurs  in  a  broadband  fashion  for  frequencies  lower  than  u>p.  In  Sec.  4.2  we  address 
the  second  scattering  process:  spontaneous  Raman  scattering  (SRS).  SRS  involves 
either  the  creation  (for  the  Stokes  process)  or  annihilation  (for  the  anti-Stokes)  of  an 
optical  phonon  in  the  nonlinear  medium.  SRS  is  also  found  to  occur  over  a  broad 
bandwidth.  In  Sec.  4.3  we  will  give  a  summary  of  the  noise  processes  and  outline 
some  constraints  for  designing  QFC  devices  with  low  noise.  The  theoretical  descrip¬ 
tions  throughout  this  chapter  will  be  semiclassical,  where  our  goal  is  to  understand 
the  factors  influencing  noise  performance.  Devices  for  low-noise  upconversion  will 
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be  discussed  in  Chapter  5,  and  an  application  for  low-noise  downconversion  will  be 
discussed  in  Chapter  6. 

We  also  note  that  in  addition  to  the  scattering  processes  discussed  here,  which 
are  fundamental  in  origin,  there  are  additional  technical  noise  sources  that  must  be 
overcome  in  any  system  implementation.  These  include  noise  photons,  for  example, 
emitted  directly  from  the  pump  source  used  in  an  experiment.  We  have  studied  the 
noise  emission  from  both  erbium-  and  thulium-doped  laser  systems,  which  will  be 
addressed  in  Chapter  5.  One  must  also  take  precautions  by  shielding  the  devices 
from  environmental  stray  light  that  can  also  manifest  as  noise  photons. 

Our  work  on  characterization  of  SPDC  noise  and  its  enhancement  due  to  QPM 
grating  disorder  has  been  published  in  [48],  and  measurements  of  second- harmonic 
generation  and  domain  characterization  studies  have  been  published  in  [75].  A  study 
of  Raman  scattering  was  included  in  our  publication  [67]. 

4.1  QPM  random  duty-cycle  errors:  definitions, 
characterization,  and  nonlinear  optical  effects 

QPM  gratings  are  subject  to  a  wide  variety  of  fabrication  errors,  which  can  influence 
the  nonlinear  optical  performance  of  a  device.  One  such  error  type,  random  duty- 
cycle  (RDC)  errors,  are  illustrated  schematically  in  Fig.  4.2.  We  consider  a  grating 
with  period  Aq,  such  that  the  domains  are  inverted  every  l  =  Ag/2,  where,  for 
first-order  QPM,  l  —  Lc  —  n / Ah' .  In  a  device  with  random  duty-cycle  errors,  the 
position  of  the  n’th  domain  boundary  zn  deviates  from  its  ideal  position  by  an  amount 
5zn  =  zn  —  znfi,  where  the  ideal  domain  boundary  position  znio  =  nl.  The  deviation 
of  the  length  of  the  nth  domain  is  thus  given  as  5ln  —  ln  —  l,  where  ln  =  zn+\  —  zn.  In  a 
device  with  RDC  errors  (and  no  random  period  errors),  the  expectation  value  (Sln)  is 
independent  of  the  domain  number  n,  corresponding  to  a  preservation  of  long-range 
order  inherent  in  the  fabrication  of  the  grating  by  lithographic  methods.  We  assume 
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Figure  4.2:  Schematic  of  QPM  domains  showing  random  duty-cycle  errors. 


a  Gaussian  probability  distribution  for  the  ln: 


y/Trrai 


exp 


1  fL-l 


2  \  (Ti 


(4.2) 


where  we  have  defined  the  domain-length  variance  as  of. 


4.1.1  Theory  of  SHG  in  disordered  QPM  gratings 

Following  [57],  we  consider  SHG  in  the  limit  of  low  conversion  efficiency.  In  the 
following,  the  subscript  1  indicates  the  fundamental  frequency  (FF)  and  2  indicates 
the  second-harmonic  (SH).  In  an  ideal  first-order  QPM  grating,  the  intrinsic  wave- 
vector  mismatch  A k'  =  k2  —  2 k\  is  compensated  by  a  periodic  reversal  in  the  sign 
of  the  nonlinear  coefficient  d(z)  every  coherence  length  lc  =  n/Ak'.  For  a  device  of 
length  L,  the  SH  held  at  z  =  L  relative  to  a  device  with  an  ideal  first-order  QPM 
grating  is 

a2(L)  =  [  g(z)eiAk'z  dz,  (4.3) 

Jo 

where  we  have  introduced  the  normalization  g(z)  =  (7r/2L)d(z)/do  such  that  the 
Fourier  transform  g(Ak')  represents  the  SHG  relative  to  an  ideal  QPM  grating,  as  we 
have  previously  shown  in  Sec.  3.1  [75]. 

For  a  QPM  grating  with  random  duty-cycle  (RDC)  errors,  we  can  write  g(z)  as 
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follows: 


9(z) 


7 r 

2L 


N 

Ef-1)”11 

n=  1 


where  II ((x  —  x0)/w)  is  a  shifted  and  scaled  rectangle  function 


(4.4) 


x  <  x0  —  w/2, 

xq  —  w/2  <  x  <  x0  +  w/2,  (4.5) 

x  >  xq  +  w/2. 

We  are  interested  in  the  quantity  (\g(Ak')\2) ,  the  efficiency,  relative  to  an  ideal  grat¬ 
ing,  of  SHG  at  a  phase  mismatch  Ah',  where  the  expectation  value  is  over  an  ensemble 
of  gratings  with  randomly  perturbed,  independent  and  identically  distributed  (HD) 
zn.  For  deviations  from  phasematching  small  compared  to  an  inverse  domain  length, 
i.e.  |Afc|  =  |  A  A:'  —  n/l\  <C  l~x,  we  find 

(|9(Afc)|2)  =  /(< T^sinc2  (/ADj  +  (4.6) 

where  the  function  /(cq)  =  exp(— TT2a/ /2l2).  A  complete  derivation  of  Eq.  4.6,  valid 
not  only  for  SHG  but  for  difference-frequency  generation,  is  given  in  Appendix  A. 

The  relative  SHG  efficiency  (|^(Afc)|2)  calculated  in  Eq.  (4.6)  consists  of  two 
components.  The  first  represents  the  sinc2(AA:L/2)  tuning  behavior  of  a  QPM  device 
with  an  efficiency  reduced  by  the  amount  /(crj),  an  effect  discussed  in  [57].  The 
second  component  is  a  flat  QPM  pedestal  of  height  (1  —  f(ai))/N  that  is  independent 
of  A k.  We  note  that  as  the  grating  disorder  a  is  reduced,  /(cq)  approaches  1,  and 
the  device  regains  the  characteristics  of  an  ideal  grating.  Conversely,  as  the  grating 
disorder  is  increased,  one  encounters  the  behavior  of  the  random  QPM  devices  studied 
in  [76,  77],  in  which  the  growth  of  the  SH  field  scales  as  L1/2,  i.e.  the  field  grows  in 
a  random-walk  fashion.  The  pedestal  height  can  be  estimated  to  within  numerical 
factors  by  a  simple  physical  argument:  on  average,  the  portion  of  a  domain  in  which 
the  SH  evolution  occurs  randomly  is  cq  =  eq/7.  Therefore,  the  intensity  of  the  random 
component  of  the  SHG  scales  as  (d/4/iV)2,  which  when  divided  by  the  N2  dependence 
of  ideal  SHG  intensity,  gives  a’f/N. 


CHAPTER  4.  NOISE  PROCESSES  IN  QFC  DEVICES 


37 


Ak  L/tu 

Figure  4.3:  Simulation  and  analytic  theory  of  SHG  in  a  QPM  grating  with  domain 
disorder.  Parameters:  di  =  0.2,  L  —  1  cm,  ensemble  average  over  100  disordered 
gratings. 

To  test  the  analytical  theory,  we  performed  both  numerical  simulations  and  ex¬ 
periments.  Our  experiments  on  SHG  in  an  RPE  PPLN  waveguide  will  be  discussed 
in  Sec.  4.1.3.  Numerically,  we  constructed  an  ensemble  of  100  disordered  QPM  grat¬ 
ings,  with  domain  lengths  drawn  from  a  Gaussian  distribution  and  with  normalized 
pseudo-random  duty-cycle  error  &i  =  20%.  The  simulation  was  for  a  device  of  length 
L  =  1  cm  with  a  poling  period  of  Aq  —  16  /im,  suitable  for  1.55-/im  SHG,  and  the 
total  number  of  domains  was  N  =  L/^Aq/2)  =  1250.  The  SHG  tuning  behavior  is 
calculated  by  a  direct  integration  of  Eq.  4.3  for  each  grating.  The  SHG  tuning  be¬ 
havior  for  each  grating  is  normalized  to  the  performance  of  an  ideal  grating,  and  an 
ensemble  average  is  computed.  This  ensemble  average  is  plotted  in  Fig.  4.3  (dashed 
red  curve)  along  with  the  prediction  of  the  analytical  theory  (solid  blue  curve),  show¬ 
ing  excellent  agreement  on  the  height  of  the  QPM  pedestal.  For  frequencies  far  from 
the  QPM  peak,  the  pedestal  results  in  elevated  generation  efficiency  when  compared 
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with  the  ideal  QPM  grating  (solid  black).  Also  shown  are  the  results  for  a  single 
disordered  QPM  grating  (thin  solid  green).  We  note  that  while  the  ensemble  average 
of  the  pedestal  is  white,  an  individual  grating  can  exhibit  large  fluctuations  in  effi¬ 
ciency  far  from  the  QPM  peak.  Properties  of  the  statistics  of  the  pedestal  height  are 
discussed  in  Appendix  A. 

4.1.2  Characterization  of  RDC  errors 

A  direct  study  of  random  duty-cycle  errors  via  microscopy  or  other  means  would 
enable  an  estimate  of  the  expected  QPM  pedestal  height  (|//|2)  in  real  devices.  We 
analyzed  the  severity  of  random  duty-cycle  errors  in  an  electric-field  periodically  poled 
LiNbCh  sample.  We  poled  a  76-mm- diameter  0.5-mm-thick  lithium  niobate  wafer  via 
the  techniques  described  in  Sec.  3.1.  To  decorate  the  inverted  domains,  we  etched 
the  original  +£  surface  of  the  wafer  in  hydrofluoric  acid  (HF)  for  20  min.  HF  etches 
the  +z  domain  orientation  at  a  higher  rate  than  the  —  z  orientation,  and  thus  the 
etch  provides  topographic  contrast  between  the  two  domain  orientations,  which  can 
be  observed  optically.  We  chose  to  measure  the  domain  statistics  on  the  original  +z 
surface  of  the  wafer  as  this  is  the  side  upon  which  the  waveguides  are  formed.  As 
the  domains  do  not  propagate  perfectly  straight  through  the  substrate  during  poling, 
the  domain  positions  on  the  +z  surface  will  provide  the  best  representation  of  the 
domain  disorder  “experienced”  by  a  signal  propagating  through  a  RPE  waveguide. 

We  next  sought  a  technique  whereby  we  could  accurately  measure  the  positions 
of  a  large  number  of  ferroelectric  domain  boundaries,  such  that  statistics  of  their  po¬ 
sitions  and  the  domain  sizes  could  be  obtained.  An  effective  measurement  technique 
would  provide  lateral  resolution  below  0.1  /im,  while  providing  enough  field  of  view 
(>  100  /irn)  to  measure  a  large  number  of  domains  simultaneously.  The  first  mea¬ 
surement  technique  we  explored  was  atomic  force  microscopy  (AFM)  [78],  which  has 
excellent  spatial  resolution  due  to  the  extremely  sharp  AFM  tip  (a  few  nm,  easily) 
but  limited  field  of  view  due  to  the  fact  that  spatial  scanning  is  done  using  piezo¬ 
electric  actuators.  A  second  technique  widely  used  to  measure  topographic  contrast 
is  surface  prohlometry,  in  which  a  metal  stylus  is  dragged  over  a  surface  and  the 
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surface  topography  is  recorded.  Profilometry  can  give  very  long  linear  travel  ranges 
(hundreds  of  /im  or  even  a  few  mm)  with  very  high  (sub-nm)  depth  resolution.  The 
disadvantage  of  such  techniques  is  that  the  stylus  tip  diameter  is  usually  on  the  order 
of  10  n m,  on  the  order  of  a  typical  domain  size  for  a  PPLN  device.  As  we  were  seeking 
to  characterize  domains  with  l  «  5  /im,  profilometry  was  not  an  appropriate  choice, 
but  could  be  used  for  very-long-pitch  QPM  gratings. 

The  most  effective  domain  characterization  techniques  were  based  on  optical  mi¬ 
croscopy.  One  can  view  the  domain  pattern  directly  using  a  high-magnification  ob¬ 
jective  lens.  Provided  any  imaging  abberations  can  be  corrected  (for  example,  by 
imaging  a  lithographic  mask  as  a  calibration  sample),  using  image  processing  meth¬ 
ods  one  can  accurately  determine  the  domain  boundary  positions.  Microscopy  has 
been  used  since  earliest  days  of  periodic  poling  to  measure  statistics  of  domains  sizes 
[79].  While  microscopy  is  effective,  our  best  results  were  obtained  using  Fizeau  (Zygo) 
interferometry,  which  can  more  easily  provide  quantitative  information  on  the  depth 
as  it  is  measured  by  interfering  the  wavefront  off  the  sample  to  be  characterized  with 
the  perfectly  flat  wavefront  from  a  reference  surface  [80]. 

We  obtained  Zygo  micrographs  at  7  positions  across  a  76-mm-diameter  PPLN 
wafer.  A  sample  image  is  shown  in  Fig.  4.4.  Using  a  50  x  objective,  we  have  a  field  of 
view  of  150  x  113  /mi  in  a  640  x  480  pixel  image.  To  determine  the  domain  size,  we 
first  projected  each  domain  row  onto  an  axis  parallel  to  Kg,  which  formed  a  line  trace 
of  surface  height  for  each  domain.  We  numerically  defined  the  domain  boundary  as 
the  position  at  which  the  surface  height  is  midway  between  the  minimum  (in  the  +z 
areas)  and  maximum  (on  the  —  z  poled  islands).  After  analyzing  all  the  images,  we 
formed  a  histogram  of  the  inverted  domain  sizes:  the  results  are  shown  in  Fig.  4.5. 

The  wafer  we  analyzed  was  poled  with  a  grating  period  A^  =  10.22  /mi.  We 
observed  a  mean  inverted-domain  width  of  5.90  /im,  which  for  the  poling  which  cor¬ 
responds  to  a  duty  cycle  of  58%.  From  the  histogram  we  also  calculate  cq  =  0.41  /mi, 
which  corresponds  to  normalized  RDC  error  of  cq  =  8.0%.  The  solid  line  in  Fig.  4.5 
is  a  fit  to  a  Gaussian  distribution,  which  shows  good  agreement  with  the  observed 
histogram,  supporting  the  assumptions  made  in  the  development  of  the  analytical  the¬ 
ory.  The  only  significant  deviation  from  Gaussian  statistics  is  a  longer  tail  towards 
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Figure  4.4:  Surface  topography  of  etched  PPLN  wafer  measured  using  Zygo  interfer¬ 
ometry,  showing  inverted  ferroelectric  domains  with  random  duty-cycle  errors.  Color 
bar  scale  in  /mi. 


larger  domains.  To  test  whether  the  tail  has  a  significant  influence  on  the  behavior 
of  the  QPM  pedestal  we  performed  SHG  simulations  with  domain  length  statistics 
drawn  from  the  distribution  we  observed  in  Fig.  4.5  and  compared  the  results  to  those 
done  using  a  Gaussian  distribution  with  the  same  average  duty  cycle.  The  computed 
efficiency  far  from  the  QPM  peak  agreed  to  within  approximately  2%  using  both  sets 
of  domain  statistics. 


4.1.3  Observation  of  QPM  pedestal 

We  next  sought  to  characterize  the  effects  of  RDC  errors  on  SFG  far  from  phasematch¬ 
ing  in  a  real  PPLN  sample.  Because  the  theory  and  simulations  described  above  both 
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Figure  4.5:  Histogram  of  inverted  domain  sizes  l  for  a  PPLN  wafer  with  Ac  = 
10.22  /mi,  with  Gaussian  fit.  (1)  =  5.90  /mi,  with  cr t  =  0.41  /mi. 

represent  plane-wave  physics,  we  performed  experiments  using  guided-wave  interac¬ 
tions,  which  are  mathematically  isomorphic  to  plane-wave  interactions.  We  fabricated 
reverse-proton-exchanged  PPLN  waveguide  devices  with  59-mm-long  QPM  gratings. 
The  waveguides  incorporate  mode  Liters  and  adiabatic  tapers  to  facilitate  launch¬ 
ing  the  FF  fundamental  spatial  mode  [30].  The  output  facet  of  the  waveguide  was 
polished  at  an  angle  of  5°  to  prevent  Fabry-Perot  interference  effects.  The  poling 
period  was  15.9  /mi,  which  resulted  in  a  first-order  QPM  peak  at  1533.75  nm  at  a 
temperature  of  25°  C,  as  seen  in  Fig.  4.7  Inset  1.  The  experimental  setup  is  shown 
schematically  in  Fig.  4.6.  In  the  experiment,  we  used  two  external  cavity  diode  lasers 
to  cover  the  90-nm  measurement  range.  The  FF  light  was  modulated  with  a  me¬ 
chanical  chopper  and  was  was  coupled  in  and  out  of  the  waveguide  using  aspheric 
lenses  (/  =  8  mm),  and  the  outgoing  SH  and  residual  FF  light  were  separated  us¬ 
ing  a  dichroic  mirror.  The  SH  light  at  the  chopper  frequency  was  detected  using  an 
unbiased  Si  photodiode  connected  to  the  current  input  of  a  lock-in  amplifier,  and 
the  outgoing  FF  light  was  monitored  using  a  calibrated  InGaAs  detector.  The  FF 
wavelength  was  stepped  in  increments  of  0.01  nm,  and  the  sensitivity  of  the  lock-in 
amplifier  was  adjusted  depending  on  the  SH  signal  level  to  enable  a  high  dynamic 
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Figure  4.6:  Experimental  setup  for  observation  of  SHG  QPM  pedestal. 


range  measurement. 

The  normalized  SHG  efficiency  (|g(A)|2)  was  computed  from  the  measured  SH 
and  FF  data,  and  is  plotted  as  the  dashed  line  in  Fig.  4.7.  When  compared  with  the 
phasematching  behavior  of  an  ideal  grating  (solid  black  curve),  at  large  detunings 
from  the  QPM  peak  one  observes  a  QPM  pedestal.  Although  the  data  appear  noisy, 
ripples  in  the  generation  efficiency  are  in  fact  well- resolved  by  the  lock-in  measurement 
technique,  as  demonstrated  by  Inset  2,  which  shows  an  expansion  of  the  data  from 
1555  to  1560  nm.  Similar  ripples  are  observed  in  the  simulations,  and  the  observed 
spectral  features  are  repeatable  over  a  series  of  measurements  of  an  individual  device. 
The  solid  red  line  is  the  average  efficiency  far  from  phasematching:  the  pedestal  height 
(|g(A)|J)  =  5.2  x  10-6,  which,  according  to  Eq.  (4.6),  corresponds  to  a  RDC  error  of 
at  =  8.9%. 

Also  observed  in  Fig.  4.7  is  a  parasitic  QPM  peak  at  1586  nm,  at  a  level  27  dB 
below  the  main  QPM  peak.  We  believe  this  smaller  peak  to  be  due  to  an  interaction 
involving  a  higher-order  spatial  mode  of  the  waveguide.  To  investigate  this  possibility, 
we  used  the  predictions  of  the  waveguide  mode  solver  that  has  been  described  in  [56] . 
For  a  primary  QPM  peak  (involving  the  TMoo  modes  at  both  the  FF  and  SH)  at 
1534  nm,  the  simulations  predict  a  QPM  peak  for  a  sum-frequency  generation  process 
involving  one  TMoo  and  one  TMoi  mode  at  1586  nm,  which  is  coincident  with  the 
observed  parasitic  QPM  feature.  We  therefore  hypothesize  that  a  small  amount  of 
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Figure  4.7:  Measured  relative  SHG  efficiency  (|g(A)|2)  (blue  dashed)  plotted  against 
ideal  sinc2(AA;L/2)  QPM  tuning  curve  (solid  black),  showing  efficiency  pedestal  far 
from  the  QPM  peak.  Inset  1:  measured  SHG  tuning  curve  plotted  on  linear  scale. 
Inset  2:  zoom  of  measured  data  between  1555  and  1560  nm. 

power  was  coupled  into  the  TM0i  mode  at  the  input  of  the  waveguide  due  to  an 
imperfect  mode  filter  or  adiabatic  taper.  Because  the  parasitic  QPM  peak  is  27  dB 
below  the  main  peak,  it  does  not  affect  the  observed  results  for  the  QPM  pedestal 
away  from  the  parasitic  peak:  any  QPM  pedestal  for  this  process  would  also  be  at  a 
level  27  dB  below  the  QPM  pedestal  for  the  fundamental  spatial-mode  interaction. 
We  also  note  that  at  wavelengths  bluer  than  the  primary  QPM  peak  at  1534  nm  we 
observed  several  additional  QPM  peaks  due  to  SHG  of  higher-order  spatial  modes  at 
the  second-harmonic  frequency. 

4.1.4  Domain-disorder-enhanced  parametric  fluorescence 

In  a  quantum  frequency  conversion  device,  the  QPM  pedestal  described  here  can  have 
several  effects.  The  first  is  that,  since  a  strong  pump  at  Xp  must  be  used,  there  will 
be  parasitic  second-harmonic  radiation  created  at  Ap/2.  Since  this  parasitic  process 
produces  radiation  at  a  known  frequency  that  is  generally  far  removed  from  either 
the  signal  or  target  frequencies  it  can  generally  be  filtered  straightforwardly. 
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Figure  4.8:  Frequencies  and  wavelengths  involved  in  downconversion  of  a  weak  signal 
at  A2  =  633  nm  to  a  target  wavelength  Ai  =  1.56  pm.  The  strong  pump  at  Xp  = 
1064  nm  is  also  undergoing  parametric  fluorescence,  creating  noise  when  the  signal 
photons  overlap  spectrally  with  the  detection  bandwidth  near  Ai- 


However,  an  additional  effect  of  the  QPM  disorder  is  in  the  process  of  sponta¬ 
neous  parametric  fluorescence.  Parametric  fluorescence  is  the  spontaneous  splitting 
of  a  pump  photon  into  a  lower  energy  signal/idler  pair  with  energies  usig  +cUidi  =  c op. 
It  is  a  phasematched  process,  and  has  been  studied  in  both  the  microwave  [43]  and 
optical  domains  [81]  since  the  1960s.  The  mathematical  expressions  describing  low- 
gain  parametric  amplification  (a  classical  approximation  to  parametric  fluorescence) 
are  the  same  as  for  SHG  (i.e.  in  both  cases  there  is  no  pump  depletion,  and  the  phase¬ 
matching  behavior  is  governed  by  the  same  Fourier  transform  relation  of  Eq.  (4.3)). 
However,  its  effects  as  a  parasitic  process  are  to  generate  a  “white  noise” -like  collec¬ 
tion  of  photons  at  all  frequencies  00  <  cop  for  which  both  the  signal  and  idler  photons 
are  within  the  transparency  range  of  the  material.  While  the  existence  of  the  noise 
photons  throughout  the  entire  transparency  window  has  not  been  verified,  here  we 
describe  an  experiment  in  which  a  noise  signal  was  observed  in  the  1 . 5- /irn  band  using 
a  pump  at  Ap  =  1.06  /irn,  in  a  waveguide  in  which  this  parametric  fluorescence  process 
was  strongly  phase-mismatched. 

We  investigated  the  parasitic  parametric  fluorescence  noise  in  the  context  of  an 
experiment  designed  to  simulate  the  downconversion  of  single  photons  from  a  di¬ 
amond  NV  center  (with  a  zero-phonon  emission  wavelength  of  A2  =  637  nm;  NV 
center  single-photon  sources  are  described  in  more  detail  in  Sec.  7.1.3)  to  a  wave¬ 
length  Ai  =  1.56  /irn  in  the  telecommunications  band  using  a  pump  wavelength 
Ap  =  (A2  1  —  AhV  =  1064  nm.  A  schematic  of  the  frequencies  that  are  involved  in 
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Figure  4.9:  Experimental  setups  for  measurement  of  weak-signal  downconversion  from 
633  nm  to  1.56  pm.  For  the  measurement  of  SPDC  noise,  we  either  (b)  used  a 
linear-mode  InGaAs  APD  and  a  lock-in  detection  technique,  or,  (c)  measured  the 
upconversion  of  noise  photons  to  633  nm  using  a  Geiger-mode  Si  APD. 


this  process  are  shown  schematically  in  Fig.  4.8.  In  addition  to  pumping  the  con¬ 
version  process  between  the  signal  at  A2  and  the  target  at  Ai,  the  strong  pump  can 
also  undergo  parametric  fluorescence,  indicated  by  the  “QPM  pedestal”.  For  pump 
photons  splitting  into  a  photon  at  the  target  wavelength  Ai,  the  corresponding  idler 
wavelength  is  A;  =  3.3  pm. 

A  schematic  of  our  experimental  setup  is  shown  in  Fig.  4.9.  As  a  weak  signal 
we  used  a  HeNe  laser  at  632.8  nm  in  which  the  power  was  adjusted  using  a  variable 
attenuator.  The  pump  laser  was  a  Nd:YAG  nonplanar  ring  oscillator  (NPRO,  [82]) 
with  maximum  output  power  of  0.75  W.  The  pump  power  was  controlled  using  a 
half-wave  plate  (HWP)  and  polarizing  beam  splitter  (PBS)  and  was  combined  with 
the  signal  beam  on  a  dichroic  mirror  (DM).  The  light  was  in-  and  out-coupled  from 
the  PPLN  waveguide  using  aspheric  lenses  (/  =  8  mm)  with  AR  coating  for  the 
strong  pump.  We  estimate  a  waveguide  coupling  efficiency  of  40%  (5%)  at  the  pump 
(signal)  wavelength.  Overall  coupling  efficiency  for  the  signal  was  not  a  concern  for 
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this  characterization  experiment,  and  as  such  no  special  effort  was  made  to  maximize 
this  quantity.  Additionally,  the  waveguide  mode  filter  for  this  device  was  designed 
to  be  single-mode  for  the  1064-nm  pump;  a  mode  filter  that  was  single-mode  at 
A2  =  633  nm  would  be  in  the  cutoff  regime  at  Ap. 

Using  the  variable  attenuator,  we  adjusted  the  signal  power  such  that  P2(0)  = 
50  11W  was  coupled  into  the  waveguide.  We  then  measured  the  generated  target- 
wavelength  power  at  the  waveguide  output  Pi(L)  as  a  function  of  the  pump  power  Pp 
by  coupling  the  light  at  the  output  of  the  waveguide  into  an  optical  spectrum  ana¬ 
lyzer  (OSA).  By  normalizing  to  the  input  signal  power,  we  calculated  the  conversion 
efficiency  and  plotted  the  results  as  squares  in  Fig.  4.10;  the  solid  curve  is  a  fit  to 
Eq.  (4.1),  showing  good  agreement.  The  maximum  conversion  efficiency  was  found  to 
be  approximately  80%  and  was  reached  at  Pp  =  Pmax  =  78  mW.  At  Pmax  we  observed 
13  dB  (95%)  depletion  of  the  input  signal,  which  we  suspect  was  limited  by  residual 
higher-order  spatial  mode  content  at  A2.  The  conversion  efficiency  was  also  limited 
by  propagation  losses,  which  were  measured  (estimated)  to  be  0.2  (0.1)  dB  cm-1  at 

Ai  (A2). 

Following  measurement  of  the  conversion  efficiency,  we  sought  to  measure  the 
parametric  fluorescence  noise.  Because  of  the  difficulty  of  counting  photons  in  the 
1.55  jum  spectral  region,  this  task  posed  a  measurement  challenge.  We  measured  the 
noise  in  two  ways,  shown  schematically  in  Fig.  4.9(b)  and  (c).  In  short-wavelength- 
pumped  upconversion  single-photon  detectors  pumped  with  1.06  /mi  lasers,  noise 
counts  of  the  level  of  106  counts/s  were  observed,  nominally  within  a  narrow  (0.5  nm) 
nm  spectral  window  around  1.55  /mi  [29].  Therefore,  it  seemed  logical  that  if  we 
deliberately  set  the  collection  bandwidth  substantially  larger  than  this,  we  might 
be  able  to  detect  the  macroscopic  numbers  of  photons  on  a  sensitive  linear-mode 
detector.  In  this  setup,  shown  in  Fig.  4.9(b),  the  polarization  of  the  light  exiting  the 
waveguide  was  first  rotated  from  .s  to  p  using  a  1550-nm  HWP.  The  light  was  then 
sent  though  a  prism  at  Brewster’s  angle  and  sent  through  a  plane-parallel  dielectric 
mirror  with  high  reflectivity  at  1064  nm  (HRP)  to  remove  the  pump  light  from  the 
signal.  It  was  then  coupled  into  a  single-mode  optical  fiber  and  sent  to  a  linear-mode 
InGaAs  APD.  To  achieve  high  sensitivity,  the  pump  was  modulated  at  ~300  Hz  using 
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Figure  4.10:  Measured  DFG  conversion  efficiency  (squares)  and  fit  to  Eq.  (4.1)  (solid 
blue),  showing  maximum  conversion  at  Pmax  =  78  mW. 


an  optical  chopper,  and  the  photocurrent  of  the  InGaAs  APD  was  detected  using  a 
lock-in  amplifier. 

We  calibrated  the  measurement  bandwidth  by  sending  the  light  from  a  tunable 
C-band  ECDL  through  the  waveguide  and  into  the  measurement  setup:  the  collection 
bandwidth  was  measured  to  be  44  nm  (1/e2  full-width).  An  absolute  calibration  of 
the  responsivity  of  the  detection  setup  was  done  by  inserting  a  1.56- /mi  probe  beam 
of  known  optical  power  into  the  filtering  setup.  Plotted  as  dots  in  Fig.  4.10  are  data 
recorded  by  monitoring  the  lock-in  signal  on  a  computer  while  the  pump  power  was 
varied  using  the  HWP.  The  data  are  well  fit  by  a  linear  dependence  on  pump  power  as 
expected  for  a  spontaneous  scattering  process.  We  find  a  noise-photon  spectral  density 
(NPSD)  of  1.45  x  106  s_1  nm-1  at  Pp  =  Pmax.  We  note  that  any  noise  produced  in 
waveguides  without  a  QPM  grating  was  below  the  estimated  105  s^1  nm-1  sensitivity 
of  the  detection  system. 

An  expression  for  the  generation  rate  of  photon  pairs  via  SPDC  in  a  waveguide 
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has  been  given  by  Fiorentino  and  coworkers  [83].  They  find 

d JV,  Pr  2  ,(AkL\ 

dAi  eoniniUpXi A*  Ai  \  2  J 


(4.7) 


where  Aj  is  the  interaction  area  given  as  Aj  =  @~2.  Evaluating  this  expression  for  the 
expected  NPSD  for  our  experimental  configuration  with  A kL  «  —1.8  x  104  gives  a 
rate  of  approximately  102  photons  s~4  nnr1,  about  4  orders  of  magnitude  lower  than 
the  value  observed  in  the  experiment.  We  note  here  that  additional  complexities  arise 
in  our  experimental  configuration  due  to  the  fact  that  the  idler  wave  at  3.3  pm  is  in 
the  cutoff  regime  of  the  waveguide,  such  that  radiation  modes  must  be  considered.  As 
a  result,  Ai  for  this  process  cannot  be  calculated  in  the  way  described  in  Chapter  2 
for  interactions  involving  only  bound  modes  of  the  structure.  A  detailed  discussion  of 
this  effect  is  given  in  Appendix  B:  using  the  waveguide  mode  completeness  relations, 
when  there  are  no  bound  modes  at  the  idler  wavelength,  the  mode-overlap  integral  0 
takes  a  different  form  than  in  Eq.  (2.10).  We  show  in  Appendix  B  that  in  the  absence 
of  bound  modes  at  the  idler  frequency, 

@cer  =  =  JJ  \Mx,y)Mx,y)\2  dx  (4.8) 

As  we  have  seen  above,  the  QPM  pedestal  dominates  the  efficiency  of  a  process  when 
|  AkL\  >  1007T.  We  then  replace  the  sine2  (AkL/2)  in  Eq.  (4.7)  with  the  QPM  pedestal 
height 

»L  =  (1 9(A*  -►  oo)|2)  =  4  (l  -  (4.9) 

and  calculate  the  NPSD  for  an  observed  grating  disorder,  obtaining  a  value  cq  ~  21%. 
This  value  is  higher  than  that  observed  by  direct  imaging  of  a  comparably  fabricated 
device,  but  we  have  not  investigated  the  uniformity  of  cq  across  an  ensemble  of  devices. 
Unfortunately  due  to  the  destructive  nature  of  RDC  error  measurements  using  Zygo 
interferometry  we  were  unable  to  measure  cq  directly  for  this  device. 

Any  noise  photons  generated  via  parametric  fluorescence  within  the  acceptance 
bandwidth  5 A  around  the  target  wavelength  of  Ai  of  the  device  are  indistinguishable 
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Figure  4.11:  Measurement  of  upconverted  parametric  fluorescence  noise  (circles)  and 
fit  to  semiclassical  theory  (solid). 

from  target  photons  and  will  be  efficiently  upconverted  via  SFG  to  the  signal  wave¬ 
length  at  A2.  We  can  therefore  measure  the  effects  of  the  parametric  fluorescence  noise 
via  measuring  this  upconverted  noise  signal  at  633  nm  using  a  Si  photon  counter.  The 
experimental  setup  we  used  for  this  measurement  is  shown  in  Fig.  4.9(c).  The  light 
exiting  the  waveguide  was  passed  through  a  HWP  for  633  nm  to  rotate  the  polariza¬ 
tion  from  s  to  p,  and  was  then  sent  through  a  Brewster-angled  prism  to  separate  the 
upconverted  noise  photons  at  633  nm  from  the  parasitic  radiation  at  532  nm  due  to 
SHG  of  the  pump  (which  had  an  optical  power  of  approximately  50  11W  due  to  RDC- 
error-enhanced  SHG).  Following  the  prism,  an  additional  long-pass  filter  was  used 
and  the  light  was  then  focused  onto  a  Si  APD  (PerkinElmer  SPCM-AQR-14).  Our 
measurement  of  the  upconverted  noise  rate  is  shown  as  dots  in  Fig.  4.11.  It  can  be 
seen  that  the  upconverted  noise  grows  quadratically  at  first,  followed  by  linear  growth 
with  pump  power.  This  behavior  can  be  understood  by  a  simple  physical  argument: 
at  low  pump  powers,  both  the  noise-photon  generation  rate  and  the  upconversion 
efficiency  are  linear  in  pump  power,  leading  to  quadratic  scaling.  At  higher  pump 
powers,  the  noise  photon  generation  rate  remains  linear  but  the  conversion  efficiency 
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flattens,  saturates,  and  rolls  over,  such  that  the  total  upconversion  rate  is  effectively 
linear. 

We  have  constructed  a  semiclassical  model  for  this  parasitic  upconversion  process 
in  QPM  media  with  RDC  errors.  The  coupled  wave  equations  for  the  combined  SPDC 
and  SFG  processes  are  given  as: 

=  iKNa*lapelAkNZ  (4.10a) 

az 

=  ftia*a2e* Aksz  +  inNa*apelAkNZ  (4.10b) 

=  iKSapaie~'lAksz  (4.10c) 

where  and  Ak's  are  the  coupling  coefficient  and  phase  mismatch  for  the  SFG 
process,  and  kn  and  Ak'N  are  for  the  noise  (parametric  fluorescence)  process.  To 
study  these  combined  interactions  in  a  material  with  RDC  errors,  our  approach  is 
to  numerically  construct  a  grating  with  quasi-random  duty-cycle  errors,  and  solve 
Eqns.  (4.10),  and  take  an  ensemble  average  over  many  such  disordered  gratings.  The 
boundary  conditions  used  to  simulate  Eqs.  (4.10)  were  a  small  seed  a8(0)  =  cpo  and 
ai(0)  =  0.  The  results  of  this  simulation  are  shown  as  the  solid  curve  in  Fig.  4.11. 
The  simulated  curve  has  one  free  parameter,  an  overall  scale  factor,  but  the  particular 
shape  has  been  determined  by  the  value  Pmax  =  78  mW. 

4.1.5  RDC  error  tolerances 

We  can  combine  the  rate  of  SPDC  of  Eq.  (4.7)  and  the  QPM  pedestal  height  of 
Eq.  (4.9)  to  derive  tolerances  on  the  maximum  allowable  random  duty-cycle  error 
to  achieve  a  desired  conversion  efficiency  t]  and  NPSD  dNi/d\\.  Writing  r/  as  a 
function  of  Pp  using  Eq.  (4.1)  gives 

Pp  = (sin-1^)2.  (4.11) 

'/nor-^ 

Inserting  this  result  into  the  NPSD  expression  of  Eq.  (4.7)  and  rearranging,  we  are 
led  to  the  following  result  following  a  Taylor  series  expansion  (to  second  order  in  c?i/l) 
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Figure  4.12:  Calculated  random  duty-cycle  error  di  =  ai/l  tolerances  for  achieving 
conversion  efficiency  ij  and  noise-photon  spectral  density  dN\/d\\  for  a  QFC  device 
with  domain  disorder. 

for  the  pedestal  height  from  Eq.  (4.9): 


2Nrii\i\\ 

7r2n2A2c 


(4.12) 


Here,  N  =  L/{Kq/2)  is  the  number  of  QPM  domains  in  the  waveguide,  and  r  is 
a  ratio  of  the  square  of  the  mode-overlap  integrals  for  the  SPDC  to  the  SFG/DFG 
processes.  For  a  waveguide  with  bound  modes  at  all  frequencies,  0  was  given  in 
Eq.  (2.10),  and  for  a  waveguide  with  a  radiated  idler,  the  calculation  in  Appendix  B 
gives  the  form  of  Eq.  (4.8).  For  our  experimental  parameters,  we  calculate  r  =  1.03, 
a  result  based  on  simulated  waveguide  modes  for  the  waveguide  parameters  used  in 
the  device  fabrication.  By  use  of  the  Cauchy-Schwarz  inequality  it  can  be  noticed 
that  in  the  Cerenkov  configuration  r  >  1  in  all  cases. 

We  can  visualize  the  implications  of  Eq.  (4.12)  by  making  a  contour  plot  of  values 
of  constant  di  as  a  function  of  the  conversion  efficiency  g  and  NPSD  dNi  /  dAi,  shown 
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in  Fig.  4.12  for  a  device  with  N  =  104  (corresponding  to  a  device  of  length  L  —  5  cm 
with  a  poling  period  =  10  pm)  and  operating  at  the  wavelengths  shown  in  Fig.  4.8. 
The  ideal  QFC  device  would  have  as  high  a  conversion  efficiency  and  as  low  a  NPSD 
as  possible,  corresponding  to  a  position  at  the  lower  right  corner  of  Fig.  4.12.  As 
can  be  seen  from  Eq.  (4.12),  operation  in  this  regime  requires  extraordinarily  tight 
tolerances  on  RDC  errors  during  fabrication  <q  <  0.5  %,  beyond  current  fabrication 
technology.  The  “state-of-the-art”  of  our  electric  held  poling  procedures  as  described 
in  Sec.  3.2  seems  to  yield  <q  ~  10%,  which  corresponds  to  a  noise-photon  spectral 
density  of  approximately  105  photons  s-1  nm_1,  a  value  that  seems  prohibitively 
high  for  use  in  applications  even  if  very  aggressive  optical  filtering  is  used.  Another 
consequence  of  the  result  of  Eq.  (4.12)  is  that  the  tolerances  on  cq  are  reduced  as  the 
device  gets  longer  and  N  increases.  This  behavior  results  from  the  increase  as  77  oc  L2 
in  the  conversion  efficiency  of  the  QFC  process  while  the  noise  generation  rate  (a 
spontaneous  process)  grows  as  d Ni  /  dAi  oc  L,  such  that  longer  L  is  always  favored. 

We  have  applied  the  analyses  described  here  to  a  previous  demonstration  of  1064- 
nm-pumped  quantum  frequency  conversion:  an  upconversion  single-photon  detector 
demonstrated  by  Albota  and  Wong  at  MIT  in  2004  [29].  That  experiment  used  a 
bulk  PPLN  crystal  as  the  conversion  medium,  requiring  a  redefinition  of  r  to  include 
the  effects  of  Gaussian  beam  propagation.  Working  in  the  near-held  Gaussian  ap¬ 
proximation  (in  which  the  diffraction  of  a  beam  through  the  crystal  is  assumed  to 
be  negligible)  one  finds  that  r  can  be  replaced  by  a  ratio  of  the  g-factors  used  to 
describe  the  overlap  of  a  held  with  its  driving  polarization  [84] .  We  find  that  a  value 
a  1  =  28%  would  explain  the  observed  noise  level  in  that  experiment.  While  this  value 
of  cq  is  large,  it  is  also  consistent  with  the  value  one  would  calculate  to  account  for 
the  reduced  conversion  efficiency  reported  in  Ref.  [29]  with  respect  to  theory. 

4.2  Spontaneous  Raman  scattering 

Raman  scattering  involves  the  scattering  of  a  photon  oh  an  optical  phonon  in  a  mate¬ 
rial  to  either  redshift  the  photon  by  creating  a  phonon  (in  so-called  Stokes  scattering) 
or  blueshifting  by  absorbing  a  phonon  (anti-Stokes  scattering).  Fig.  4.13  diagrams 
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Stokes  anti-Stokes 


Figure  4.13:  Diagram  of  spontaneous  Raman  scattering  in  either  the  Stokes  (in  which 
a  photon  at  frequency  cu  redshifts  to  u  —  Aco  via  the  creation  of  a  phonon  of  frequency 
Acn)  or  anti-Stokes  configuration  (in  which  a  photon  is  blueshifted  to  cu  +  Au  via  the 
annihilation  of  a  phonon). 

each  of  these  two  processes.  Raman  scattering  was  the  first  inelastic  scattering  pro¬ 
cess  observed  with  photons;  it  was  observed  by  C.  V.  Raman  in  1928  by  sending  a 
beam  of  focused  sunlight  through  a  liquid  and  observing  the  scattering  by  eye  [85]. 
Since  then,  Raman  scattering  has  been  used  extensively  as  a  material  characterization 
method  as  the  phonon  frequencies  and  polarization  response  gives  a  fingerprint  of  the 
material. 

Raman  scattering  in  crystals  has  been  studied  by  Loudon  [86] .  Raman  scattering 
in  LiNb03  appears  to  have  been  first  observed  by  Schaufcle  and  Weber  in  1966  [87]. 
Lithium  niobate  is  a  trigonal  crystal  with  3m  symmetry:  the  conventional  hexagonal 
unit  cell  contains  6  formula  units  [88].  The  Raman  spectrum  has  been  measured 
many  times,  with  a  focus  put  on  identifying  and  characterizing  the  peaks  of  the 
Raman  spectrum.  However,  in  quantum  frequency  converters,  it  seems  obvious  that 
Raman  peaks  should  be  avoided.  To  our  knowledge,  there  has  not  been  a  detailed 
investigation  of  the  character  of  SRS  both  between  and  at  frequency  shifts  past  the 
Raman  peaks,  which  we  describe  in  the  next  section. 
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4.2.1  Raman  spectroscopy  of  H+:LiNbC>3 

We  sought  a  detailed  characterization  of  the  Raman  scattering  behavior  of  LiNbC>3 
for  frequency  shifts  away  from  any  of  the  dominant  phonon  peaks.  In  a  Raman 
spectroscopy  experiment,  the  basic  quantity  measured  is  the  differential  scattering 
cross  section  dcr  /  dfh  A  Raman  spectrometer  records  a  spectrum  S(Af),  which 
is  a  number  of  photocounts  at  a  frequency  difference  A /  =  /0  —  fs,  where  /o  is 
the  frequency  of  the  laser  excitation  and  fs  is  the  Stokes  frequency.  The  measured 
spectrum  and  the  differential  cross  section,  which  can  be  used  to  calculate  the  gain 
in  a  stimulated  Raman  scattering  experiment,  are  related  empirically  as 


(4.13) 


where  ps  is  the  number  density  of  scatterers,  Lg  is  the  length  of  the  scattering  vol¬ 
ume,  Nq  is  the  number  of  incident  laser  photons  over  the  course  of  the  measurement, 
and  AD  is  the  solid  angle  subtended  by  the  collection  optics.  These  quantities  (other 
than  ps  which  is  a  material  property)  depend  in  detail  on  the  experimental  conditions, 
and  can  be  hard  to  calibrate  accurately.  We  are  primarily  interested  in  the  lineshape 
function  gjJ{Af)1  which  can  be  used  straightforwardly  to  compare  the  Raman  scat¬ 
tering  rate  at  different  frequency  shifts  A/.  We  have  performed  Raman  spectroscopy 
experiments  on  a  variety  of  samples  using  a  WiTec  Alpha300  S  Raman  microscope. 
The  use  of  a  Raman  micrsocope  rather  than  a  bulk  Raman  microscopy  setup  enabled 
the  investigation  of  the  effects  of  proton  doping,  which  it  was  thought  may  alter  the 
Raman  spectrum  substantially  [89]. 

Initial  measurements  were  of  an  A-cut  LiNbO.j  sample.  We  measured  a  Raman 
spectrum  in  the  x(zz)x  configuration1,  the  results  of  which  are  shown,  in  logarithmic 
units,  in  Fig.  4.14.  Two  measurements  were  taken  using  a  grating  with  1800  lines/mm 
over  overlapping  frequency  ranges,  and  the  background  counts  from  the  Peltier-cooled 
CCD  array  detector  (which  were  measured  to  be  time  invariant  to  high  statistical 

lrThe  notation  a(bc)d  is  commonly  used  in  the  Raman  scattering  literature:  a  is  excitation  di¬ 
rection,  b  and  c  are  the  polarizations  of  the  excitation  beam  and  collected  signal,  and  d  is  the 
propagation  direction  of  the  collected  signal.  All  quantities  are  defined  with  respect  to  the  crystal 


axes. 
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Figure  4.14:  Measured  Raman  spectrum  of  congruent  X-cut  LiNbOa  in  the  x{zz)x 
configuration,  and  a  fit  to  a  sum  of  Lorentzians  as  in  Eq.  (4.14),  where  the  fit  param¬ 
eters  are  provided  in  [90]. 


accuracy)  were  subtracted.  Raman  measurements  are  most  frequently  shown  on  linear 
scales  because  Raman  spectroscopists  have  been  primarily  interested  in  the  peak 
locations.  The  logarithmic  scaling  aids  in  a  visual  understanding  of  the  unexpectedly 
slow  rate  at  which  the  spectrum  rolls  of  at  large  frequency  shifts. 

The  spectral  response  of  an  individual  Raman  scattering  mode  is  well  described 
by  a  Lorentzian  distribution  [12].  It  is  generally  assumed  that  materials  with  many 
Raman-active  phonons  such  as  crystals  can  be  well- fit  by  a  sum  of  Lorentzians  [86] . 
We  have  fit  the  complex  Raman  susceptibility  inferred  from  the  measured  spectrum 
via  the  Kramers-Kronig  relation  to  such  a  sum  of  Lorentzians 


XrW)  = 


Nl 

3=0  D 


(A/)2  +  2z7iA  / 


(4.14) 


where  a  table  of  fit  parameters  can  be  found  in  our  publication  [91].  In  addition  to 
the  peaks  referenced  in  [91],  we  insert  the  peak  at  885  cm-1  “by  hand”  as  the  fitting 
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routines  did  not  converge  properly  when  it  was  included.  In  the  vicinity  of  the  peaks, 
for  frequency  shifts  smaller  than  approximately  A /  <  1000  cm-1,  the  superposition 
of  Lorentzians  accurately  describes  the  Raman  scattering  behavior.  However,  for 
large  frequency  shifts  A /  >  1000  cm-1,  the  measured  spectrum  S(Af)  exceeds  the 
Lorentzian  predicition  significantly,  in  an  effect  we  have  termed  a  high-frequency 
“Raman  pedestal,”  in  analogy  with  the  QPM  pedestal  described  in  Sec.  4.1. 

We  next  sought  to  determine  the  origin  of  the  Raman  pedestal.  To  the  best  of  our 
knowledge,  this  large-shift  pedestal  has  not  been  analyzed  elsewhere  in  the  literature 
on  LiNb03.  We  initially  suspected  it  was  due  to  the  high  defect  concentration  of 
congruently- melting  LiNb03.  When  grown  from  a  melt,  LiNb03  does  not  contain 
a  1:1  stoichiometric  ratio  of  Li:Nb,  but  rather  has  a  concentration  ratio  [Li] /([Li]  + 
[Nb])  =  0.485.  This  non-stoichiometry  results  in  a  high  degree  of  crystal  disorder  due 
to  the  large  number  of  point  defects  [92],  By  comparing  a  congruent  crystal  to  one 
with  stoichiometric  composition,  it  is  therefore  possible  to  test  this  hypothesis.  Since 
stoichiometric  LiNb03  is  not  readily  available,  we  compared  congruent-composition 
LiTa03  against  stoichiometric  LiTa03.  The  stoichiometric  LiTa03  showed  a  similar 
shelf  as  congruent  LiTa03,  suggesting  that  disorder  is  not  important  to  the  existence 
of  the  shelf.  The  effects  of  stoichiometric  versus  congruent  composition  were  found  to 
be  Raman  line-narrowing,  which  has  been  described  in  the  literature  [93,  94],  without 
the  addition  of  a  pedestal. 

We  also  suspected  the  pedestal  may  be  due  to  the  proton  exchange  process.  As 
discussed  in  Chapter  3,  proton  exchange  removes  Li+  ions  and  replaces  them  with 
H+  ions,  raising  the  refractive  index  for  e-polarized  light.  The  Raman  scattering 
strength  of  H+:LiNb03  has  been  found  to  be  substantially  diminished,  which  partially 
explains  the  so-called  “dead  layer”  near  the  surface  of  a  PE  waveguide  in  which  x('2> 
vanishes  [89].  We  characterized  both  unannealed  and  annealed  PE  planar  waveguides 
formed  on  both  X-  and  Z- cut  substrates,  as  well  as  RPE  planar  waveguides  on  Z- cut 
substrates.  While  we  have  not  endeavored  to  obtain  a  detailed  picture  of  the  effects 
of  proton  doping  on  Raman  behavior,  we  did  not  see  any  qualitative  differences  in  the 
large-frequency-shift  behavior  that  would  indicate  that  proton  exchange  is  the  cause 
of  the  Raman  pedestal. 


CHAPTER  4.  NOISE  PROCESSES  IN  QFC  DEVICES 


57 


In  the  absence  of  independent  confirmation  of  the  Raman  pedestal  by  direct  spec¬ 
troscopic  measurement,  measurements  of  noise  behavior  in  a  quantum  frequency  con¬ 
version  device  can  provide  indirect  insight  on  this  effect.  The  theory  of  noise  gen¬ 
eration  by  spontaneous  Raman  scattering  is  described  in  the  next  section  and  an 
experiment  studying  the  noise  behavior  of  a  long-wavelength-pumped  upconversion 
single-photon  detector  for  the  1.55-pm  band  is  described  in  Sec  5.3. 


4.2.2  Raman  noise  in  QFC  devices 

Here,  we  seek  to  describe  the  physics  of  spontaneous  Raman  scattering  in  a  wave¬ 
guide  geometry  that  can  be  used  to  predict  the  spectral  dependence  of  noise  in  a 
quantum  frequency  converter  in  which  the  strong  pump  is  undergoing  Raman  scat¬ 
tering.  Describing  the  Raman  scattering  process  using  a  formalism  that  includes  a 
Raman  susceptibility  \R  has  been  presented  in  [12].  We  modify  that  presentation 
slightly  to  include  the  effects  an  interaction  involving  only  a  single  spatial  mode  at 
each  frequency  in  a  more  restrictive  geometry.  An  equation  for  the  evolution  of  pho¬ 
ton  held  envelopes  (defined  as  in  Chapter  2)  with  a  laser  frequency  co0  and  scattered 
frequency  cos,  where  Au  =  cos  —  co0,  is: 


do* 

dz 


,3ujsZ0xr(Auj) 
4 nsnLc 


P0@Rh(Au)as, 


(4.15) 


where  we  have  written  the  evolution  equation  in  terms  of  the  pump  power  P0  = 
hu;o|oo|2>  the  Raman  overlap  integral 


Qr  = 


l^o(^,  2/)|2|^s(a;,  2/)|-  dx  d y, 


(4.16) 


and  h(Au>)  is  a  factor  that  describes  the  occupation  of  the  relevant  phonon  modes 
and  has  different  forms  for  Stokes  and  anti-Stokes  scattering. 

A  proper  quantum  mechanical  treatment  of  Raman  scattering  (in  which  both  the 
light  held  and  vibrations  are  quantized)  is  required  to  deduce  the  form  of  h(Aco)  and 
will  not  be  presented  here  [95].  The  quantity  h(Au)  takes  into  account  the  Raman 
signals  by  summing  across  all  the  vibrational  energy  levels  of  the  oscillating  modes 
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in  the  system.  For  Stokes  and  anti-Stokes  scattering,  respectively,  we  find 

h( Auj  <  0)  =  1  +  (n„(|Ao;|))  (4.17a) 

h(Auj  >  0)  =  (n„(|Acn|))  (4.17b) 

where  the  quantity  (nv  ( |  Aw  | ) )  is  the  average  excitation  of  a  harmonic  oscillator  and 
is  obtained  via  the  Boltzmann  distribution  as 


(n„(|Au;|))  = 


eh\Auj\/kT  _ 


(4.18) 


where  T  is  the  temperature  of  the  system. 

To  calculate  the  spontaneous  scattering  rate,  we  define  Ns  =  |as|2  and  notice 


d Ns  _  „da.  da* 

Ah  ~  dz7  +“*  dT' 


(4.19) 


Combining  Eq.  (4.15)  and  (4.19),  and  inserting  an  “extra  photon”  in  the  usual  way 
to  describe  spontaneous  emission  processes,  we  are  led  to  a  photon-number  evolution 
equation 


diV, 


3ojsZ0 


XR,i(Auj)h(Auj)PoOR(Ns  +  1), 


(4.20) 


dz  2nsnoc/ 

where  Xru  is  the  imaginary  component  of  the  full  complex  Raman  susceptibility  xr- 
To  calculate  the  spontaneous  scattering  rate,  we  assume  77,(0)  <C  1  and  find,  after 
propagation  through  a  device  of  length  L,  a  noise-photon  spectral  density 


dNs  _  67r2cZoPoh(Au)0RXR,i(Auj)L 


dA 


A  3snsn0 


(4.21) 


Evaluating  the  spontaneous  scattering  rate  dNs  /  dA  as  a  function  of  Auj  for  a  pump 
power  Pp  =  100  mW  (assuming  a  pump  wavelength  As  =  1.55  pm)  and  a  device  length 
L  —  5  cm  gives  the  results  shown  in  Fig.  4.15.  The  curves  in  Fig.  4.15  are  calculated 
using  the  measured  Raman  spectra  described  in  Sec.  4.2.1  for  large  frequency  shifts. 
We  were  unable  to  measure  the  Raman  scattering  signal  for  frequency  shifts  smaller 
than  approximately  200  cm-1  due  to  the  limitations  of  the  notch  filter  used  in  the 
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Pump  wavelength  kp  for  \  =  1 .55  pm 
1.3  1.4  1.5  P  1.6  S  1.7  1.8  1.9  2 


Frequency  shift  (co1  -  cop)/(2tt) 

Figure  4.15:  Spontaneous  Raman  scattering  noise  photon  spectral  density  calculated 
from  measured  Raman  spectrum,  plotted  against  signal-pump  frequency  shift  Aoj  = 
uji  —  up  and  equivalently  against  the  pump  wavelength  Xp  used  for  an  input  signal 
wavelength  Ai  in  an  upconverter.  Results  calculated  from  our  Raman  spectroscopy 
measurements  are  plotted  in  blue,  and  the  low-frequency-shift  data  plotted  in  green 
are  taken  from  Ref.  [96]. 


Raman  microscope.  The  low-frequency-shift  data  of  Fig.  4.15,  plotted  in  green,  have 
been  copied  from  the  work  of  Surovtsev  et  al.  and  have  been  adjusted  to  give  the 
same  Raman  signal  for  the  252-cm_1  peak  [96]. 

The  data  of  Fig.  4.15  clearly  show  major  difference  between  Stokes  (Aca  <  0) 
and  anti-Stokes  (Acu  >  0)  scattering.  As  the  pump  frequency  becomes  substantially 
lower  than  the  signal  (corresponding  to  Ac <j  increasing),  or  correspondingly  a  very 
long  pump  wavelength,  anti-Stokes  Raman  scattering  decreases  dramatically.  The 
effects  of  anti-Stokes  Raman  scattering  on  an  upconversion  single-photon  detector  for 
1.55-/im-band  signals  will  be  discussed  in  Chapter  5. 
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4.3  Design  implications 

In  the  previous  two  sections  we  have  discussed  the  two  primary  noise  sources  in  quan¬ 
tum  frequency  conversion  devices.  To  the  best  of  our  knowledge,  no  other  spontaneous 
scattering  processes  have  been  found  to  play  a  major  role  in  the  noise  performance 
of  QFC  devices.  Typical  QPM-grating  disorder  parameters  ai/l  ~  10%  yield  NPSDs 
that  are  in  general  too  high  for  useful  QFC  systems.  This  fact  alone  basically  elim¬ 
inates  the  possibility  of  using  pump  wavelengths  Xp  shorter  than  the  input  signal 
Ai .  The  analysis  of  SRS  from  the  previous  section  shows  that  Raman  scattering  is 
a  major  noise  source  (in  the  Stokes  configuration)  even  for  pump-signal  frequency 
differences  well  past  the  main  Raman  peak  at  632  cm-1.  Anti-Stokes  scattering  is 
found  to  be  a  significant  noise  source  for  0  <  Aoj/2tt  <  1000  cm-1.  However,  due  to 
its  dependence  on  the  Boltzmann  ratio,  anti-Stokes  scattering  is  reduced  significantly 
for  large  positive  Aw. 

Thus,  in  most  situations,  very  low  noise  performance  will  only  be  achieved  when 
Acu  >  1500  cm-1.  Such  a  system,  for  downconversion  of  single-photons  from  quan¬ 
tum  dots,  will  be  described  in  Chapter  6.  Another  approach  for  achieving  low-noise 
performance,  however,  is  optical  filtering.  In  most  circumstances,  an  input  signal  will 
be  spectrally  narrow  compared  to  either  the  acceptance  bandwidth  of  the  converter 
or  any  “naive”  optical  filtering  at  the  output.  We  describe  one  experiment  in  Sec.  7.2 
in  which,  through  using  an  optical  filter  substantially  narrower  than  the  upconversion 
acceptance  bandwidth,  the  noise  can  be  reduced  significantly  (by  a  factor  of  about 
50)  compared  with  the  naive  filtering  case. 

We  also  note  that  for  spontaneous  scattering  processes,  the  signal-to-noise  ratio  is 
generally  improved  with  longer  devices.  As  the  conversion  efficiency  of  a  device  scales 
as  L 2  while  the  generation  of  noise  photons  scales  as  L ,  this  behavior  is  understood 
quite  simply.  We  have  not  explored  the  use  of  devices  longer  than  approximately  6.7- 
cm,  roughly  the  maximum  chip  length  available  by  processing  76-mm  (3”)  wafers.  In 
future  work,  the  improved  efficiency  from  fabricating  on  4”  wafers  may  enable  further 
reduction  of  noise. 


Chapter  5 

Upconversion-assisted 
single-photon  detectors 


Single-photon  detectors  are  important  in  a  variety  of  applications,  ranging  from  bio¬ 
chemistry  to  sensing  to  communications.  This  chapter  presents  work  on  applications 
of  quantum  frequency  converters  towards  improving  single-photon  detectors  based 
on  avalanche  photodiodes  (APDs).  Our  work  on  a  long- wavelength  pump  source  for 
upconversion  detectors  (based  on  a  monolithic  PPLN  optical  parametric  oscillator 
(OPO))  was  published  in  [97],  and  the  performance  of  a  1.55-/im-band  upconversion 
detector  was  published  in  [67]. 

5.1  Single-photon  detection  with  avalanche  photo¬ 
diodes 

The  objective  of  a  single-photon  detector  (SPD)  is  to  produce  a  macroscopic  electrical 
pulse  upon  the  absorption  of  a  single  optical  photon.  Photons  of  visible  or  near-IR 
frequency  have  energies  in  the  range  of  10~19  J  (0.1  aJ),  so  the  devices  must  be 
extremely  sensitive  to  incoming  optical  radiation.  A  recent  review  due  to  Hadfield 
describes  the  state-of-the  art  for  SPDs  across  a  wide  range  of  device  technologies  [22], 
The  three  most  important  performance  parameters  of  single  photon  detectors 
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are  the  photon  detection  efficiency  (PDE)  pu,  the  dark-count  rate  (DCR)  D,  and 
the  timing  jitter  At.  The  PDE  is  simply  the  probability  that  an  incoming  single 
photon  will  register  an  output  pulse,  which  is  a  function  of  the  wavelength  of  the 
incoming  photon.  The  DCR  is  the  rate  at  which  the  detector  produces  output  pulses 
when  it  is  known  that  there  are  no  incident  photons.  The  timing  jitter  describes 
the  spread  of  arrival  times  of  the  electrical  pulses  for  a  train  of  well-timed  optical 
pulses.  Additionally,  the  spectral  range  of  a  SPD  is  of  critical  importance  as  well, 
as  it  determines  which  wavelengths  can  be  used  in  system  applications.  Two  other 
important  parameters,  which  we  will  not  discuss  here,  are  the  ability  to  resolve  photon 
number  (which  can  be  important  for  linear-optical  quantum  computation  [98]),  and 
the  dead  time  td  (which  will  be  discussed  in  Chapter  7)  which  determines  the  maximal 
flux  of  incident  photons  before  the  detector  saturates. 

The  first  widely  used  SPDs  were  based  on  photomultiplier  tubes,  which  are  still 
the  primary  detectors  used  for  the  UV  spectral  range.  The  most  established  detectors 
for  the  visible  and  near-IR  are  based  on  semiconductor  avalanche  photodiodes.  Addi¬ 
tional  emerging  detector  technologies  are  based  on  either  superconducting  nanowires 
[99]  or  superconducting  transition-edge  sensors  [100],  but  both  of  these  technologies 
require  cryogens  and  are  viewed  to  be  somewhat  impractical  for  field  use  in,  for  ex¬ 
ample,  quantum  key  distribution  systems.  Here,  we  focus  on  detection  with  APDs. 

Avalanche  photodiode  SPDs  consist  of  a  semiconductor  p-n  (or  p-i-n )  junction 
reverse-biased  above  the  breakdown  voltage,  in  so-called  “Geiger  mode”  operation. 
Photocarriers  produced  by  the  absorption  of  a  photon  undergo  avalanche  gain  trigger¬ 
ing  macroscopic  breakdown  of  the  junction.  Following  the  avalanche,  the  current  must 
be  stopped  either  through  the  use  of  active  or  passive  quenching  schemes  [101,  102], 
There  are  two  commercial  device  types  based  on  Si,  which  are  referred  to  as  “thick- 
junction”  or  “thin-j unction”  types.  The  thick-junction  detectors  are  optimized  for  the 
highest  PDE  and  have  a  thick  active  region,  which  results  in  a  large  timing  jitter  due 
to  variability  in  the  carrier  extraction  time  depending  on  where  in  the  junction  the 
absorption  event  took  place  [103].  The  performance  of  these  detectors  with  respect 
to  recent  reports  on  InGaAs  APDs  is  summarized  in  Table  5.1. 

The  primary  drawback  with  InGaAs  APDs  is  the  effect  of  afterpulsing.  After  the 
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Tabic  5.1:  Summary  of  avalanche  photodiode  single-photon  detector  performance  for 
both  the  Si  and  InGaAs  material  systems. 


Si  (thick) 

Si  (thin) 

InGaAs  (Self- 
differencing  [104]) 

InGaAs  (Free- 
running  [105]) 

Peak  PDE  (%) 

70 

49 

10 

4.5 

Spectral  range  (/mi) 

0. 4-1.0 

0.4-0. 8 

1. 1-1.6 

1. 1-1.6 

Dark-count  rate  (s_1) 

100 

20 

1.6  x  104 

5  x  104 

Timing  jitter  At  (ps) 

400 

35 

55 

450 

Max  rate  ry,1  (MHz) 

10 

10 

100 

1.1 

registration  of  an  avalanche  event,  some  portion  of  the  photocarriers  occupy  intra¬ 
band  trap  states,  and  the  decay  of  these  trap  states  can  excite  additional  avalanches 
[106].  The  primary  way  around  this  effect  is  to  gate  the  detector;  that  is,  to  bias  in 
Geiger  mode  only  for  a  short  duration  (typically  0.5  to  10  ns),  followed  by  turning  off 
the  detector  such  that  afterpulses  cannot  undergo  avalanche  gain.  Researchers  have 
developed  sophisticated  electronic  gating  and  quenching  schemes,  including  sine-wave 
biasing  [107]  and  self-differencing  [104]  which  enable  high  clock  rates  and  improved 
performance  over  the  initial  generation  of  devices.  As  material  quality  has  improved, 
free-running  InGaAs  APDs  have  been  developed,  but,  as  is  seen  in  Table  5.1,  the 
DCR  is  still  very  high. 

The  idea  of  upconversion  detection  is  that  if  a  single  photon  with  an  energy 
lower  than  the  Si  bandgap  (wavelength  longer  than  1  /im)  can  be  upconverted  with 
high  efficiency  to  an  energy  larger  than  the  Si  bandgap,  one  could  then  combine 
the  efficient  upconverter  with  a  Si  APD  for  improved  single-photon  detection.  A 
summary  of  experimental  achievements  in  upconversion  SP  detection  has  been  given 
in  Chapter  2. 


5.2  Long- wavelength  pump  source  based  on  mono¬ 

lithic  PPLN  OPO 

Upconversion  single-photon  detection  is  only  a  viable  technology  if  the  upconversion 
can  be  done  with  little  added  noise.  Therefore,  as  discussed  in  Chapter  4,  a  pump 
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Figure  5.1:  Schematic  diagram  of  monolithic  PPLN  OPO  used  for  long-wavelength- 
purnped  1.5-pm  upconversion  single  photon  detector.  The  OPO  was  pumped  at 
1064  nm  by  a  Yb:fiber  oscillator  and  amplifier. 


wavelength  Xp  substantially  longer  than  the  input  signal  wavelength  Ai  must  be  used. 
For  input  signals  in  the  1.55-pm  band,  this  therefore  means  a  pump  wavelength 
Ap  >  1.8  H m,  which  is  in  a  spectral  region  where  commercial  laser  sources  with 
the  desired  properties  (tunable,  single-frequency,  output  power  up  to  1  W)  are  not 
available. 

We  therefore  developed  a  pump  source  for  long-wavelength-pumped  quantum  fre¬ 
quency  converters  based  on  a  monolithic  PPLN  optical  parametric  oscillator  (OPO) 
[97].  In  our  OPO,  the  path  of  the  resonant  signal  wave  is  entirely  within  the  PPLN 
QPM  grating  [108].  This  approach  has  several  advantages  including  mechanical  and 
thermal  stability,  simplicity  of  operation  and  alignment,  and  the  possibility  of  low 
losses  since  the  monolithic  resonator  has  no  lossy  internal  interfaces,  and,  impor¬ 
tantly  for  mid-IR  operation,  has  no  (possibly  absorptive)  air  path.  This  monolithic 
configuration  was  accomplished  by  applying  an  angled  spherical  polish  and  a  highly 
reflective  (HR)  coating  to  two  of  the  facets  of  the  PPLN  crystal,  and  using  total 
internal  reflection  on  one  side  of  the  crystal  to  form  a  closed  resonant  path  [108],  see 
Fig.  (5.1).  To  ensure  that  the  resonator  axis  was  within  the  1-mm-thick  crystal  aper¬ 
ture,  the  angles  of  each  curved  facet  was  tested  optically  during  the  OPO  fabrication 
process  via  a  reflection  measurement  with  a  HeNe  laser. 

In  our  OPO,  the  mirror  curvature  R  =  40  mm,  and  the  crystal  has  dimensions 
52x5x1  mm.  The  PPLN  crystals  were  poled  by  Crystal  Technology,  Inc.,  with  a 
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Figure  5.2:  Monolithic  0P0  signal  wavelength  versus  temperature  (squares)  and  fit 
to  LiNbC>3  Sellmeier  equation  [109]. 

QPM  period  of  Ac  —  30.9  /irn,  and  were  fabricated  into  resonators  by  the  Ginzton 
Lab  Crystal  Shop.  We  resonate  the  short-wave  signal,  which  is  tuned  around  1.8  /irn. 
The  coating  for  the  curved  input  crystal  facet  is  highly  reflective  (HR)  between  1.7 
and  2.0  /irn  (reflectance  R  >  99.9%),  highly  transmissive  (HT)  between  2.4  and  2.9 
/irn  (R  <  5%),  and  has  R  ~  5%  at  1.064  /irn.  The  coating  for  the  output  facet  has 
R  ~  99.5%  between  1.7  and  1.9  /irn,  R  <  5%  between  2.4  and  2.9  /irn,  and  R  ~  5% 
at  1.064  /irn.  The  focusing  parameter  is  defined  as  £,■  =  L/(kjWj0)  (j  G  {i,  s,  p}), 
where  L  is  the  QPM  grating  length,  kj  =  ‘lixrij/Xj  is  the  wave- vector  at  wavelength 
A  j  ,  and  Wjo  is  the  1/e2  intensity  radius  of  wave  j  at  the  beam  waist.  We  focused  the 
pump  to  slightly  looser  than  confocal  (£p  «  0.9),  while  the  cold-cavity  signal  mode 
defined  by  the  geometry  of  the  PPLN  crystal  is  focused  slightly  tighter  than  confocal 

(&  «  1-3). 

A  complete  description  of  the  measurements  of  the  performance  of  the  monolithic 
OPO  has  been  given  in  [97].  The  OPO  had  a  threshold  pump  power  of  approximately 
1  W,  and  at  a  maximum  pump  power  of  7.0  W,  we  observed  a  signal  output  power 
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of  0.98  W  with  a  pump  depletion  of  80%.  Here,  we  provide  data  on  the  tuning  of  the 
OPO  signal  wavelength  As  versus  the  device  temperature  T.  Our  results  are  shown 
in  Fig.  5.2  (squares)  along  with  a  fit  to  the  dispersion  relation  of  LiNbOs  (solid) 
[109].  The  fit  includes  a  +3°  C  oven  temperature  offset  to  match  theory.  By  tuning 
the  device  temperature  from  135  to  200°  C,  the  signal  wavelength  tuned  from  1.76 
to  1.94  /mi.  We  used  a  Michelson- interferometer-based  wavelength  meter  and  optical 
spectrum  analyzer,  which  gave  an  instrument- limited  linewidth  of  0.04  nm  (3.66  GHz). 
This  linewidth  is  approximately  2.7  times  the  free-spectral  range  of  the  OPO  cavity 
(1.34  GHz).  However,  we  are  confident  that  the  OPO  runs  in  a  single  longitudinal 
mode  at  pump  powers  near  threshold.  Single-mode  operation  can  be  inferred  from 
the  high  signal  depletion  we  measured  in  the  upconversion  detector  described  in  the 
next  section;  if  multiple  longitudinal  modes  were  oscillating,  mode  beating  would 
cause  the  observed  depletion  to  reach  a  maximum  value  of  approximately  5  dB  [110], 
whereas  we  observed  a  signal  depletion  of  41  dB. 

We  note  the  device  design  was  suboptimal  in  that  the  cavity  was  somewhat  un¬ 
dercoupled  compared  to  the  losses.  Based  on  our  estimate  of  the  absorption  losses, 
the  available  pump  power  of  8  W,  and  90%  pump  depletion,  an  output  coupler  of 
around  4%  would  lead  to  optimal  signal  output  power  of  approximately  3.4  W.  Due 
to  the  undercoupling,  we  observed  several  additional  effects  which  caused  changes  to 
the  signal  power  and  its  spectrum.  These  effects  included  an  optical  parametric  am¬ 
plification  process  involving  a  backwards-propagating  THz-frequency  wave,  as  well 
as  stimulated  Raman  scattering  of  the  intracavity  signal.  These  effects  (discussed 
in  detail  in  [97],  in  addition  to  an  intrinsic  modulation  instability  of  singly  resonant 
OPOs  [111],  resulted  in  our  needing  to  pump  the  OPO  close  to  threshold  (<  1.5x) 
in  order  for  the  OPO  to  run  stably  in  a  single  line. 
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Figure  5.3:  Experimental  setup  for  demonstration  of  long-wavelength-pumped  up  con¬ 
version  single-photon  detection.  The  single-photon-level  signal  from  a  tunable  C-band 
ECDL  is  combined  with  the  long-wavelength  pump  in  an  off-chip  WDM  coupler.  The 
upconverted  signal  is  outcoupled  into  free  space,  is  filtered  using  a  dichroic  mirror 
(DM),  prism,  and  band-pass  filter  (BPF),  and  focused  onto  a  Si  APD. 


5.3  Long- wavelength-pumped  upconversion  detec¬ 

tor  at  1.55  //rr l 

We  developed  an  upconversion  single-photon  detector  for  the  1.5- /mi  telecommunica¬ 
tions  band  based  on  an  RPE  PPLN  waveguide.  In  this  section,  we  first  describe  the 
device  and  its  classical-signal  characterization,  followed  by  a  summary  of  the  upcon¬ 
version  detector  system  performance  and  observed  noise  properties.  The  system  and 
its  performance  have  been  published  in  [67]. 


5.3.1  Experimental  setup  and  characterization 

The  experimental  setup  for  our  long-wavelength-pumped  upconversion  detector  is 
shown  in  Fig.  5.3.  We  fabricated  RPE  PPLN  waveguides  designed  for  sum-frequency 
generation  of  an  1850-nm  pump  and  1550-nm  signal.  The  waveguides  had  total  length 
L  —  52  mm,  and  were  fabricated  with  a  poling  period  Aq  =  18.4  or  18.8  /im.  The 
waveguides  had  at  their  inputs  spatial-mode  Liters  designed  to  match  the  mode  size  of 
SMF-28  optical  fiber,  and  were  fiber-pigtailed  with  coupling  losses  of  approximately 
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Figure  5.4:  TrnDFA  noise  spectrum  between  1450  and  1580  nrn,  measured  with  an 
optical  spectrum  analyzer  with  a  resolution  bandwidth  of  2  nrn,  showing  the  TrnDFA 
pump  line  at  1567  nm  and  shelf  of  spontaneous  emission  down  to  approximately 
1470  nm. 

0.8  dB.  We  measured  propagation  losses  below  0.2  dB  cm-1  using  the  Fabry-Perot 
technique  (see  Sec.  3.3).  The  waveguide  facets  were  flat-polished  and  antireflection 
coated  to  eliminate  interference  effects  and  improve  the  sum- frequency  (SF)  through¬ 
put.  The  waveguides  were  held  in  a  temperature-controlled  oven. 

The  tunable  pump  source  we  used  for  this  experiment  was  based  on  the  monolithic 
PPLN  optical  parametric  oscillator  described  in  the  previous  section.  For  stability,  we 
operated  the  OPO  close  to  threshold,  and  used  the  narrow  linewidth  signal  wave  to 
seed  a  Tm-doped  fiber  amplifier  (TrnDFA,  IPG  Photonics)  which  could  produce  up  to 
800  mW  of  output  power.  Before  inserting  the  long-pass  filters  into  the  pump  path, 
a  strong  noise  signal  overwhelmed  the  upconversion  of  the  signal  light,  saturating 
the  Si  APD.  We  investigated  the  noise  spectrum  of  the  TrnDFA  by  separating  the 
long-  and  short- wavelength  components  using  a  C/L-band  micro-optic  WDM  and 
coupling  the  C-band  components  into  an  optical  spectrum  analyzer.  The  TrnDFA 
noise  spectrum  is  shown  in  Fig.  5.4  for  a  3-mW  seed  at  1800  nm  and  output  power  of 
700  mW.  The  TrnDFA  is  pumped  by  an  Er:hber  laser  at  1567  nm  [112],  but  a  shelf 
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of  spontaneous  emission  noise  is  observed  at  all  wavelengths  down  to  approximately 
1470  nm.  Another  popular  pump  scheme  for  Tin  lasers  is  to  pump  at  790  nm; 
a  subsequent  cross-relaxation  process  enables  the  generation  of  two  laser  photons 
for  each  pump  photon.  Investigation  of  cross-relaxation  Tm  lasers  also  shows  C- 
band  spontaneous  emission  [113].  For  low-noise  upconversion,  the  pump  light  and 
spontaneous  emission  must  be  filtered.  Previous  work  on  upconversion  using  Tm  laser 
systems  used  fiber-optic  WDMs  as  filters  [32],  Here,  we  use  a  series  of  three  free-space 
long-pass  filters  placed  into  a  fiber  bench  to  achieve  the  very  high  extinction  (>100 
dB)  needed  to  eliminate  noise  due  to  SFG  of  noise  photons  due  to  pump  processes 
with  the  strong  long-wavelength  pump.  With  our  devices  tuned  to  a  phasematched 
wavelength  of  1555  nm,  the  phasematching  acceptance  bandwidth  of  the  upconverter 
is  itself  a  very  effective  filter  of  the  1567-nm  TmDFA  pump;  we  observe  50  dB  rejection 
of  upconversion  of  the  1567-nm  noise  photons  versus  1555  nm  signal  photons,  which 
is  consistent  with  a  prediction  based  on  a  domain- disorder-induced  QPM  pedestal 

[75]. 

The  1550-nm-band  signal  was  produced  by  a  tunable  external-cavity  diode  laser 
and  a  series  of  calibrated  attenuators.  As  PPLN  RPE  waveguides  support  only  TM- 
polarized  modes,  polarization  controllers  were  used  to  rotate  the  polarization  of  both 
the  signal  and  pump  beams.  If  necessary,  polarization- independent  upconversion 
detectors  as  in  [114]  could  be  built  in  the  future.  The  long- wavelength  pump  and 
single-photon-level  signal  were  combined  off-chip  in  a  micro-optic  C/L-band  WDM 
(Oplink)  found  to  have  low  insertion  loss  for  both  the  signal  and  pump  wavelengths. 
For  measurements  of  the  detection  efficiency,  the  signal  level  was  set  to  106  photons 
s^1  at  the  entrance  of  the  WDM.  The  upconverted  SF  radiation  at  0J2  was  collected 
by  an  aspheric  lens  (/  =  8  mm).  The  pump  and  signal  were  separated  from  the  SF  by 
a  dichroic  mirror  and  the  pump  power  was  monitored  using  a  thermal  power  meter. 
The  upconverted  output  was  further  filtered  using  a  Brewster-angle  prism  and  50- 
nm-bandwidth  optical  bandpass  filter  to  reject  parasitic  second  harmonic  radiation  of 
the  strong  pump.  The  SF  light  was  focused  onto  a  Si  single-photon-counting  module 
(Perkin-Elmcr  SPCM-AQR-13).  The  transmission  of  the  optical  system  between  the 
waveguide  output  facet  and  the  SPCM  was  92%. 
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Figure  5.5:  Measured  signal-wavelength  tuning  (dots)  and  calculation  from  LiNbC>3 
Scllmeier  coefficients  (solid),  where  Xp  =  1800.0  nm  and  =  18.4  //rn. 

We  first  characterized  the  phasematching  of  the  upconversion  waveguides.  The 
pump  wavelength  Xp  was  fixed  at  1800  nm  and  the  signal  wavelength  was  swept. 
For  the  device  used  in  subsequent  experiments  with  =  18.8  fin i,  we  observed 
a  phasematching  FWHM  of  0.57  nm,  which  matches  the  expected  width  calculated 
from  the  LiNbC>3  dispersion  relations  [109]  and  the  47-nnn  QPM  grating  length. 
We  also  measured  the  temperature  tuning  of  the  phasematching,  which  is  shown  in 
Fig.  5.5.  Experimental  data  for  a  device  with  A^  =  18.4  //rn  are  represented  as  dots 
and  the  solid  curve  is  a  theoretical  calculation  based  on  a  temperature-dependent 
Scllmeier  relation  for  LiNbOs  [109].  The  phasematching  temperature  tuning  rate  for 
Xp  =  1800  nm  was  0.27  nm/°C. 

We  measured  the  conversion  efficiency  of  the  waveguide  as  a  function  of  both 
pump  power  Pp  and  signal  wavelength  Ai  for  a  constant  classical-level  input  signal 
power  of  20  /iW;  the  results  are  shown  in  Fig.  5.6(a).  We  reach  a  maximum  inter¬ 
nal  conversion  efficiency  rj  of  86%  which  is  consistent  with  our  observed  propagation 
losses  of  0.17  dB  cm-1  at  Ai  and  a  value  of  0.1  dB  cm-1  at  A2  and  our  observed  41  dB 
(99.99%)  signal  depletion.  Measurements  at  24  pump  powers  between  0  and  350  mW 
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Figure  5.6:  (a)  Experimentally  observed  and  (b)  numerically  predicted  conversion 
efficiency  versus  pump  power  and  signal  wavelength.  Simulation  parameters:  au  (a^) 
=  0.17  (0.1)  dB  cm-1,  and  Pmax  =  151  mW;  (c)  internal  conversion  efficiency  of 
PPLN  waveguide  (Ai  =  1554  nm,  A2  =  834  nm). 
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have  been  interpolated  to  provide  a  smooth  plot.  Figure  5.6(b)  shows  a  theoretical 
prediction  of  the  internal  conversion  efficiency  calculated  by  numerically  integrating 
the  coupled-wave  equations  for  SFG  (as  in  Chapter  2)  with  our  experimental  param¬ 
eters.  There  is  good  agreement  of  the  major  features  and  the  observed  theoretical 
and  experimental  conversion  efficiencies  match  very  well.  Figure  5.6(c)  shows  a  slice 
of  the  two-dimensional  data  at  the  QPM  peak.  The  data  (squares)  are  well  fit  by 
Eq.  (2.19),  with  r/nor  =  68%  W'1  cnr2  and  thus  Pmax  =  151  mW. 

5.3.2  Single-photon  performance 

We  next  characterized  the  performance  of  the  detection  system  by  attenuating  the 
signal  laser  to  a  level  of  106  photons  s^1  at  the  input  of  the  WDM  using  a  series 
of  calibrated  fiber-optic  attenuators.  After  blocking  the  pump  noise  with  long-pass 
filters,  we  could  measure  the  photon  detection  efficiency  (PDE)  and  noise  count  rate 
(NCR)  of  our  upconversion  detector.  Measurements  were  made  using  several  pump 
wavelengths  between  1796  and  1859  nm.  A  plot  of  the  PDE  (left  axes)  and  NCR 
(right  axes)  versus  pump  power  for  a  pump  wavelength  Xp  =  1810  nm  is  shown  in 
Fig.  5.7.  At  Pp  =  Pmax  =  151  mW  we  attain  a  maximum  PDE  of  37%.  The  measured 
conversion  efficiency  matches  our  expected  conversion  efficiency  by  combining  the 
losses  of  each  component,  as  shown  in  Table  5.2.  While  the  SPCM  datasheet  quotes 
an  expected  PDE  of  55%  at  834  nm,  we  find  that  a  value  of  58%  is  needed  to  explain 
our  observed  system  detection  efficiency.  With  technologically  feasible  improvements 
to  the  setup  and  waveguide,  and  a  detector  with  PDE  rjv  =  58%,  we  should  be  able 
to  obtain  an  internal  conversion  efficiency  of  94%  and  a  maximum  PDE  of  48%. 

As  seen  in  Fig.  5.7,  there  are  still  appreciable  noise  counts  (approximately  7  x  103 
counts/s)  for  Xp  =  1810  nm  and  Pp  =  Pmax-  This  result  contrasts  with  earlier  work  in¬ 
volving  an  1810-nm-pumped  upconversion  detector  [32],  In  [32],  there  was  insufficient 
pump  power  (only  60  mW  was  coupled  into  the  waveguide)  to  reach  Pmax,  and,  per¬ 
haps  more  importantly,  the  SFG  output  was  passed  through  a  monochromator  before 
being  routed  to  the  SPCM,  significantly  narrowing  the  spectrum  of  detected  pho¬ 
tons.  We  estimate  the  bandwidth  of  our  output  filtering  system  to  be  approximately 
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Figure  5.7:  Photon  detection  efficiency  (PDE,  circles)  and  noise  count  rate  (squares) 
of  upconversion  SPD  with  Ap  =  1.81  //rn.  Solid  blue  is  a  fit  to  PDE  = 

77max  sin2  ^7r / 2  sj Pp / Pmaxj ,  where  77max  =  37%  and  Pmax  =  151  rnW.  The  noise  count 
rate  at  peak  PDE  is  approximately  7  x  103  counts/s. 

20  nm.  This  was  narrow  enough  to  reject  any  pump  second-harmonic  generation,  but 
not  narrow  enough  to  filter  the  upconverted  signal  (with  linewidth  likely  well  below 
1  GHz)  from  any  upconverted  SRS  photons  (which  will  fill  the  74-GHz  (0.57-nm) 
upconversion  acceptance  bandwidth). 

In  Sec.  4.2,  we  discussed  noise  in  QFC  devices  due  to  spontaneous  Raman  scat¬ 
tering  (SRS).  We  varied  the  pump  wavelength  Ap  between  1796  and  1859  nm,  while 
keeping  the  device  temperature  constant  at  T  =  30  °C.  As  a  result,  the  signal  wave¬ 
length  tuned  from  Ai  =  1557  to  1511  nm.  For  each  pump-signal  frequency  difference 
A /  =  (cni  —  u}p)/2n,  we  measured  the  PDE  and  NCR  of  the  SPD  as  a  function  of  Pp. 
Fig.  5.8  shows  a  plot  of  the  NCR  at  Pp  =  Pmax  for  each  A/.  We  observe  a  significant 
decrease  in  the  NCR  as  the  frequency  difference  A /  grows.  For  the  largest  frequency 
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Tabic  5.2:  Loss  and  transmission  of  upconversion  detector  components.  The  column 
with  feasible  transmission  values  is  based  on  idealized  optical  components  and  im¬ 
proved  waveguides  with  pigtailing  losses  of  0.5  dB  and  propagation  losses  of  0.1  dB 
cm-1 .  _ 


Component 

Loss  (dB) 

Transmission 

Feasible  T 

Wavelength  combiner 

0.45 

0.97 

0.99 

Fiber  pigtailing 

0.8 

0.83 

0.9 

Conversion  efficiency 

0.66 

0.86 

0.94 

Optical  system 

0.36 

0.92 

0.99 

SPCM PDE 

2.4 

0.58 

W 

Total: 

0.37 

0.83  rju 

separation  A /  =  1241  cm-1,  our  noise  count  rate  was  approximately  1000  counts/s. 
In  Fig.  5.8,  the  intrinsic  DCR  of  the  SPD  used  in  this  experiment,  500  counts/s,  was 
subtracted,  and  as  such,  the  experimental  data  represent  only  the  noise  counts  due 
to  spontaneous  scattering  or  other  technical  noise  sources  such  as  stray  light. 

The  solid  curve  in  Fig.  5.8  is  a  calculation  of  the  expected  NCR  due  to  spontaneous 
Raman  scattering  rate  for  this  system.  The  agreement  of  experimental  and  theoretical 
results  is  quite  strong.  The  calculation  uses  the  Raman  susceptibility  (imaginary 
component)  XR,i  given  in  [91],  uses  a  Raman  effective  area  Aeff,R  =  20  /./m2  (calculated 
using  the  modal  fields  from  RPE  waveguide  simulations  and  assumed  to  be  constant 
over  the  spectral  range  shown  in  Fig.  5.8,  valid  to  within  5%),  and  uses  a  pump 
power  Pp  =  Pmax  =  151  mW.  The  only  unknown  in  the  calculation  was  the  effective 
conversion  efficiency  of  the  Raman  noise  photons.  The  value  used  in  the  calculation 
in  Fig.  5.8  was  72%.  This  value  is  higher  than  the  expected  conversion  efficiency, 
which  was  calculated  to  be  approximately  40-50%  by  a  semiclassical  simulation  of 
upconversion  of  noise  photons  generated  in  a  distributed  fashion  across  the  device. 
This  is  perhaps  unsurprising,  as  it  is  known  that  spontaneous  Raman  scattering  in  the 
input  fiber  is  also  a  contributing  source  of  noise  photons  for  upconversion  detectors 
[30]. 

The  noise  performance  of  the  1.5-/im  SPD  could  be  improved  using  two  major 
techniques.  The  first  technique  is  the  reduction  of  the  device  temperature.  As  dis¬ 
cussed  in  Sec.  4.2,  anti-Stokes  Raman  scattering  requires  thermally  excited  phonons 
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Pump  wavelength  Lp  (pm)  for  A,1  =  1 .55  pm 
1.8  1.85  1.9  1.95  2 


Frequency  shift  (co1  -  cop)/2n  (cm'1) 


Figure  5.8:  Noise  count  rate  of  1.5-pm  upconversion  SPD  versus  pump-signal  fre¬ 
quency  difference  (up  —  u}p)/2ir  and,  equivalently,  pump  wavelength  Xp  needed  for  a 
1.55-pm  signal  (dots).  The  blue  curve  is  a  calculation  of  the  expected  noise-count  rate 
based  on  known  LiNbOa  material  properties  (see  Sec.  4.2)  and  system  parameters. 

to  be  present  in  the  substrate.  The  Boltzmann  factor  exp(— HAuj/kT)  describing 
this  thermal  population  decreases  as  T  is  reduced.  For  this  experiment,  the  device 
operated  at  a  temperature  of  30°  C.  With  a  reduction  of  the  device  temperature  to 
—50°  C,  the  noise  count  rate  would  be  reduced  by  a  factor  of  8,  to  a  level  roughly 
equivalent  to  specified  DCRs  on  high  quality  thick-junction  Si  APDs. 

The  second  technique  by  which  the  noise  could  be  further  reduced  is  spectral  fil¬ 
tering.  We  note  that  in  our  experiment  no  special  effort  was  made  to  spectrally  filter 
the  upconverted  signal  beyond  rejection  of  parasitic  pump  SHG.  In  an  earlier  ex¬ 
periment  demonstrating  low-noise  upconversion  detection  pumped  at  1810  nm,  noise 
performance  near  the  DCR  of  the  SPD  was  achieved  by  filtering  the  upconverted 
light  using  a  monochromator,  at  significant  expense  to  the  system  throughput  [32], 
In  most  applications,  the  signal  of  interest  will  have  substantially  smaller  bandwidth 
than  the  acceptance  bandwidth  of  the  upconverter  (74  GHz  (0.57  nm)  for  5-cm-long 
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RPE  PPLN  waveguides  at  1.55  /mi).  Because  upconverted  SRS  noise  will  fill  the 
acceptance  bandwidth  while  the  signal  will  presumably  be  narrowband,  the  noise  can 
be  significantly  reduced  by  spectral  filtering.  High-diffraction-efficiency  holographic 
gratings  with  minimal  losses  have  been  used  as  post- up  conversion  filters  in  a  recent 
publication  [42],  In  Chapter  7,  we  will  discuss  the  use  of  a  volume  Bragg  grating 
filter,  which  was  able  to  reduce  the  noise  in  a  1.5-/tm- pumped  upconversion  SPD  for 
signals  at  1.3  /an  to  a  level  of  approximately  500  counts/s,  which  is  to  our  knowledge 
the  lowest  value  reported  for  these  wavelengths. 

5.4  Upconversion  detection  of  l-//m  radiation 

It  is  expected  that  for  very  large  pump-signal  frequency  differences  A /  the  noise- 
count  rate  should  decrease  due  to  the  lower  thermal  occupation  ratio  for  Raman 
phonons.  We  demonstrated  a  low-noise  upconversion  single-photon  detector  for  the 
l-/nn  spectral  region  through  the  use  of  a  1.5-/nn-band  pump. 

Single-photon  detectors  for  the  l-/nn  spectral  band  are  in  an  awkward  position 
between  the  regimes  of  applicability  of  direct  detection  via  either  Si  or  InGaAs  APDs. 
However,  they  are  important  for  several  applications.  Although  they  would  have 
natural  applicability  in  remote  sensing  (due  to  the  prevalence  of  l-/tm  lasers  for  these 
systems),  we  constructed  a  l-/tm  upconversion  detector  for  use  in  an  experiment 
designed  to  demonstrate  the  feasibility  of  P  neutral-donor  bound  excitons  in  Si  as  a 
platform  for  quantum  computation  [115]. 

Our  l-/un  upconversion  SPD  was  based  on  a  6.7-cm-long  RPE  PPLN  waveguide 
with  a  poling  period  Aq  =  10.22  /on.  The  primary  spectral  range  of  interest  was 
near  1080  nm,  but  by  a  combination  of  temperature  tuning  (between  30  and  120°  C) 
and  pump  wavelength  tuning  (between  1.53  and  1.57  /an)  we  demonstrated  spectral 
coverage  between  1053  and  1087  nm.  The  device  achieved  a  maximum  depletion  of 
approximately  20  dB,  most  likely  limited  by  higher-order  spatial  mode  content  at  Ai 
as  no  directional  coupler  was  used  at  the  device  input.  Propagation  losses  at  1.55  /an 
were  measured  to  be  0.12  dB/cm.  The  maximum  internal  conversion  efficiency  was 
81%,  at  a  pump  power  Pp  =  Pmax  ~  50  mW.  Pmax  was  substantially  lower  than 
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Figure  5.9:  Measured  system  photon  detection  efficiency  (circles)  and  noise  count 
rate  (squares)  for  upconversion  SPD  with  Xp  =  1.54  /mi  and  Ai  =  1078  nm,  showing 
very  low  noise  counts  due  to  a  large  pump-signal  frequency  difference. 


for  the  1.5-yum  upconversion  detector  discussed  above  due  to  both  the  longer  device 
length  and  the  scaling  of  the  normalized  efficiency  discussed  in  Chapter  2  with  optical 
frequencies. 

We  constructed  an  upconversion  SPD  using  a  similar  setup  to  that  shown  in 
Fig.  5.3,  where  in  this  case  the  C-band  pump  radiation  was  produced  by  a  tunable 
telecom-band  ECDL  and  an  Erbium-doped  fiber  amplifier  (EDFA).  The  EDFA  had 
spectral  components  of  its  emission  in  the  range  of  1.0  to  1.1  /im  which  initially  dom¬ 
inated  the  noise  performance  of  the  system.  A  tunable  band-pass  filter  (bandwidth 
0.25  nm)  and  two  fiber-optic  WDMs  were  used  to  remove  the  1- /mi-band  spectral 
components  of  the  EDFA  emission  before  combining  the  pump  with  the  signal  ra¬ 
diation  in  an  off-chip  WDM.  The  target  radiation  at  A2  ~  630  nm  was  filtered  and 
routed  to  a  Si  APD  for  detection. 
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System  performance  measurements  were  made  with  two  commercial  single-photon 
detectors.  Preliminary  measurements  were  done  with  a  Perkin-Elmer  SPCM-AQR14 
detector  (PDE  at  A2  =  630  nm:  70%,  D  =  200  counts/s),  and  we  observed  a  sys¬ 
tem  detection  efficiency  of  23%,  with  a  total  NCR  (including  detector  DCR)  of  500 
counts/s.  For  an  ultra-low-NCR,  we  also  used  a  thin-junction  Si  APD  (id  Quan- 
tique  idlOO:  PDE  at  A2  ~  20%,  D  =  20  counts/s).  Our  measured  PDE  and  NCR 
are  shown  in  Fig.  5.9,  along  with  a  fit  to  the  usual  sin2(-)  relation  for  the  PDE.  At 
Pp  =  Pul ax  =  50  mW,  we  observed  a  maximum  system  detection  efficiency  of  5.2% 
and  a  NCR  of  50  counts/s.  The  main  contributions  to  the  low  system  PDE  were 
the  input  coupling  loss  of  1.1  dB  ( T  =  0.77),  output  collection  and  filtering  losses  of 
1.4  dB  (T  =  0.73),  the  quoted  photon  detection  efficiency  of  the  APD  of  20%,  and 
Fresnel  losses  of  0.7  dB  (T  =  0.86)  due  to  the  uncoated  waveguide  facet.  System 
improvements,  and  the  use  of  the  MPD  (PDM  series)  Si  APD  with  40%  detection 
efficiency  at  630  nm  would  result  in  a  system  PDE  of  27%.  The  results  of  this  sec¬ 
tion  demonstrate  that  with  judicious  choice  of  pump  wavelength  Xp  with  respect  to 
the  input  signal  Ai,  upconversion  detection  (and  QFC  in  general)  can  be  done  with 
extremely  low  added  noise,  even  with  no  special  attempts  at  output  filtering.  In  the 
next  chapter,  we  describe  a  long-wavelength-pumped  downconverter  for  the  interface 
of  a  solid-state  single-photon  emitter  at  910  nm  with  the  1 . 5- //m  telecom  band. 


Chapter  6 

Downconversion  of  single  photons 
from  quantum  dots 


Sources  of  single  photons  have  numerous  applications  in  quantum  optics  and  quan¬ 
tum  information.  The  best-developed  single-photon  (SP)  sources  emit  at  wavelengths 
in  the  visible  or  near-visible  spectral  region,  while  many  applications,  particularly  in 
quantum  communications,  would  benefit  from  sources  in  the  1.3-  or  1.55  /irn  telecom¬ 
munications  bands.  In  this  chapter  we  review  the  physics  of  quantum  dots  and  de¬ 
scribe  their  use  as  single-photon  sources  and  spin  qubits.  We  then  describe  an  experi¬ 
ment  in  which  single-photons  emitted  from  an  InAs  quantum  dot  at  A2  ~  910  nm  are 
downconverted  to  Ai  =  1.56  /irn  via  interaction  with  a  strong  pump  at  Xp  =  2.16  /irn 
in  a  PPLN  waveguide.  For  the  first  time,  to  our  knowledge,  we  are  able  to  downcon- 
vert  single  photons  from  a  quantum  emitter  to  the  1 . 5- /irn  telecommunications  band 
while  preserving  the  single  photon  character  and  demonstrating  coherent  control  of 
the  stationary  quantum  two-level  system. 

6.1  Quantum  dot  single-photon  sources 

The  earliest  single-photon  emitters  at  optical  frequencies  were  based  on  heralded 
single  photons:  first  from  atomic  cascades  [116]  and  later  by  SPDC  [117].  A  true 
triggered  (rather  than  heralded)  SP  source  would  be  due  to  the  excitation  of  a  single 
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atom  or  atom-like  system  with  a  single  valence  electron.  The  first  demonstration  of  a 
triggered  SP  source  including  a  measurement  of  the  photon  correlations  G^(t)  was 
in  1999  [118]  and  was  based  on  the  fluorescence  of  an  isolated  single  molecule  at  a 
temperature  of  1.8  K. 

Quantum  dots  are  islands  of  narrow-bandgap  semiconductor  embedded  in  a  larger- 
bandgap  semiconductor  matrix,  an  energy  schematic  of  which  is  shown  in  Fig.  6.1.  In 
the  InAs/GaAs  material  system,  if  a  layer  of  InAs  is  grown  by  molecular  beam  epitaxy 
(MBE,  [119])  on  a  GaAs  substrate,  due  to  the  lattice  mismatch  between  InAs  and 
GaAs,  after  the  growth  of  a  conformal  wetting  layer  of  approximately  2-nm  thickness, 
the  InAs  forms  islands  with  diameters  between  20  and  40  nm  and  heights  of  between 
4  and  7  nm,  depending  on  the  growth  conditions  [120].  If  the  growth  flux  is  low  and 
the  growth  is  stopped  soon  after  the  wetting  layer  forms,  the  quantum  dots  (QDs) 
can  be  kept  to  a  low  area  density  of  less  than  1  QD//im2;  the  low  QD  density  allows 
the  possibility  of  isolating  single  dots  in  a  confocal  microscope  setup.  An  attractive 
potential  results  in  the  confinement  of  electrons  and  holes  in  these  quantum  dots  and 
results  in  quantized  energy  levels  (as  in  the  canonical  “particle-in-a-box”  potential): 
the  bound  states  of  electrons  and  holes  can  then  be  calculated  as  the  eigensolutions 
to  Schrodinger’s  equation.  An  attractive  feature  of  InAs/GaAs  quantum  dots  is  that 
there  appear  to  be  realistic  possibilities  for  the  construction  of  quantum  repeaters 
based  on  qubits  consisting  of  either  a  single  electron  spin  [5]  or  single  hole  spin  [121] 
manipulated  by  ultrafast  optical  techniques. 

The  full  details  of  the  physics  of  the  level  structure  of  quantum  dots  and  its 
optical  transitions  are  beyond  the  scope  of  this  dissertation,  but  have  been  studied  in 
detail  and  can  be  reviewed  in  Santori’s  recent  monograph  [120].  For  our  purposes,  it 
will  suffice  to  describe  a  quantum  dot  as  an  idealized  quantum  two-level  system,  as 
shown  in  Fig.  6.1.  The  ground  state  |0)  consists  of  a  single  electron  in  the  conduction 
band.  There  exists  an  optical  transition  into  the  excited  state,  which  consists  of  an 
additional  exciton  (electron-hole  pair).  For  the  quantum  dot  we  studied,  the  energy 
difference  between  |0)  and  |1)  was  E1  —  E0  =  hw2  ~  1.3  eV,  corresponding  to  a 
transition  wavelength  A2  ~  910  nm. 
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Figure  6.1:  (a)  Band  structure  of  a  quantum  dot  in  its  ground  and  excited  states;  (b) 
Two-level  picture  of  same,  showing  optical  transition  at  A2  ~  910  nm. 


There  are  several  approaches  to  exciting  quantum  dots  for  SP  emission.  The  ini¬ 
tial  approach  taken  historically  was  to  use  incoherent  excitation  by  a  femtosecond 
Ti:  Sapphire  laser  above  the  GaAs  bandgap  at  liquid  helium  temperatures  (approxi¬ 
mately  780  nm)  [122,  123].  This  method  of  excitation  has  the  advantage  that  carrier 
generation  is  very  efficient  but  can  tend  to  result  in  unwanted  excitations  of  other 
nearby  dots  or  of  higher-energy  excited  states  of  the  QD.  A  potentially  cleaner  ap¬ 
proach  is  to  use  an  excitation  pulse  on  resonance  with  the  desired  emission  line  [124] . 
To  facilitate  spatially  filtering  the  quantum  dot  emission  from  the  strong  laser  exci¬ 
tation,  the  excitation  pulse  was  incident  from  a  steep  angle.  High  exctinction  can 
also  be  achieved  using  crossed  polarizers.  A  third  approach  has  been  termed  quasi¬ 
resonant  excitation,  in  which,  as  in  resonant  excitation,  one  can  select  one  quantum 
dot  out  of  a  sea  of  dots,  but  in  which  the  excitation  is  done  via  a  higher-level  excited 
state  (not  shown  in  Fig.  6.1,  and  not  yet  well  characterized)  so  that  the  excitation 
pulse  can  be  spectrally  filtered  from  the  quantum  dot  emission  [125]. 

For  initial  characterization  of  the  quantum  dot,  we  used  quasi- resonant  excitation 
at  895  nm.  To  determine  whether  the  emission  is  from  a  single  emitter  or  an  ensemble, 
the  standard  technique  is  to  measure  an  intensity  autocorrelation  g^(r)  in  a  Hanbury 
Brown-Twiss  intensity  interferometer  [126],  shown  schematically  in  Fig.  6.2(a).  The 
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Figure  6.2:  (a)  Measurement  of  photon  statistics  using  a  beamsplitter  and  two  detec¬ 
tors;  (b)  intensity  autocorrelation  g^2\r)  of  emission  from  the  quantum  dot  used  in 
the  downconversion  experiment. 


QD  is  excited  with  a  repetition  period  of  tr  =  13  ns.  If  it  emits  a  stream  of  single 
photons,  one  should  never  observe  coincidences  between  two  detectors  placed  on  the 
output  ports  of  a  50/50  beamsplitter.  The  measured  coincidence  data  for  the  quantum 
dot  used  in  this  experiment  is  shown  in  Fig.  6.2(b),  in  which  we  see  a  very  small  peak 
at  r  =  0  when  compared  with  the  surrounding  “accidental”  peaks  corresponding  to 
photons  emitted  from  different  excitation  events.  From  computing  the  central  peak 
area  and  comparing  it  to  the  average  areas  of  the  accidental  peaks,  we  compute  a 
value  c/2)(0)  =  0.13. 
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6.2  Classical  characterization  of  the  downconver- 
sion  system 

In  contrast  to  the  work  described  in  the  previous  chapter  on  upconversion,  in  this 
experiment  we  make  use  of  an  additional  feature  of  nonlinear  optical  interactions: 
they  can  be  used  to  rapidly  gate  the  detection  of  a  certain  event  with  higher  time 
resolution  than  would  be  achievable  by  direct  detection.  This  concept  is  effectively 
the  same  as  that  used  to  measure  femtosecond  optical  phenomena  (by  optical  gating) 
when  detectors  with  femtosecond  response  times  are  not  currently  technologically 
feasible  [127].  In  this  system,  when  the  quantum  dot  is  placed  in  a  magnetic  field 
B  =  6  T,  the  spin  of  the  electron  precesses  with  a  period  T  =  30  ps  [5].  If  one  hopes 
to  temporally  resolve  events  on  the  time  scale  of  this  precession,  it  is  not  possible 
to  use  conventional  single-photon  detectors,  as  these  have  timing  jitters  of  typically 
At  =  50  ps  or  more.  However,  the  gating  of  the  single  photon  emitted  by  the  quantum 
dot  with  a  short  optical  pulse  can  be  done  with  time  resolutions  of  a  few  picoseconds 
or  better,  limited  only  by  group-velocity  dispersion  effects. 

A  system-level  schematic  diagram  of  the  experiment  is  shown  in  Fig.  6.3.  The 
quantum  dot  is  excited  on  resonance  either  with  a  2-ps  pulse  from  the  Ti:Sapphire 
laser  shown  in  the  diagram  or  with  a  100-ps  pulse  from  a  second  Ti:Sapphire  laser  in 
which  the  repetition  rate  //?  has  been  locked  to  the  first  through  active  control.  The 
clock  signal  /#  =  76  MHz  from  the  2-ps  mode-locked  Ti: Sapphire  laser  serves  as  the 
master  clock  for  the  experiment. 

To  drive  the  pulsed  frequency  downconversion  process,  a  short-pulse  source  at 
Xp  ~  2.2-//1U  source  was  needed.  We  derived  this  source  from  the  master  TkSapphire 
2-ps  oscillator  through  the  use  of  a  DFG  process  in  a  5-cm-long  bulk  PPLN  (or 
MgO:PPLN)  crystal.  Approximately  300  mW  of  average  power  at  910  nm  was  mixed 
with  a  C-band  signal  to  generate  ps-long  pulses  at  2.2  /irn.  In  initial  experiments, 
the  C-band  source  consisted  of  a  cw  signal  from  an  ECDL  which  was  amplified  using 
an  EDFA.  At  its  maximum  drive  current,  the  EDFA  produced  a  cw  output  power 
of  approximately  2  W,  and  we  generated  a  2.2-/im  signal  with  a  peak  power  of  ap¬ 
proximately  16  W.  To  obtain  higher  peak  power  at  Xp,  we  also  utilized  a  strategy 
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Figure  6.3:  Experimental  setup  for  downconversion  of  single  photons  from  quantum 
dots.  Abbreviations:  HWP,  half-wave  plate;  PBS,  polarizing  beam  splitter;  SNSPD, 
superconducting  nanowire  single-photon  detector. 


where  the  ED  FA  seed  signal  was  pulsed  at  fn  with  a  short  optical  pulse  width.  This 
modulation  was  achieved  by  using  the  TLSapphire  electrical  clock  signal  to  synchro¬ 
nize  a  pulse-pattern  generator  (PPG).  The  10-GHz  PPG  produced  300-ps  electrical 
pulses  which  were  used  to  drive  an  electro-optic  intensity  modulator.  The  modulated 
optical  signal  was  then  first  amplified  in  a  low-power  EDFA  (preamplifier)  designed 
with  a  low  input  saturation  power  of  -23  dBm.  The  output  of  the  preamplifier  had 
an  average  power  of  approximately  10  mW  and  was  used  to  seed  the  the  same  power 
amplifier  as  had  been  used  with  cw  C-band  light.  The  pulsed  seed  enabled  substan¬ 
tially  higher  peak  powers  in  the  C-band,  and  therefore  produced  much  higher  peak 
powers  at  2.2  pm,  as  also  shown  in  Fig.  6.4. 

The  crystals  used  in  the  bulk  DFG  source  had  a  length  of  L  =  5  cm.  This  length  is 
substantially  longer  than  the  group-velocity  walkoff  length  between  the  input  signals 
Tgvm  =  T0/5V12  =  1.07  cm,  where  To  =  2  ps  is  the  TLSapphire  pulsewidth  and 
5vij  =  uj1  —  wy7,  where  tp  is  the  group  velocity  at  A*  with  i  G  {p,  1,  2}.  As  such, 
the  pulsewidth  of  the  generated  pulse  at  2.2  pm  is  determined  by  the  group  velocity 
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Figure  6.4:  Bulk  DFG  results  for  generation  of  pump  radiation  at  \p  ~  2.2  /mi.  For 
the  “Pulsed  EDFA  Seed”  case,  the  radiation  from  the  ECDL  was  modulated  with  a 
pulse  width  of  300  ps  to  increase  the  peak  power  extracted  from  the  EDFA. 

walkoff,  and  has  been  numerically  simulated  to  be  8.6  ps  full-width  at  l/e2-intensity 
maximum  using  a  split-step  Fourier  simulation  technique  [128].  The  fact  that  the  time 
resolution  of  the  experiment  can  be  determined  by  the  pulsewidth  of  the  2.2-/im  gating 
pulse  provides  a  convenient  control  variable.  We  have  procured  MgO:PPLN  crystals 
with  QPM  gratings  of  three  different  lengths  (1  cm,  3  cm,  and  5  cm)  positioned 
adjacently  on  the  same  5-cm  substrates  to  allow  control  over  the  time  resolution. 

Again  referencing  Fig.  6.3,  the  generated  2.2-/mn  short  pulses  are  sent  through 
a  variable  optical  delay  line  and  then  coupled  into  optical  fiber.  In  the  quantum 
dot  path,  the  single-photon  emission  is  collected  with  a  polarization  perpendicular  to 
the  input  excitation  to  enable  polarization  filtering  of  the  strong  resonant  excitation. 
After  spectral  filtering  using  a  set  of  diffraction  gratings  and  filters,  the  single-photon 
emission  at  A2  ~  910  nm  is  fiber  coupled  using  HI-780  fiber,  with  approximately  60% 
coupling  efficiency. 

Owing  to  the  complexity  of  combining  two  very  dissimilar  wavelengths  in  the 
fundamental  spatial  mode  of  a  single  waveguide,  the  waveguide  design  used  in  this 
experiment  incorporated  a  directional  coupler.  Light  from  the  quantum  dot  at  A2  and 
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from  the  pulsed  source  at  Xp  =  2.2  /irn  enter  the  device  from  separate  fibers  into  two 
separate  waveguides.  The  two  fibers  are  terminated  in  a  silicon  V-groove  array  and  are 
pigtailed  to  the  chip  input.  The  QD  light  at  A2  goes  through  an  adiabatic  taper  and 
an  S-bend  to  bring  it  into  proximity  with  the  waveguide  carrying  the  pump  light  at  Xp 
where  they  enter  the  two  input  ports  of  a  directional  coupler.  The  directional  coupler 
was  designed  such  that  the  pump  light  was  coupled  over  to  the  adjacent  waveguide 
while  the  coupling  for  the  quantum  dot  light  was  low  (parameters:  waveguide  channel 
width:  5.4  /irn,  waveguide  center-to-center  spacing:  10  /irn,  directional  coupler  length: 
0.95  mm).  The  combined  pump  and  signal  light  then  entered  the  QPM  mixing  region, 
which  was  poled  with  =  21.9  /irn  and  had  length  Lq  =  4  cm.  At  the  output  of 
the  mixing  region,  the  converted  light  at  Ai  =  1550  nm  was  sent  through  another 
adiabatic  taper  and  mode  filter  and  coupled  into  an  optical  fiber  (SMF-28)  pigtailed 
to  the  exit  facet  of  the  waveguide.  Input  coupling  and  propagation  losses  at  910 
nm  were  observed  to  be  1.5  dB,  and  the  output  coupling  loss  at  1550  nm  was  also 
1.5  dB.  The  observed  coupling  loss  and  waveguide  propagation  loss  at  2.2  /irn  was 
approximately  7  dB.  There  is  an  OH"  absorption  feature  due  to  a  vibration-libration 
resonance  at  a  wavelength  of  2250  nm  [129].  The  intentional  proton  doping  to  create 
the  waveguide  enhances  this  absorption.  It  is  suspected  that  the  short-wavelength 
tail  of  this  absorption  feature  is  responsible  for  the  relatively  high  losses  observed 
when  compared  with  the  NIR  quantum  dot  emission  which  is  very  far  detuned  from 
the  absorption  resonance. 

We  characterized  temporal  resolution  of  the  conversion  by  performing  a  DFG 
cross-correlation  in  the  PPLN  waveguide.  Classical  pulses  from  the  2-ps  Ti:Sapphire 
laser  with  peak  powers  of  a  few  mW  were  sent  into  the  waveguide  along  with  the 
2.2-/im  pump  pulses  and  we  measured  the  generated  DFG  signal  at  1550  nm  as  a 
function  of  the  2.2-/mr  pulse  delay.  Our  results  are  shown  in  Fig.  6.5,  along  with  a 
simulation  based  on  a  split-step  Fourier  method.  From  this  data,  we  observe  that  for 
2-ps  Ti:  Sapphire  pump  pulses,  the  approximate  conversion  time  window  is  6-10  ps. 
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Figure  6.5:  Measured  DFG  radiation  at  Ai  =  1550  nm  as  a  function  of  delay  time 
between  2.2-/im  pump  pulse  and  910-nrn  probe  pulse  (dots)  and  split-step  Fourier 
simulation  (solid). 

6.3  Single-photon  experiments 

Following  classical-signal-level  characterization  of  the  conversion  processes,  we  at¬ 
tempted  single-photon  downconversion  experiments.  Again  referring  to  Fig.  6.3,  we 
utilized  a  superconducting  nanowire  single-photon  detector  (SNSPD)  to  detect  the 
down-converted  single  photons.  The  detector  was  packaged  with  a  SMF-28  fiber  optic 
input  and  was  housed  in  a  3He  cryostat.  The  detector  was  operated  at  a  temperature 
of  approximately  2  K  with  a  voltage  bias  of  3.3  V  [130].  The  SNSPD  had  a  detection 
efficiency  of  approximately  14%  and  a  dark-count  rate  of  approximately  40  counts/s. 

We  briefly  comment  on  the  filtering  techniques  used  after  the  waveguide  to  sepa¬ 
rate  the  outgoing  downconverted  signal  from  other  unwanted  spectral  components  of 
radiation.  One  difficulty  of  using  SNSPDs  is  that  they  are  sensitive  to  effectively  all 
wavelengths  of  light  shorter  than  the  mid-infrared.  There  are  multiple  intrinsic  and 
technical  noise  sources  that  must  be  filtered.  At  the  output  of  the  waveguide,  we  used 
a  fiber-optic  circulator  and  a  fiber  Bragg  grating  (FBG)  with  a  reflection  bandwidth  of 
approximately  2  nm  to  attempt  to  eliminate  all  but  the  converted  quantum-dot  single 
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photons.  The  loss  of  this  setup  was  approximately  1.5  dB  (dominated  by  the  circula¬ 
tor)  and  the  extinction  ratio  was  approximately  60  dB,  limited  by  the  directivity  of 
the  circulator.  Some  amount  of  both  pump  light  at  Xp  =  2.2  /im  along  with  its  second 
harmonic  at  1100  nm  were  therefore  present  at  the  output  port  of  the  circulator.  To 
filter  this  radiation  further,  we  employed  two  strategies:  tight-radius  bending  of  the 
optical  fiber  to  remove  the  long-wavelength  pump  via  radiation  losses  and  a  1350-nm 
long-pass  filter  on  a  Si  substrate  inserted  into  a  fiber  U-bench  to  remove  both  the 
parasitic  second  harmonic  light  and  any  unconverted  910  light  from  the  quantum  dot 
setup.  The  loss  of  the  U-bench  and  long-pass  filter  was  approximately  1.5  dB.  We 
observed  a  saturation  of  the  conversion  efficiency  at  a  2.2-/im  peak  pump  power  of 
approximately  3  W  in  the  input  fiber.  With  these  filtering  mechanisms  in  place,  at 
this  power  level  we  observed  noise  count  rates  due  to  spontaneous  scattering  or  other 
leakage  processes  below  1  count/s.  We  note  that  if  the  narrowband  FBG  filter  was 
removed,  a  significant  noise  count  rate  (in  the  range  of  approximately  100  counts/s) 
was  observed,  which  we  suspect  are  due  to  broadband  spontaneous  anti-Stokes  Raman 
scattering. 

Throughout  this  experiment,  we  perform  measurements  using  the  time-correlated 
single-photon  counting  technique  (TCSPC).  In  TCSPC,  the  arrival  time  of  the  electri¬ 
cal  pulse  from  a  single-photon  detector  is  compared  to  a  regular  trigger  signal,  in  our 
case,  from  the  clock  output  of  the  Ti:  Sapphire  master  oscillator.  Over  many  events, 
a  histogram  of  arrival  times  is  formed,  resulting  in  a  time-ensemble  measurement  of 
a  temporal  waveform  containing  less  than  one  photon  per  clock  cycle  on  average. 

We  estimate  the  incident  photon  rate  entering  the  PPLN  waveguide  from  the 
quantum  dot  at  a  level  of  approximately  6  x  104  photons/s.  Using  a  TCSPC  mea¬ 
surement,  we  characterized  the  spontaneous  emission  decay  time  of  the  quantum  dot 
emission  and  found  it  to  be  well-described  by  a  single  exponential  with  an  approxi¬ 
mately  600  ps  time  constant.  We  then  turned  on  the  conversion  pulse,  which  sampled 
an  approximately  10-ps  time  window  of  the  quantum  dot  emission.  Experimental  re¬ 
sults  of  the  pulsed  conversion  of  the  quantum  dot  signal  as  a  function  of  the  time 
delay  between  the  quantum  dot  excitation  pulse  and  the  conversion  pump  pulse  are 
shown  in  Fig.  6.6.  At  early  times,  we  observe  a  relatively  rapid  rise  in  the  converted 
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Figure  6.6:  Measured  net  photon  counts  at  1550  nm  as  the  radiation  from  the  quan¬ 
tum  dot  is  downconverted  and  the  delay  between  the  quantum  dot  excitation  and 
the  conversion  pump  pulse  is  varied.  The  temporally  un correlated  dark  counts  of  the 
SNSPD  have  been  subtracted. 

counts  as  the  quantum  dot  is  excited  by  a  100-ps-long  excitation  pulse.  At  later 
times,  we  see  that  the  converted  light  exhibits  a  tail  of  emission,  corresponding  to  the 
600-ps  spontaneous  emission  lifetime  of  the  quantum  dot.  The  dark-count  level  of 
the  SNSPD  (counts  which  are  uncorrelated  in  time)  in  Fig.  6.6  has  been  subtracted, 
and  the  measurement  at  each  time  delay  was  integrated  for  five  minutes.  The  solid 
curves  represent  Gaussian  fits  to  the  data  at  each  time  delay,  and  all  have  temporal 
widths  in  the  range  of  100  ps,  the  SNSPD  timing  jitter.  The  dashed  solid  curve  is  an 
exponential  fit  to  the  peaks  of  the  count  distributions  for  each  time  delay,  and  has  a 
time  constant  of  600  ps.  The  signal  count  rate  at  the  emission  peak  (at  a  delay  of 
4.7  ps),  integrated  over  the  Gaussian  instrument  response  of  the  SNSPD,  is  approx¬ 
imately  20  counts/s.  This  count  rate  is  consistent  with  the  measured  rate  of  single 
photons  into  the  HI-780  input  fiber  of  approximately  6  x  104  photons/s,  a  temporal 
overlap  factor  of  approximately  1.7%,  and  the  loss  sources  already  delineated. 

The  measurement  described  above  of  tracing  out  the  temporal  shape  of  the  emitted 
waveform  from  the  sample,  while  exhibiting  the  exponential  spontaneous  emission  tail, 
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Figure  6.7:  Measured  photon  counts  at  910  nm  (blue  dots)  or  downconverted  to  1550 
nm  (green  squares)  as  a  function  of  excitation  power,  showing  Rabi  oscillations. 


does  not  demonstrate  the  coherent  manipulation  of  a  quantum  two  level  system.  To 
demonstrate  this,  we  did  a  Rabi  oscillation  experiment,  whereby  the  photon  count 
rate  was  measured  as  a  function  of  the  power  of  the  optical  excitation.  If  (referring 
to  Fig.  6.1)  the  quantum  dot  starts  in  state  1 0) ,  then  as  the  excitation  field  E(t)  = 
EA(t)elW2t  (where  we  have  defined  EA{t)  as  the  slowly- varying  field  envelope)  is  varied 
the  state  is  rotated  to  a  superposition  state: 

W)  =  cos  (0  1°)  +  sin  (Q  I1)  (6-1) 


where  the  rotation  angle 


6  =  k  EA{t')  d t' 


(6.2) 


where  the  coupling  constant  k  =  2 /i/ft,  where  /i  is  the  dipole  moment  of  the  system 
[131].  From  Eq.  (6.2),  we  note  that  the  rotation  angle  is  expected  to  be  proportional 
to  y/p.  If,  at  the  end  of  the  coherent  manipulation  via  the  resonant  pulse  E(t),  the 
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system  is  in  state  \rjl)  given  by  Eq.  (6.1),  we  will  measure  the  projection  |(,0|1)|2  = 


sin2  (0/2)  because  after  the  pulse  is  turned  off  the  upper  state  will  spontaneously 


decay  via  the  emission  of  a  single  photon,  which  can  then  be  measured. 

Our  experimental  results  are  shown  in  Fig.  6.7.  The  blue  dots  are  the  recorded 
photon  counts  of  910-nm  photons  emitted  from  the  quantum  dot  without  downcon- 
version.  As  we  can  see,  they  rise  from  a  low  value  when  the  excitation  power  is  very 
low  and  undergo  a  series  of  oscillations  as  the  state  has  been  coherently  manipulated 
through  rotations  of  run,  where  the  maxima  correspond  to  m  odd  and  the  minima 
correspond  to  m  even.  The  solid  blue  curve  in  Fig.  6.7  is  a  fit  to  an  empirical  relation 
of  the  form: 


(6.3) 


where  we  have  included  a  DC  offset  due  to  dark  or  leakage  counts,  a  background  level 
rising  as  a  power  law,  and  an  exponentially  damped  oscillatory  term.  We  then  turned 
on  the  conversion  pump  at  a  fixed  time  delay  and  measured  the  converted  counts  as  a 
function  of  the  quantum  dot  excitation  power.  The  results,  scaled  to  the  same  vertical 
scale,  as  the  non-converted  quantum  dot  light,  are  shown  as  green  squares.  We  see 
the  same  oscillations  in  the  count  rate  as  the  excitation  power  is  varied,  confirming 
that  the  downconverter  is  preserving  quantum  properties  of  the  source. 

An  important  aspect  to  prove  that  the  quantum  frequency  downconverter  main¬ 
tains  the  quantum  state  of  the  light  emitted  from  the  quantum  dot  is  the  investiga¬ 
tion  of  photon  statistics.  This  experimental  task  proved  to  be  challenging  as  the  low 
count  rates  (dominated  by  the  low  temporal  overlap  of  the  conversion  pulse  with  the 
quantum  dot  spontaneous  emission  waveform)  necessitated  very  long  integrations  of 
several  hours.  Again  referring  to  the  experimental  setup  in  Fig.  6.3,  we  now  seek  to 
measure  correlations  in  detected  events  between  the  Si  APD  which  is  sampling  the 
emitted  quantum  dot  before  the  waveguide  and  the  converted  photons  at  1550  nm 
measured  following  downconversion. 

We  used  the  temporal  resolution  of  the  detection  setup  to  eliminate  the  calculation 
of  coincidence  counts  from  the  Si  APD  and  the  SNSPD  at  times  that  do  not  corre¬ 
spond  to  the  proper  phase  of  the  clock  cycle  for  when  the  converted  photons  appear. 
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Figure  6.8:  Measured  coincidence  counts  as  a  function  of  time  delay  between  downcon- 
verted  photons  counted  on  the  SNSPD  and  non-converted  photons  measured  before 
the  waveguide  on  a  Si  APD.  We  observe  an  antibunching  dip  corresponding  to  a  value 
(/2)(0)  =  0.17,  showing  preservation  of  the  single-photon  character  of  the  quantum 
dot  emission  after  downconversion. 

Our  coincidence  data  is  shown  in  Fig.  6.8.  At  zero  time  delay,  we  observe  a  dip  in  the 
coincidences  corresponding  to  antibunched  statistics.  Computing  the  area  of  the  dip 
and  comparing  with  the  surrounding  accidental  peaks  gives  a  value  g^2\ 0)  =  0.17, 
which  is  only  slightly  higher  than  the  value  observed  with  quasi-resonant  excitation 
and  no  downconversion:  0.13.  We  should  note  that  neither  here  nor  in  the  non- 
converted  5,(-2')(t)  measurement  shown  in  Fig.  6.2(b)  has  any  background  subtraction 
been  done.  Any  detector  dark  counts  appearing  inside  the  time  window  of  the  SNSPD 
in  which  coincidences  may  be  registered  will  add  to  the  noise.  If  these  dark  counts 
were  subtracted  we  believe  that  the  measured  g^2\r)  would  probably  be  somewhat 
improved  when  compared  with  the  nonconverted  data  of  Fig.  6.2.  This  improvement 
would  be  due  to  the  fact  that  the  conversion  process  effectively  provides  filtering  in 
both  frequency  (due  to  the  limited  acceptance  bandwidth  of  the  downconverter  of 
approximately  0.2  nm)  and  the  strong  temporal  filtering  due  to  the  pulsed  conversion 
process,  where  no  temporal  filtering  was  done  in  the  experiment  without  conversion. 


CHAPTER  6.  SINGLE-PHOTON  DOWNCONVERSION 


93 


6.4  Summary 

To  summarize  the  work  described  in  this  chapter,  we  have  developed  a  waveguide 
mixer  and  short-pulsed  2.2-pm  pump  source  for  efficient,  temporally  resolved  down- 
conversion  of  single  photons  from  InAs/GaAs  quantum  dots  into  the  telecommu¬ 
nications  band.  We  have  confirmed  the  quantum  nature  of  the  emitter  through  a 
Rabi  oscillation  experiment,  whereby  the  oscillatory  behavior  observed  in  the  count 
rate  as  a  function  of  excitation  power  was  also  observed  upon  downconversion  of  the 
photons.  We  have  also  confirmed  the  single-photon  character  of  the  emission  by  mea¬ 
suring  antibunching  of  the  counts  in  a  cross-correlation  experiment.  The  use  of  a  very 
long-wavelength  pump  source  enabled  low-noise  frequency  conversion:  our  measured 
noise  count  rates  (counts  due  to  the  presence  of  the  2.2- /im  pump)  were  estimated  to 
be  approximately  1  count  per  second. 

To  the  best  of  our  knowledge,  this  experiment  represents  the  first  demonstration 
of  the  downconversion  of  a  stationary  quantum  emitter  to  the  1.5  /mi  telecommuni¬ 
cations  band.  There  have  been  two  recent  demonstrations  of  downconversion  to  the 
1.3-/tm  telecom  band.  The  first  experiment  was  a  component  of  a  quantum  memory 
based  on  trapped  8'Rb  atoms  where  the  downconversion  was  achieved  via  four-wave 
mixing  in  cold  Rb  vapor  [132] .  The  second  experiment  was  a  recent  demonstration  of 
the  downconversion  of  quantum  dot  single  photons  at  A2  =  710  nm  to  Ai  =  1310  nm 
using  a  strong  pump  at  1550  nm  [133].  However,  the  quantum  dots  used  in  this  ex¬ 
periment  were  based  on  InP/GalnP  and  there  have  been  no  demonstrations  of  spin 
control  in  this  system,  but  these  quantum  dots  are  suitable  for  electrically  triggered 
single-photon  generation  [134].  There  has  also  been  recent  work  on  a  1.5-/im  telecom- 
band  single-photon  stource  based  on  InAs/InP  quantum  dots  [135],  but  there  have 
been  no  attempts  at  spin  control  of  quantum  dots  with  C-band  transitions. 

For  the  experiments  described  in  this  chapter,  the  use  of  a  pulsed  pump  for  the 
downconversion  process  was  not  necessary  or  even  perhaps  advisable  if  demonstration 
of  low-noise  downconversion  was  the  only  goal.  The  ultimate  goal  of  this  experiment, 
however,  was  the  demonstration  of  entanglement  between  the  quantum  dot  electron 
spin  and  the  polarization  of  the  emitted  photon.  For  a  quantum  dot  in  a  magnetic 
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field,  the  Zeeman  effect  lifts  the  energy  degeneracy  of  the  two  spin  states  of  the 
electron,  which  we  denote  as  |  t)  and  |  J,),  where  the  energy  difference  corresponds  to 
a  frequency  5e  =  (E^  —  E±)/h  ~  2tt  x  30  GHz.  Therefore,  the  entangled  state  can  be 
written 

|T)  =  —  (|  t)  0  \iH-,u2  +  5e)  +  |  |)  <8)  \V ;a;2))  (6.4) 

From  Eq.  (6.4),  we  see  that  there  is  entanglement  not  only  between  the  photon  polar¬ 
ization  and  the  spin  state,  but  also  the  photon  frequency  and  spin  state.  Therefore, 
the  photon  frequency  provides  “which-path”  information  about  the  state  of  the  spin, 
and  therefore,  would  destroy  the  entanglement  unless  the  frequency  information  can 
be  erased.  In  a  pioneering  spin-photon  entanglement  experiment,  it  was  realized  that 
a  time-resolved  detection  (with  time  resolution  At  -C  l/5e  provides  such  a  a  quantum 
eraser  [4],  This  erasure  was  possible  in  the  experiment  of  Togan  and  coworkers  be¬ 
cause  in  the  NV-center  system  the  energy  splitting  (in  that  system,  provided  by  strain 
rather  than  Zeeman  splitting)  was  very  low  (approximately  300  MHz)  which  meant 
that  standard  single-photon  detectors  provided  sufficient  time  resolution  to  achieve 
the  quantum  erasure  effect.  In  our  experiment,  however,  where  5e  ~  2tt  x  30  GHz, 
a  time  resolution  of  approximately  30  ps  is  needed,  which  is  lower  than  is  available 
with  commercial  single-photon  detectors.  Hence,  we  employ  the  pulsed  frequency 
conversion  technique. 

We  have  recently  verified  the  entanglement  between  the  quantum  dot  spin  state 
and  the  downconverted  telecom-band  photon  [136].  Combining  spin-photon  entangle¬ 
ment  with  a  two-photon  interference  experiment  should  enable  entanglement  gener¬ 
ation  between  remote  quantum  dot  spins.  The  fact  that  the  spin-spin  entanglement 
will  be  mediated  by  telecom  band  single  photons  could  potentially  enable  entangled 
states  of  stationary  qubits  at  distances  of  tens  of  km,  compared  to  the  relatively  mod¬ 
est  distances  of  approximately  1  m  thus  far  achieved  in  the  entanglement  between 
two  trapped  ions  [137]. 


Chapter  7 

Advanced  quantum  frequency 
conversion  devices 


In  this  chapter,  we  move  past  devices  for  one-step  frequency  up-  and  downconver- 
sion  and  address  some  novel  device  designs  that  can  bring  additional  functionality  to 
single-photon  frequency  conversion  systems.  Broadly  speaking,  we  can  use  similar  de¬ 
vice  engineering  techniques  as  are  found  in  classical  optical  signal  processing  systems 
[9]  and  refocus  them  for  the  quantum  realm.  The  primary  applications  we  address 
are  to  enable  low-noise  frequency  conversion  between  widely  disparate  frequencies, 
and  to  enable  detection  and  processing  of  high-data-rate  modulated  signals. 

7.1  Cascaded  quantum  frequency  conversion  and 
applications 

In  this  section,  we  will  describe  aspects  of  the  theory  and  design  constraints  of  incor¬ 
porating  two  QPM  gratings  within  a  single  waveguide  for  cascaded  conversion.  With 
such  a  device  we  achieved  quasi-phasematching  of  two  SFG  processes,  enabling  87% 
conversion  efficiency  of  a  signal  from  1550  nm  to  570  nm  [65].  Sec.  7.1.2  describes 
an  application  for  a  high-speed  SPD  with  a  total  system  timing  jitter  of  66.9  ps. 
Sec.  7.1.3  presents  potential  applications  for  downconversion  of  single-photon  sources 
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Figure  7.1:  Schematics  of  (a)  frequencies  involved  with  and  (b)  PPLN  waveguide 
design  for  two-stage  upconversion  for  low-timing-jitter  single-photon  detection;  (c) 
PDE  of  commercial  low-jitter  Si  SPDs  versus  wavelength,  manufactured  by  MPD 
and  id  Quantique. 


in  the  visible  spectral  range  to  the  1.55  /irn  telecom  band  [48]. 


7.1.1  Cascaded  conversion  theory  and  device  design 

First,  we  briefly  develop  the  theory  of  a  cascaded  frequency  conversion  device.  We 
consider  here  the  case  of  cascaded  SFG,  which  has  been  demonstrated  and  will  be 
described  in  the  next  section.  We  consider  a  device  for  the  upconversion  of  a  signal 
at  cui  to  C02  via  a  two-stage  upconversion  process  using  a  single  optical  pump  at  cup, 
a  schematic  of  which  is  shown  in  Fig.  7.1(a).  However,  the  description  is  generic 
and  could  equally  well  be  applied  to  cascaded  DFG,  an  application  of  which  will  be 
described  in  Sec.  7.1.3.  Here,  the  input  signal  at  coi  is  Erst  summed  with  a  long- 
wavelength  pump  at  cop  to  produce  radiation  at  cuint  =  uj\  +  ojp  in  a  process  we  term 
SFG-A.  A  second  SFG  process  (SFG-B,  here,  within  the  same  PPLN  waveguide)  is 
then  used  to  sum  the  light  at  ow  with  the  same  pump  to  produce  the  target  radiation 
at  U2  =  c^int  +  ujp  =  ui\  +  2 ujp. 

We  consider  a  QPM  grating  with  Fourier  components  at  Ka  =  A =  kint  —  ki  —  kp 
and  Kb  =  A ks  =  k2  —  k\rA  —  kp.  The  amplitudes  of  each  frequency  component, 
which  we  refer  to  as  qa{z)  and  gs(z)j  may  vary  along  the  propagation  distance  z 
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of  the  device.  Ideally,  they  obey  the  constraint  |^(z)| 2  +  \gB(z)\ 2  =  1,  which  is  a 
manifestation  of  Parseval’s  identity  for  the  Fourier  transform  of  a  nonlinear  coefficient 
d(z)  whose  magnitude  (but  not  necessarily  whose  phase)  is  independent  of  z.  A  two- 
component  QPM  grating  could  be  fabricated  using,  e.g.,  a  phase-modulated  structure1 
[138]  or  an  optical  superlattice  approach  [139]. 

The  coupled  amplitude  equations  for  normalized  held  envelopes  in  a  two-stage 
frequency  conversion  process,  in  the  limit  of  an  undepleted  pump  and  negligible  prop¬ 
agation  loss,  are  given  by 

iVBgB{z)a-mtelAkBZ ,  (7.1a) 

i^BgB(z)a2e~lAkBZ  -  i-fAgA(z)a1elAkAZ,  (7.1b) 

i'yAgA(z)aiQte~lAkAZ ,  (7.1c) 

where  7 j  =  ^VnorjPp,  where  r]noi,j  is  the  normalized  conversion  efficiency  of  process 
j  G  {A,  B},  and  [|afc(z)|~]  =  photons/second  [67].  Analytical  solutions  of  Eqs.  (7.1) 
are  possible  for  either  an  abrupt  transition  where  gA^z)  =  1  for  0  <  z  <  LA  or  where 
the  Fourier  amplitudes  are  constant  with  z:  gA{z)  =  go  and  gB(z )  =  (1  —  gl)1^2 ■  We 
focus  here  on  the  case  of  sequential  gratings,  the  simplest  structure  to  design  and 
fabricate.  Assuming  incident  radiation  only  at  07,  and  perfect  QPM  for  both  SFG 
processes,  the  conversion  efficiency  rj2  (r/int)  to  frequency  co2  (tAnt)  are  given  as 

m  =  =  sin2b aLa]  sin2 [ye (T  -  LA)\  (7.2a) 

Vint  =  =  sin2  [7^]  cos2 [7b (A  -  LA)\,  (7.2b) 

from  which  we  see  that  complete  conversion  from  uj\  to  oo2  is  possible  only  when 
7 aLa  =  7b(T  —  La )  =  7t/2.  This  constraint  sets  the  optimal  grating  transition 
location  LA/L  =  ^B/{lA  +  7b)- 

Although  the  fabricated  devices  described  in  the  following  section  have  used  the 
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1A  detailed  discussion  of  phase-modulated  QPM  design  is  given  in  Sec.  7.2. 


CHAPTER  7.  ADVANCED  QFC  DEVICES 


98 


Figure  7.2:  Photon-flux  evolution  in  a  two-component  grating  optimized  for  efficient 
cascaded  npconversion. 


sequential  grating  approach  solved  in  Eqs.  (7.2),  it  is  interesting  to  consider  more 
sophisticated  QPM  designs  in  which  gA{z)  and  gs {z)  are  spatially  varying.  A  useful 
model  function  that  can  be  used  to  smoothly  interpolate  between  the  sequential-  and 
simultaneous-grating  approaches  is  given  as: 


9a{z)  =  2-“ 


1  —  tanh 


{  z-  La 

V  w 


(7.3) 


With  the  parameters  a  =  0.33  and  W  =  0.183P,  we  find  that  complete  conversion 
is  achieved  at  a  pump  power  22.1%  lower  than  for  sequential  gratings  and  28.9% 
lower  than  for  simultaneous  gratings.  The  photon  flux  evolution  and  grating  weight 
\gA(z)\2  for  an  optimized  grating  is  shown  in  Fig.  7.2.  In  future  work  we  may  pursue 
a  more  sophisticated  grating  design,  as  any  reduction  of  Pmax  also  yields  a  reduction 
in  noise.  We  have  not  proven  that  Eq.  (7.3)  produces  a  global  minimum  Pmax,  so 
further  reductions  in  Pmax  may  be  possible. 
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7.1.2  Cascaded  upconversion  for  high-speed  single-photon 
detection 

The  SPD  figure-of- merit  H  =  r/^/ (D  At)  for  the  upconversion  detector  described 
in  Chapter  5  was  limited  by  the  large  value  of  At  &  400  ps  of  the  thick-junction 
APD  used  in  the  experiment.  A  thick-junction  APD  was  used  because  the  thin- 
j unction  devices  are  optimized  for  blueshifted  wavelengths,  as  shown  in  Fig.  7.1(c). 
To  upconvert  a  signal  to  a  near-peak  wavelength  (close  to  600  nm)  in  a  single  step 
would  result  in  the  use  of  a  short-wavelength  pump,  which  results  in  high  noise: 
Tliew  et  al.  observed  a  noise  count  rate  of  3  x  105  counts/s  in  an  early  experiment 
pumped  at  980  nm  [33].  To  bypass  this  problem  we  use  the  cascaded  frequency 
conversion  approach.  A  long  wavelength  pump  at  ojp  is  summed  with  the  input  signal 
at  uj\  to  produce  radiation  at  uj-int  —  uj\  +  ojp.  A  second  SFG  process  (SFG-B)  is 
used  to  sum  the  light  at  wint  with  the  same  pump  to  produce  the  target  radiation  at 
0J2  =  oJint  +  (jJp  =  +  2cOp.  For  an  input  signal  at  Ai  =  1.55  /irn  and  Xp  =  1.8  /an, 

~  830  nm  and  A2  ~  570  nm.  This  two-step  process  produces  target  radiation 
near  the  peak  PDE  for  low-jitter  APDs. 

We  fabricated  a  two-step  upconverter  based  on  the  sequential  QPM  grating  ap¬ 
proach  described  above.  The  waveguide  had  total  length  L  =  5.2  cm,  and  incor¬ 
porated  mode  filters  and  adiabatic  tapers  to  spatially  mode-match  the  input  of  the 
waveguide  to  SMF-28  optical  fiber.  The  waveguide  was  fiber-pigtailcd  with  coupling 
loss  of  0.5  dB,  and  had  propagation  loss  of  0.12  dB  cnr1.  The  schematic  of  our  exper¬ 
iment  is  shown  in  Fig.  7.3.  Pump  light  at  1801  nm  from  a  monolithic  PPLN  optical 
parametric  oscillator  [97],  amplified  with  a  Tm-doped  fiber  amplifier,  was  combined 
with  a  1.55-/U1U  signal  using  a  fiber-optic  WDM.  The  waveguide  had  QPM  periods 
A^  =  18.55  /im  and  A#  =  8.25  /im,  with  La/ L  =  0.64.  For  initial  characterization  of 
the  devices,  the  signal  source  was  a  cw  tunable  external- cavity  diode  laser. 

We  first  characterized  the  conversion  efficiency  of  the  device  versus  signal  wave¬ 
length  and  temperature.  In  order  to  achieve  efficient  conversion,  both  SFG-A  and 
SFG-B  must  be  tuned  near  their  QPM  peaks.  This  tuning  was  done  by  noticing  that 
SFG-A  and  SFG-B  temperature-tune  at  slightly  different  rates.  Fig  7.4(a)  shows  the 
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Figure  7.3:  Experimental  setup  for  cascaded  upconversion  of  single  photons.  Abbre¬ 
viations:  PC,  polarization  controller;  LPF,  long-pass  filter;  VATT,  variable  optical 
attenuator;  MLL,  mode-locked  laser;  BPF,  band-pass  filter. 

experimental  results:  for  a  pump  wavelength  Xp  =  1801.3  nm,  quasi-phasematching 
for  both  SFG  processes  was  observed  for  Ai  =  1548.2  nm  and  T  =  76°  C. 

We  next  measured  the  conversion  efficiency  of  a  weak  classical  input  signal  at  uj i. 
When  SFG-A  and  SFG-B  are  both  tuned  to  their  QPM  peaks,  the  conversion  effi¬ 
ciencies  r/2  and  r][nt  are  plotted  in  Fig.  7.5,  along  with  a  fit  to  Eq.  (7.2).  Experimental 
values  of  the  conversion  efficiency  are  defined  with  respect  to  the  signal  power  at 
the  input  of  the  waveguide.  The  efficiency  for  the  cascaded  process  is  limited  by  the 
observed  propagation  loss  and  by  a  slightly  non-optimal  value  of  L a  for  the  device. 

We  demonstrated  an  application  of  this  device  for  high  clock-rate  single-photon 
detection  of  a  pulsed  light  source.  A  pulsed  1550-nm  signal  was  generated  using  a 
10-GHz  actively  mode-locked  laser  with  pulse  length  of  10  ps,  and  using  an  electro¬ 
optic  intensity  modulator  to  downsample  the  pulse  train  to  1  GHz.  This  signal 
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Figure  7.4:  Measured  cascaded  SFG  at  A2  =  570  nm  as  a  function  of  signal  wavelength 
Ai  and  device  temperature  T,  showing  dual  QPM  at  Ai  =  1548.2  nm  and  T  =  76°  C. 
Green  and  red  curves  are  the  individual  QPM  trajectories  for  SFG-A  and  SFG-B 
respectively;  the  efficiency  of  the  two-step  process  is  maximized  where  these  curves 
intersect. 

was  attenuated  to  a  level  of  0.01  photons  per  pulse  and  upconverted  using  the  two- 
step  upconversion  waveguide  described  above.  The  output  signal  at  0J2  was  detected 
using  a  Si  APD  (id  Quantique  idlOO).  Measurements  were  performed  with  Pp  = 
100  mW,  corresponding  to  a  background  dark-count  rate  of  approximately  850  counts 
s"1  which  was  attributed  to  two-step  upconversion  of  noise  photons  generated  by 
anti-Stokes  spontaneous  Raman  scattering  of  the  1.8-/im  pump  [67].  Time-correlated 
single-photon  counting  measurements  were  carried  out  using  a  time-interval  analyzer 
(TIA)  started  by  the  clock  output  of  the  pattern  generator  and  stopped  by  the  APD 
output  electrical  pulses.  Our  results  are  plotted  in  Fig.  7.6:  we  observe  a  FWHM 
timing  jitter  of  66.9  ps,  which  is  close  to  the  50-ps  specification  provided  by  the 
manufacturer.  Additional  sources  of  timing  uncertainty  are  due  to  jitter  within  the 
TIA  of  approximately  25  ps,  and  noise  in  the  triggering  electronics.  We  do  not  suspect 
the  introduction  of  additional  timing  jitter  by  the  upconversion  process. 

To  conclude  this  section,  we  have  demonstrated  87%  conversion  efficiency  of  a 
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Figure  7.5:  Measured  internal  conversion  efficiencies  r]2  and  r]m t  versus  pump  power 
Pp  and  fit  to  Eq.  (7.2). 

telecom-band  signal  to  the  green  spectral  region  via  a  two-step  upconversion  process 
using  a  single  1.8- /mi  pump.  Using  a  picosecond  pulsed  single-photon-level  signal  we 
have  demonstrated  single-photon  detection  timing  resolution  below  70  ps,  which  may 
enable  higher-data-rate  quantum  communications  and  quantum  key  distribution.  If, 
using  a  narrow  band  filter  such  as  a  volume  Bragg  grating  (described  in  Sec.  7.2), 
and  the  state-of-the-art  PDM  series  Si  SPAD,  we  could  achieve  1.55  /mi  single-photon 
detection  with  system  efficiencies  of  up  to  35%  at  noise  count  levels  below  100  s_1 
and  timing  jitter  below  40  ps,  potentially  enabling  QKD  with  modulation  rates  of 
greater  than  10  GHz  over  distances  longer  than  250  km. 

7.1.3  Two-step  downconversion  of  quantum  light  sources 

In  Chapter  6,  we  discussed  frequency  conversion  devices  for  downconverting  single¬ 
photons  from  InAs/InGaAs  quantum  dots  from  their  emission  wavelength  near  A2  = 
900  nm  to  a  target  wavelength  Ai  in  the  telecommunications  band  near  1550  nm.  By 
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Figure  7.6:  Timing  histogram  of  10-ps-pulsed  1.55-/mi  signal  using  two-step  up  con¬ 
verter  and  idlOO  Si  APD,  showing  system  timing  jitter  of  66.9  ps. 

energy  conservation,  this  process  requires  a  strong  pump  wave  at  up  =  uj2  —  uji,  which 
yields  Xp  =  2.15  /im.  If  one  wishes  to  use  a  pump  frequency  up  <  cni,  one  is  therefore 
limited  to  input  signal  frequencies  uj2  <  2cj1  . 

A  great  deal  of  research  is  being  devoted  to  quantum  emitters  at  frequencies  in  the 
visible  and  near- visible  that  do  not  satisfy  this  constraint.  One  notable  example  is  the 
diamond  nitrogen-vacancy  (NV)  color  center,  an  intrinsic  diamond  impurity  which 
when  localized  is  a  single  quantum  system  [140].  There  has  been  significant  progress 
in  fabricating  nanostructures  to  enhance  the  radiative  properties  of  NV  centers  [141, 
142],  as  well  as  demonstrations  of  spin-photon  entanglement  [4]  and  coherent  control 
of  single  spins  [143].  Less  well-developed  but  still  promising  defect  centers  include 
silicon- vacancy  as  well  as  Cr  defect  centers,  all  of  which  operate  near  700  nm  [140]. 
However,  if  one  seeks  to  couple  these  sources  to  the  telecommunications  infrastructure, 
it  requires  a  more  sophisticated  frequency  conversion  scheme. 

To  this  end,  we  propose  the  use  of  a  cascaded  downconversion  scheme,  first  dis¬ 
cussed  in  [48],  and  illustrated  schematically  in  Fig.  7.7(a).  When  the  input  signal 
is  at  a  frequency  co2  >  2cni,  we  consider  a  downconversion  process  in  which  a  single 
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Figure  7.7:  (a)  Schematic  of  frequencies  involved  in  two-stage  downconversion  of 
single  photons  from  a  diamond  NV  center  to  the  1.5-/im  telecom  band;  (b)  schematic 
device  design  for  two-stage  downconversion. 

pump  at  Up  is  used  to  convert  the  input  signal  to  the  target  via  an  intermediate  held 
at  cCjnt .  From  Fig.  7.7,  we  therefore  have  up  =  (u2  —  u\)/2  and  u-irit  =  (cui  +  U2)/2.  A 
schematic  PPLN  waveguide  for  this  process  is  shown  in  Fig.  7.7(b).  The  device  will 
use  an  integrated  directional  coupler  as  discussed  in  Chapter  6  to  combine  the  long- 
wavelength  pump  and  short-wavelength  single  photons  into  the  fundamental  spatial 
mode  of  the  mixing  waveguide,  where  two  sequential  DFG  processes  will  convert  the 
signal  photons  at  637  nm  to  target  photons  at  1550  nm. 

While  diamond  NV  centers  are  touted  as  potential  room-temperature  qubits,  the 
photon  emission  from  NV  centers  is  not  ideal  for  frequency  conversion  applications. 
Even  at  cryogenic  temperatures,  NV-center  emission  is  dominated  by  a  phonon  side¬ 
band  with  a  width  of  approximately  100  nm  centered  around  700  nm.  Zero-phonon 
emission  occurs  at  a  wavelength  of  637  nm  (for  NV-)  and  570  nm  (for  NV°)  but  even 
at  10  K  the  zero-phonon  emission  fraction  is  approximately  3%  [141],  At  cryogenic 
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temperatures,  the  linewidths  of  the  zero-phonon  lines  (ZPLs)  are  very  narrow,  suit¬ 
able  for  frequency  conversion,  but  near  room  temperature  they  broaden  to  several  nm, 
reducing  their  overlap  with  the  acceptance  bandwidth  of  a  down-conversion  device. 
Nevertheless,  the  fabrication  of  photonic  nanostructures  in  single-crystal  diamond  has 
enabled  a  strong  Purcell  enhancement  of  the  spontaneous  emission  rate  [142],  with 
up  to  70%  of  the  emitted  photons  in  the  ZPL.  These  recent  developments  make  the 
implementation  of  a  cascaded  frequency  downconverter  very  promising. 

7.2  Multichannel  and  time- frequency  converters 

In  Sec.  7.1.2  we  discussed  the  use  of  cascaded  frequency  upconversion  to  enable  detec¬ 
tion  by  low-timing-jitter  APDs.  Here,  we  explore  the  utility  of  using  multiple  pumps, 
either  pulsed  or  cw,  to  upconvert  an  optical  signal.  Experiments  in  this  section  have 
been  done  in  collaboration  with  Paulina  Kuo  and  Xiao  Tang  of  the  Information  Tech¬ 
nology  Laboratory  at  NIST  [144], 

7.2.1  Multichannel  upconversion  and  applications 

We  consider  a  waveguide  in  which  the  following  SFG  processes  can  be  simultaneously 
quasi-phasematched: 


(7.4a) 

(7.4b) 


Uq  +  UJpi  —  uq 
U>0  +  UJP2  =  UJ2 


(7.4c) 


Wo  +  UJpN  —  Wat, 


where,  as  in  Fig.  7.8,  we  imagine  that  a  single-photon  signal  at  cuo  is  undergoing 
simultaneous  sum-frequency  generation  with  N  strong  pumps  at  frequencies  ojp\  to 
u>pN-  Considering  now  the  simplest  (non-trivial)  case  of  N  =  2,  the  coupled  wave 
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equations  for  the  photon  fields  (normalized  as  in  Chapter  2),  are  written: 


da0 

d^ 

dax 

dT 

da2 

dT 


=  -7lfl2  -  72^3 


=  7i«! 


=  72^1 


(7.5a) 

(7.5b) 

(7.5c) 


where  the  7 j  =  ^nPPj.  The  pumps  at  ojpi  and  cop 2  are  assumed  to  be  undepleted 
(valid  for  the  weak  input  signals  we  consider),  the  processes  are  both  phasematched 
(which  can  be  achieved  using  a  phase- modulated  QPM  grating  described  below),  and 
the  propagation  losses  are  negligible. 

We  seek  to  solve  the  Eqs.  (7.5a)  and  determine  the  input-output  relation  for  the 
photon  annihilation  operators  at  the  output  of  the  device.  We  start  by  taking  the 
second  derivative 

=  -(7i2  +  7f)ao-  (7.6) 

Substituting  T  =  (7^  +  7I)1/2,  this  equation  is  solved  to  yield 


a0(z)  =  AcosIW  +  B  sinT^. 


(7.7) 


We  proceed  to  solve  the  remaining  two  equations  and  put  in  the  initial  conditions 
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a0(0),  ai(0),  and  a2(0)  to  yield  the  input-output  (z  =  L )  relation.  Defining  normalized 
coupling  constants  7 j  =  Jj/T,  we  find: 


a0(L )  =  a0(0)  cosTL  -  [7101  (0)  +  72a2(0)]  sinTL, 


(7.8a) 


ai(L)  =  7ia0(0)  sin  T L  +  07 (0)  (7^  cosTL  +  yf)  +  7i72a2(0)  (cosTL  -  1) ,  (7.8b) 
a2(L)  =  72a0(0)  smTL  +  a2(0)  (72  cosTL  +  7^)  +  7i72ai(0)  (cosTL  -  1) .  (7.8c) 

The  input-output  relations  of  Eq.  (7.8a)  can  be  easily  generalized  where  we  note  that 
for  N  pumps,  T  =  ( JV  7 |)1//2;  where  JV7?  =  1-  For  an  ./V-pump  system,  we  can 
create  probabilistic  mixtures  of  photons  at  the  target  frequencies  07  to  u>jy. 

QPM  design  and  device  characterstics 

To  fabricate  a  multichannel  reconversion  device,  we  employed  the  QPM  phase- 
modulation  technique  developed  by  Asobe  and  coworkers  [138].  The  phase  modu¬ 
lation  technique  is  one  of  several  techniques  that  have  been  developed  in  order  to 
engineer  the  phasematching  profile  of  a  device.  In  this  technique,  one  applies  a  su¬ 
perstructure  to  the  QPM  grating  by  modulating  the  phase  of  the  QPM  pattern  using 
a  phase  modulation  function  cp(z)  that  is  periodic  in  z  with  period  Aph-  The  phase 
modulation  consist  of  shifting  position  of  a  domain  boundary  zn  with  respect  to  its 
nominal  location  zn  0  of  a  periodic  grating.  In  this  formalism,  the  modulation  phase 
4>(zn)  =  27i(zn  —  zn<a)/ Ag>.  Via  Fourier  analysis,  it  can  be  shown  that  one  obtains 
QPM  peaks  when  the  relation 


(7.9) 


is  satisfied  for  integer  m.  To  choose  the  appropriate  phase  modulation  function  (p(z) 
we  minimize  a  trial  function  given  as 


T = E  -  i"‘)2 


mgS 


(7.10) 
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Figure  7.9:  (a)  Theoretical  and  (b)  measured  SFG  efficiency  for  phase-modulated 
SFG  device  for  A0  =  1319  nm. 

where  the  sum  is  over  integer  m  in  the  set  S  containing  the  desired  QPM  peak 
locations,  r/m  is  the  computed  efficiency  for  a  particular  (p(z),  and  r/^rg  is  the  desired 
efficiency.  For  the  two  channel  device  we  fabricated,  S  =  {—1,  1},  Vm=± i  =  0-5,  and 
rfmg  =  0  f°r  all  other  m.  The  use  of  a  simplex  convex  optimization  algorithm  (the 
fminsearch  function  in  Matlab)  was  successful  at  optimizing  <fi(z)  for  several  different 
functional  parametrizations  (e.g.  Fourier  series  and  power  series).  The  simple  two- 
channel  structure  did  not  require  such  a  complex  optimization  technique,  but  will 
enable  the  design  of  more  complex  devices  (including  asymmetric  about  m  =  0,  and 
unequal  peak-height)  structures  with  more  channels  [145].  A  calculated  QPM  tuning 
curve  is  shown  in  Fig.  7.9(a). 

The  channel  separation  (which  we  define  as  the  wavelength  difference  between 
the  two  QPM  peaks  for  m  =  ±1)  is  controlled  via  a  choice  of  the  phase  modulation 
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Figure  7.10:  Required  phase  modulation  period  Aph  for  SFG  between  two  pumps  with 
wavelengths  spaced  by  AA  around  1550  nm  and  a  1310-nm  signal  (circles)  with  power 
law  Aph  =  a(AA)-1  fit. 


period.  The  required  phase-modulation  period  for  a  two-channel  device  is  found  by 
rearranging  Eq.  (7.9) 


1  _  Ak'(AX)  1 

Aph  27 r  A  g 


(7.11) 


We  plot  the  result  of  Eq.  (7.11)  in  Fig.  7.10.  The  A k'  values  were  computed  using  the 
simulated  dispersion  data  for  RPE  waveguides  with  1.85  /irn  proton  exchange  depth 
(circles),  and  are  found  to  be  fit  accurately  by  a  power  law  Aph  =  a(AA)-1,  with 
a  =  37.4  mm/nm.  For  our  devices  we  chose  target  channel  separations  of  5,  10,  and 
15  nm,  which  implied  phase  modulation  periods  of  7.1,  3.94,  and  2.56  mm. 

We  fabricated  two-channel  upconversion  RPE  PPLN  waveguides  for  SFG  between 
a  1.3  /irn  signal  and  two  1.55-/im-band  pumps  using  the  techniques  discussed  in 
Chapter  3.  The  devices  were  antireflection  coated  (Foreal,  Inc.)  for  wavelengths  of 
1310  and  1550  at  the  input  between  LiNbOs  and  optical  fiber,  and  for  1310,  1550, 
and  710  nm  at  the  output  between  LiNb03  and  air.  The  devices  were  fiber  pigtailed 
as  well,  and  were  measured  to  have  total  optical  throughput  of  -1.3  dB  at  1550  nm 
and  -2.4  dB  at  1319  nm. 

For  preliminary  experiments  described  below,  we  chose  a  device  with  AA  =  15  nm. 
A  QPM  tuning  curve  is  shown  in  Fig.  7.9(b),  where  we  can  see  excellent  agreement 
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Figure  7.11:  Experimental  setup  for  multichannel  upconversion  system. 


with  theory.  The  device  has  two  QPM  channels  with  phasematching  of  1532.88  and 
1547.43  nm. 

We  incorporated  the  waveguide  into  an  upconversion  single-photon  detector,  a 
schematic  of  which  is  shown  in  Fig.  7.11.  To  pump  both  upconversion  channels 
simultaneously,  two  pump  lasers  were  used.  As  will  be  described  in  Sec.  7.2.2,  the 
pumps  were  modulated  using  electro-optic  intensity  modulators  (IM).  The  two  pumps 
were  combined  using  a  C/L-band  WDM  and  were  amplified  using  an  EDFA.  Optical 
noise  in  the  1.3  /irn  band  coming  from  the  EDFA  was  filtered  using  a  series  of  three 
1550/1310-nm  WDMs.  After  passing  through  a  tap  coupler,  the  combined  dual¬ 
channel  pump  was  then  combined  with  the  single-photon  level  signal  at  1.3  /irn  with 
another  1550/1310  WDM,  and  the  combined  pump  and  signal  radiation  was  incident 
on  the  fiber  pigtail.  Polarization  controllers  were  used  for  both  pumps  and  the  signal 
to  align  the  input  light  to  the  TM  polarization  of  the  waveguide. 
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Following  the  waveguide,  an  aspheric  lens  (AL)  coated  was  used  to  collimate  the 
output  radiation  from  the  waveguide.  To  separate  the  outgoing  pump  light  and  filter 
the  two  upconverted  channels  at  Ai  and  A2,  we  used  a  reflective  volume  Bragg  grating 
(OptiGrate,  Inc.).  The  reflective  VBG  acts  as  a  narrowband  filter  (AA  <  0.05  nm) 
where  the  reflected  band  depends  on  the  angle  of  incidence.  The  VBG  was  aligned 
such  that  the  long-wavelength  channel  at  A2  had  high  reflectivity  back  through  the 
input  facet.  The  upconverted  light  at  Ai  passed  through  the  VBG  and  was  directed 
back  in  using  a  mirror  at  a  slight  angle  so  that  high  reflection  was  achieved  back 
through  the  VBG.  The  two  upconverted  channels  were  then  directed  to  two  Perkin- 
Elmer  SPCMs,  where  a  20-nm  bandpass  filter  at  the  entrance  of  the  SPCMs  blocked 
any  stray  light  or  parasitic  pump  second  harmonic  radiation2 *.  Not  shown  in  Fig.  7.11 
are  reimaging  systems  for  each  channel  used  to  adjust  the  spot  size  incident  on  the 
aspheric  lenses  before  the  APDs.  The  VBG  performance  depends  critically  on  the 
focusing  of  the  incident  beam  (the  reflectance  is  maximized  for  a  perfectly  collimated 
beam)  and  we  thus  needed  an  additional  focusing  degree  of  freedom  when  compared 
with  the  upconversion  detectors  discussed  in  Chapter  5. 

We  next  measured  the  system  PDE  and  NCR  of  each  channel,  separately.  The  re¬ 
sults  are  shown  in  Fig.  7.12,  with  fits  of  the  PDE  to  r]tot(Pp)  =  rj0  sin2 (zr/2 \JPp/ -Pmax)  • 
We  find  a  maximum  PDE  of  32%  for  channel  2  and  25%  for  channel  1.  The  PDE  for 
channel  1  is  lower  due  the  additional  optical  path  needed  as  it  is  the  channel  trans¬ 
mitted  through  the  VBG.  A  major  improvement  of  these  results  when  compared  with 
earlier  1.3  /im-band  upconversion  detectors  is  the  very  low  noise  count  rate-here,  just 
barely  above  the  intrinsic  dark  count  rate  of  the  SPCMs,  compared  with  greater  than 
104  counts/s  in  a  previous  experiment  with  a  PPLN  waveguide  [30].  This  reduction 
in  the  noise  is  due  to  the  very  narrowband  reflectance  of  the  VBGs.  I11  this  system, 
the  noise  is  primarily  due  to  upconverted  anti-Stokes  Raman  photons  generated  by 
the  strong  pump,  which  fill  the  acceptance  bandwidth  of  the  upconverter.  Here,  the 

2  An  unintended  consequence  of  the  phase- modulated  QPM  design  was  the  appearance  of  several 

parasitic  accidental  QPM  peaks  for  SHG  of  the  pump  radiation.  We  chose  the  operating  wavelengths 

such  that  these  peaks  were  avoided,  but  even  away  from  these  peaks  there  were  still  up  to  a  few  tens 
of  microwatts  of  generated  second  harmonic  radiation  when  pumping  with  approximately  100  mW 
of  pump  light. 
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Figure  7.12:  Measured  two-channel  upconversion  detector  (a)  PDE,  and  (b)  noise 
count  rate. 

VBG  provides  a  filter  narrower  than  the  acceptance  bandwidth  such  that  most  of  the 
noise  photons  are  eliminated. 

Applications  as  a  frequency-shifting  variable  beam  splitter 

When  both  pump  frequencies  are  active  at  the  same  time,  the  dual-channel  up  con¬ 
verter  followed  by  spectral  filtering  acts  as  a  frequency  shifting  variable  beamsplitter. 
Beamsplitters  have  several  uses  in  quantum  optics;  we  demonstrated  the  use  of  the 
dual-channel  upconverter  for  two  purposes:  (1)  the  splitting  of  high- count-rate  sig¬ 
nals  to  enable  single-photon  count  rates  above  the  dead-time-limited  value  of  a  single 
detector,  and  (2)  the  measurement  of  the  second-order  correlation  function  g^2'{r)  of 
both  coherent  and  quasi-thermal  light  sources. 

As  discussed  in  Sec.  5.1,  the  maximum  count  rate  of  an  APD  is  determined  by 
the  detector  dead  time  7p  as  Rmax  =  r p1.  As  the  incident  photon  rate  R.inc  begins 
to  approach  i?max  the  detector  begins  to  saturate  at  this  maximum  count  rate  (the 
details  of  the  saturation  depend  on  the  statistics  of  the  input  light).  We  consider  a 
single  photon  at  coq  incident  at  the  input  of  the  device.  When  the  pump  powers  Pp\ 
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Figure  7.13:  Measured  count  rate  R  versus  input  signal  photon  flux  R[nc  for  a  single- 
channel  (blue  circles)  or  dual-channel  (green  squares)  device,  where  we  observe  a 
saturation  rate  of  2Rmax  for  the  two-channel  detector.  The  solid  curves  are  simulated 
counting  statistics  for  coherent  light. 

and  Pp2  are  chosen  such  that  YL  =  7t/2  and  7i  =  72,  we  find,  with  50%  probability 
that  photon  is  converted  to  07  or  uj2-  Therefore,  much  like  a  classical  beamsplitter, 
the  dual-channel  upconverter  will  probabilistically  divide  the  incident  1 . 3- /.irn  signal 
to  the  two  upconverted  channels  where  the  signal  can  be  detected  by  two  APDs. 

We  have  measured  the  saturation  of  a  Perkin-Elmer  SPCM  by  varying  Rinc  and 
measuring  the  count  rate  R  of  a  single  detector.  The  input  signal  i?inc  at  Ao  =  1302  nm 
was  controlled  via  a  programmable  optical  attenuator  from  a  level  of  41  photons/s 
to  2  x  109  photons/s.  With  only  pump  1  turned  on,  we  measured  the  count  rate  Rj 
on  both  channels  j  G  {1,2}:  the  count  rate  R\  is  shown  as  blue  circles  in  Fig.  7.13. 
The  solid  blue  curve  is  a  simulation  of  the  expected  counting  statistics  of  coherent 
light  (with  Poissonian  photon  arrival  statistics)  assuming  a  detector  dark  count  rate 
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D  =  250  s~4  and  a  dead  time  of  67  ns.  By  comparing  the  rates  R\  and  with  only 
pump  1  turned  on,  we  measured  a  channel-to- channel  crosstalk  of  -44  dB,  which  is 
likely  due  to  a  combination  of  amplified  spontaneous  emission  (ASE)  from  the  EDFA 
and  imperfect  isolation  of  the  VBG  channels. 

With  both  pumps  turned  on,  we  measured  the  summed  count  rate  Rs  =  R\  +  i?2 
after  equalizing  the  count  rates  on  each  detector,  and  at  the  same  overall  conversion 
efficiency  as  in  the  pump-l-only  case.  The  observed  Rs  versus  Ainc  is  plotted  as  green 
squares  in  Fig.  7.13,  along  with  an  analogous  simulation  of  the  counting  statistics. 
We  see  that  at  low  Rmc,  the  dual-channel  case  has  a  DCR  of  approximately  2 D,  and 
saturates  at  high  Rs  =  2 i?max.  From  this  result,  we  can  extrapolate  that  for  an  N- 
pump  upconverter,  the  system  count  rate  Rs  would  saturate  at  NRmax  at  the  expense 
of  a  DCR  of  ND. 

A  common  application  of  beamsplitters  in  quantum  optics  is  to  split  a  signal  to 
two  detectors  such  that  the  second-order  correlation  function  (r)  of  a  light  source 
can  be  measured  (see  Chapter  6  for  discussion).  We  have  used  the  dual-channel 
converter  as  an  upconverting  beamsplitter  and  have  measured  the  g^2\r)  of  both  a 
coherent  source  and  a  pseudothermal  source.  In  the  measurement,  we  used  the  time- 
tagging  mode  of  a  time-correlated  single-photon  counting  (TCSPC)  system  (PicoHarp 
300E)  to  record  the  arrival  times  of  photons  at  APD1  and  APD2  where  the  pump 
powers  Ppi  and  Pp 2  were  chosen  to  equalize  the  two  count  rates.  The  g^2\r)  was 
then  calculated  offline  upon  completion  of  a  1-minute  integration  at  approximately 
5  x  104  counts/s/channel.  The  results  for  the  coherent  source  are  shown  in  Fig.  7.14. 
As  expected,  due  to  the  Poissonian  arrival  statistics  of  the  photons  from  a  laser, 

9('2\t)  —  1  f°r  all  r- 

It  has  been  shown  that  by  passing  a  laser  beam  through  (or  reflecting  a  laser 
beam  off)  a  moving  diffuser,  one  can  create  pseudothermal  light  source  with  a  vari¬ 
able  correlation  time  [146,  147].  We  formed  a  pseudothermal  source  by  coupling  the 
signal  laser  into  free  space  and  collecting  the  scattered  light  off  a  rotating  ground- 
glass  diffuser.  By  changing  the  geometry  of  the  collection  or  the  rotation  speed  of  the 
diffuser  different  coherence  times  At  can  be  realized.  Our  results  for  the  pseudother¬ 
mal  source  are  shown  in  Fig.  7.14(b),  where  we  see  a  bunching  peak  with  a  FWHM 


CHAPTER  7.  ADVANCED  QFC  DEVICES 


115 


Delay  x  (ps) 


Delay  x  (ps) 


Figure  7.14:  Measured  second-order  correlation  functions  g^2\r)  of  (a)  coherent  and 
(b)  quasi-thermal  light  sources  using  dual-channel  upconverter. 


width  of  98  ps. 

The  results  obtained  by  measuring  a  g^ij)  using  a  dual-channel  upconverter 
are  presumably  no  different  than  would  be  obtained  via  a  standard  technique  using 
a  classical  beam-splitter.  One  benefit  is  the  ability  to  easily  and  rapidly  change 
the  “splitting  ratio”  via  control  over  the  pumping  ratios  71  and  72.  That  we  have 
observed  the  possibility  of  probabilistically  upconverting  to  one  state  or  another  has 
interesting  consequences  when  we  imagine  manipulating  nonclassical  states  of  light 
with  this  device.  If  a  single  photon  at  cuo  is  incident  on  the  device,  we  can  write  the 
input  state  as 


|^)in  =  1 100)  =  |  1)oj0  ®  |0)W1  <g>  |0)w2 

=  a0(0)t|000).  (7.12) 

After  propagating  through  the  dual-channel  device,  using  the  results  calculated  in 
Eq.  (7.8a)  we  find  that  the  output  state 


|-0)out  =  ao(-h)t|V'in)  =  7i  1 010)  +  e^rel7|001) 


(7.13) 
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which  implies  that  the  photon  has  been  converted  to  a  superposition  state  of  two  dif¬ 
ferent  frequencies.  If  the  relative  phase  (f)re\  between  the  two  pumps  can  be  controlled, 
coherent  superpositions  of  the  photon  can  be  created  in  Ai  and  A2,  which  could  be  a 
novel  qubit  in  quantum  information  processing  systems. 

7.2.2  Upconversion  with  multi-wavelength  sampling 

In  the  previous  section  we  discussed  applications  of  a  multichannel  converter  with 
the  pump  sources  on  simultaneously.  In  that  demonstration,  the  splitting  of  an  input 
signal  to  multiple  output  ports  happened  stochastically.  If,  however,  the  pumps  were 
turned  on  and  off  via  external  modulators  controlled  by  a  pattern  generator,  one  can 
deterministically  control  the  output  channel  via,  for  example,  switching  the  pumps 
on  and  off  sequentially.  An  upconversion  detector  using  multi-wavelength  sampling 
techniques  has  recently  been  demonstrated  by  Ma  and  coworkers  at  NIST  [148]. 

The  primary  motivation  for  a  multi-wavelength  sampling  upconversion  detector 
is  to  enable  the  detection  of  single-photon-level  signals  modulated  at  rates  higher 
than  the  timing  jitter  At  of  a  given  detector  would  allow.  The  timing  jitter  sets  the 
maximum  clock  rate  for  a  QKD  system,  and  with  APD  single-photon  detectors  the 
long  tails  of  the  timing  response  (see  Fig.  7.16(a))  would  lead  to  high  intersymbol 
interference  when  a  system  is  clocked  at  rates  above  fc  ~  1/At  «  1  GHz.  The  use  of 
N  switched  distinctly  colored  pumps  enables  the  subdivision  of  this  1-ns  bin  into  N 
sub-bins,  enabling  a  time-division  demultiplexing  of  a  signal  that  can  now  be  clocked 
at  Nfc. 

Our  experimental  setup  is  diagrammed  schematically  in  Fig.  7.11.  A  pulse-pattern 
generator  (PPG,  Tektronix  DTG-5274)  drives  electro-optic  intensity  modulators  (IM) 
for  the  signal  and  each  pump.  The  timing  diagram  for  the  experiment  is  shown 
Fig.  7.15.  After  combining  the  two  pumps  and  signal  in  the  waveguide,  the  resulting 
dual-channel  SFG  products  are  separated  and  filtered  by  the  VBG  and  routed  to 
two  Si  APDs  (Perkin  Elmer  SPCM-AQR14).  The  wide  channel  separation  would 
make  possible  a  filtering  scheme  either  using  a  holographic  diffraction  grating  (as  in 
[148],  but  with  significantly  shorter  optical  path  required),  or  a  Brewster  prism.  One 
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Figure  7.15:  Timing  diagram  for  detection  of  high-clock-rate  signals  via  multi¬ 
wavelength  sampling. 


difficulty  about  working  with  the  large  channel  separation  is  that  the  shape  of  the 
EDFA  gain  spectrum  must  be  carefully  considered  to  balance  the  conversion  efficiency 
between  the  two  channels. 

We  first  measured  the  timing  jitter  of  the  APD  used  in  the  experiment,  by  modu¬ 
lating  the  signal  into  a  short  pulse  every  10  clock  cycles,  and  using  the  pulsed  pump 
to  upconvert  it.  Our  data  is  shown  in  Fig.  7.16(a)  as  green  dots,  where  we  observe  a 
FWHM  jitter  of  380  ps,  and  FW1%M  timing  jitter  of  1.05  ns.  The  primary  imped¬ 
iment  to  using  this  detector  for  quantum  communication  schemes  is  the  long  timing 
jitter  tail,  which  is  due  to  carrier  diffusion  effects  [149].  Using  a  single  pump  modu¬ 
lated  at  the  same  rate  as  the  signal  (1.6  GHz),  we  recorded  a  timing  histogram  which 
is  shown  as  the  black  dots  in  Fig.  7.16(a).  The  background  level  of  approximately 
1CD3  normalized  counts  is  due  to  residual  pump  and  signal  transmission  due  to  im¬ 
perfect  extinction  of  the  EOMs.  Due  to  the  detector  timing  jitter,  a  photon  arriving 
within  one  clock  cycle  will  have  a  large  probability  of  producing  an  electrical  pulse 
that  will  register  in  an  adjacent  time  bin,  which  will  result  in  errors  in  the  quantum 
communication  protocol.  We  can  calculate  this  error  probability  by  comparing  the 
probability  of  assigning  a  count  to  the  wrong  time  bin  to  the  count  statistics  for 
a  single  pulse.  For  the  case  of  detecting  a  1.6-GHz-modulated  signal  with  a  single 
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(a)  Measured  counting  histograms  for  an  isolated  1302-nm  pulse  (green  dots)  showing 
a  long  tail  which  results  in  high  intersymbol  interference  for  a  1.6-GHz  clock  signal 
measured  using  only  one  channel  (black  dots). 
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(b)  Measured  counting  histograms  for  a  1.6-GHz  clock  signal  using  a  dual-channel  al¬ 
ternating  pump  scheme,  showing  significantly  improved  intersymbol  interference  com¬ 
pared  with  a  single-channel  detection  scheme. 


Figure  7.16:  Timing  performance  of  dual-channel  upconversion  detector. 
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detector,  we  find  an  error  probability  Pe  =  5.7%. 

We  next  implemented  a  modulation  scheme  like  that  shown  in  Fig.  7.15.  A  1302- 
nm  clock  signal  at  1.6  GHz  was  upconverted  via  an  alternating  pump  such  that  each 
detector  is  only  receiving  pulses  at  a  rate  of  0.8  GHz.  We  recorded  timing  histograms 
for  both  detectors,  the  results  of  which  are  shown  in  Fig.  7.16(b).  We  now  observe 
that  due  to  the  effects  of  downsampling,  the  error  rate  is  reduced  by  more  than  an 
order  of  magnitude:  we  calculate  a  Pe  of  0.38%  (0.41%)  for  Channel  1  (Channel  2). 
As  in  [148],  we  have  also  used  the  system  to  detect  complex  pulse  sequences. 

7.2.3  Future  work 

We  have  described  the  use  of  a  dual-channel  upconverter  for  the  purposes  of  either 
creating  frequency  superpositions  or  detecting  high-clock-rate  modulated  signals  for 
communications.  Our  experimental  efforts  have  focused  on  a  two-channel  upconverter 
made  via  phase-modulated  QPM  techniques.  Such  techniques  are  easily  adaptable  to 
more  channels  [145],  and,  as  discussed  in  [148],  with  N  pumps  one  can  increase  the 
system  clock  rate  to  N/ At. 

A  difficulty  when  considering  increasing  the  number  of  channels  in  a  multichannel 
upconversion  system  is  that  one  needs  ever  more  lasers  and  modulators  to  expand  the 
system.  A  more  modular  approach  may  be  possible  based  on  optical  pulse  shaping 
techniques  that  have  been  heavily  exploited  for  arbitrary  waveform  generation,  among 
other  application  [150].  In  this  technique,  a  high-bandwidth  pulse  train  (e.g.  fem¬ 
tosecond  pulses  from  a  modelocked  laser)  are  manipulated  in  the  Fourier  domain.  If 
the  modelocked  laser  pulses  have  been  carrier-envelope-phase  (CEP)  stabilized  [151], 
then  one  can  envision  the  creation  of  the  frequency-entangled  states  with  stable  phase 
relationships  between  the  different  components  of  the  superposition. 

The  multi-wavelength  sampling  upconverter  can  in  some  senses  be  viewed  as  a 
discrete  analogue  of  a  time-frequency  converter.  In  a  time-frequency  converter,  each 
temporal  slice  of  an  input  signal  is  linearly  mapped  to  a  specific  frequency.  This 
mapping  corresponds  to  the  imposition  of  a  quadratic  phase  in  time,  which  is  the 
space-time  analog  of  a  lens  and  is  the  central  component  of  temporal  imaging  systems 
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Figure  7.17:  Schematic  of  a  single-photon  time-frequency  converter  utilizing  SFG  of 
an  input  signal  with  a  chirped  pump  pulse  in  a  chirped  QPM  grating. 


[152],  Temporal  imaging  systems  have  most  often  been  used  to  magnify  an  ultrafast 
waveform  in  time,  such  that  it  can  be  recorded  with  a  photodiode  and  standard 
oscilloscope;  temporal  imaging  systems  with  magnifications  of  up  to  450  x  have  been 
reported  [153]. 

We  propose  a  single-photon  time-frequency  converter,  a  schematic  of  which  can 
be  seen  in  Fig.  7.17.  A  chirped  pump  (created,  e.g.  by  dispersing  a  femtosecond 
pulse  in  a  length  of  optical  fiber)  is  combined  with  an  weak  input  signal  with  rapid 
temporal  features.  A  SFG  interaction  in  a  chirped  waveguide  (with  QPM  bandwidth 
broad  enough  to  support  the  entire  chirped  pump)  maps  the  temporal  features  of 
the  input  signal  onto  frequency.  The  chief  difference  between  this  system  and  ordi¬ 
nary  temporal  imaging  systems  is  the  use  of  a  diffractive  optic  and  detector  array 
to  detect  the  upconverted  signal.  This  configuration  allows  one  to  map  each  array 
element  onto  a  particular  slice  in  time,  such  that  even  if  each  detector  element  has 
a  timing  jitter  large  compared  to  the  temporal  features  of  the  signal,  the  features 
can  be  recorded  statistically  by  recording  many  shots.  SPAD  arrays  with  as  many 
as  128x128  elements  have  been  reported  [154],  We  believe  this  system  will  enable 
quantum  communications  systems  to  begin  to  approach  the  data  rates  of  classical 
communications  systems:  the  main  impediment  of  detector  timing  jitter  is  overcome 
by  the  effective  time-to-space  conversion  provided  by  this  system. 


Chapter  8 
Conclusion 


8.1  Summary  of  contributions 

The  main  theme  of  this  dissertation  has  been  the  development,  application,  and 
analysis  of  devices  and  systems  for  frequency  conversion  of  single-photon  signals. 
Here  I  shall  give  an  overview  of  the  main  contributions  of  the  work  described  in  this 
dissertation,  as  well  as  a  few  projects  to  which  I  have  contributed  that  are  not  covered 
in  this  dissertation. 

The  analysis  of  noise  processes  in  quantum  frequency  converters  was  a  principal 
subject  of  the  work  described.  As  discussed  in  Chapter  4,  there  are  two  primary  noise 
sources  in  QFC  devices.  The  first  source  is  spontaneous  parametric  fluorescence, 
whereby  a  pump  photon  splits  into  a  signal-idler  pair  in  which  the  signal  photon 
overlaps  spectrally  with  either  the  acceptance  bandwidth  of  the  converter  or  the 
detection  window.  In  [48],  we  showed  that  random  duty-cycle  errors  in  RPE  PPLN 
waveguides  substantially  increases  the  rate  of  parametric  fluorescence.  Fabrication 
tolerances  were  derived,  and  it  was  shown  that  for  typical  experimental  configurations, 
random  duty-cycle  errors  must  be  kept  to  a  level  ajl  <  1%  in  order  to  achieve  low- 
noise  performance.  Parametric  fluorescence  noise  was  measured  in  two  ways,  both 
directly  in  the  1.5-/mi  band  by  a  InGaAs  APD,  and  via  upconversion.  In  [75],  we 
also  described  a  measurement  of  random  duty-cycle  errors  in  PPLN  devices  via  direct 
imaging,  and  obtained  a  result  (cq/Z  ~  8%)  that  was  consistent  with  a  measurement 
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of  the  QPM  pedestal  via  a  careful  measurement  of  SHG  in  a  waveguide. 

The  second  noise  source  in  QFC  devices  is  spontaneous  Raman  scattering.  In 
Sec.  4.2,  we  presented  an  analysis  of  spontaneous  Raman  scattering  in  LiNbOs  sam¬ 
ples  and  waveguides.  The  Raman  spectrum  of  LiNbOs  was  found  to  decay  much  more 
slowly  at  high  frequency  shifts  than  would  be  expected  based  on  a  sum  of  Lorentzians, 
and  these  measurements  were  confirmed  by  noise  measurements  in  an  upconversion 
detector  [67].  However,  a  synthesis  of  the  analyses  of  both  noise  sources  indicates 
that  choosing  a  pump  wavelength  substantially  longer  than  the  signal  wavelength 
(or  equivalently,  up  <  u i)  is  highly  effective  at  reducing  the  impact  of  these  noise 
processes. 

The  primary  technological  interest  of  QFC  has  been  upconversion  single-photon 
detection.  Single-photon  detectors  for  the  1.5-/im  telecom  band  based  on  InGaAs 
APDs  suffer  severe  drawbacks  when  compared  with  Si  APDs  for  visible  and  near- 
visible  light.  Upconversion  of  1.5-/2111  signals  allows  detection  using  Si,  but  preliminary 
implementations  have  exhibited  noise  count  rates  well  above  105  counts/s  due  to 
spontaneous  scattering  [29,  30].  In  Chapter  5  and  in  [67],  we  demosntrated  lower- 
noise  single-photon  detection  using  a  long-wavelength  pump  at  Xp  ~  1850  nm.  Our 
system  exhibited  greater  than  99.99%  depletion  efficiency,  86%  net  internal  conversion 
efficiency  and  a  system  detection  efficiency  of  37%,  with  noise  counts  as  low  as  103 
counts/s.  The  spectral  dependence  of  the  noise  count  rate  was  found  to  be  consistent 
with  noise-photon  generation  by  spontaneous  anti-Stokes  Raman  scattering.  We  also 
demonstrated,  in  Sec.  5.4,  an  ultra-low-noise  upconversion  single-photon  detector  for 
signals  in  the  1.08-/mi  band,  which  are  important  for  proposed  quantum  information 
processing  schemes  based  on  neutral  donors  in  Si.  This  detector  used  a  1.55-/im-band 
pump  and  achieved  internal  conversion  efficiencies  of  81%  and  very  low  noise  count 
rates  of  approximately  30  counts/s,  due  to  to  the  very  large  frequency  difference 
between  the  signal  and  pump. 

While  it  had  been  realized  in  the  first  papers  on  QFC  that  quantum-state  down- 
conversion  was  possible,  these  devices  were  not  developed  until  much  more  recently. 
Most  candidate  systems  for  implementing  quantum  network  nodes  operate  with  op¬ 
tical  transitions  in  the  visible  or  near-visible  spectral  range,  therefore  requiring  a 
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frequency  conversion  technique  to  enable  long-haul  quantum  communications  based 
on  transmission  through  optical  fiber  at  1.5  /irn.  A  particularly  promising  quantum 
node  candidate  consists  of  an  electron  spin  in  an  InAs/GaAs  self-assembled  quantum 
dot.  In  Chapter  6  we  developed  a  donwconversion  quantum  interface,  whereby  the  sin¬ 
gle  photons  emitted  from  the  quantum  dot  at  A2  =  910  nm  were  mixed  with  a  strong 
pump  at  A p  =  2.2  /irn  to  produce  a  target  held  at  Ai  =  (A2  1  —  A'1)  1  =  1.56  /irn. 
We  generated  a  ps-pulsed  2. 2- /irn  source  based  on  DFG  in  a  bulk  PPLN  crystal, 
and  combined  the  2. 2- /irn  pulses  with  the  quantum-dot  single  photons  in  an  RPE 
PPLN  waveguide  for  downconversion.  The  resulting  single  photons  at  1.5 6- /irn  were 
detected  using  a  superconducting  single-photon  detector.  The  use  of  picosecond 
conversion  pulses  shaped  the  single-photon  emission,  as  observed  by  time-resolved 
measurements.  Due  to  the  large  frequency  difference  between  the  pump  and  target 
frequencies  in  this  experiment,  noise  counts  were  kept  to  a  level  below  approximately 
1  count/s.  A  measurement  of  the  intensity  cross-correlation  of  the  downconverted 
signal  showed  a  near  perfect  preservation  of  the  single-photon  character  of  the  quan¬ 
tum  dot  emission.  Additionally,  we  showed  that  the  downconversion  interface  was 
compatible  with  coherent  control  of  the  quantum  dot  spin  via  the  observation  of  Rabi 
oscillations.  These  characteristics  enabled,  to  our  knowledge  for  the  first  time,  gen¬ 
eration  and  verification  of  entanglement  between  the  solid-state  matter  qubit  (the 
quantum  dot  spin)  and  a  telecom-band  single  photon,  in  results  not  described  in  this 
dissertation.  We  believe  that  the  downconversion  technique  provides  a  route  towards 
long-distance  quantum  communications  based  on  quantum  dots. 

Chapter  7  described  two  novel  device  concepts  moving  past  the  one-step  up-  and 
downconversion  of  photons  usually  considered  to  be  the  limitations  of  quantum  fre¬ 
quency  conversion.  Since  with  RPE  PPLN  waveguides  it  is  possible  to  achieve  near¬ 
unity  conversion  efficiency  with  relatively  modest  pump  powers  of  ~  100  mW,  it  is 
also  possible  to  incorporate  two  conversion  steps  within  the  same  waveguide,  using 
two  QPM  gratings,  to  achieve  larger  frequency  excursions  while  still  maintaining  low- 
noise  operation.  In  the  first  demonstration  of  this  concept  [65],  we  used  a  1.8- /irn 
pump  to  upconvert  a  1.55-/nn  input  signal  to  the  green  spectral  range  in  a  two-stage 
conversion  process,  which  enabled  efficient,  low-noise,  and  low-timing  jitter  detection 


CHAPTER  8.  CONCLUSION 


124 


of  a  pulsed  1.5-/irn  signal  through  the  use  of  a  thin-junction  Si  APD  single-photon 
detector.  We  also  proposed  the  use  of  a  cascaded  frequency  converter  as  a  technique 
for  downconverting  single  photons  from  NV  centers  in  diamond  (at  a  wavelength 
of  A2  =  637  11111)  to  the  1.5  /1111  telecom  band  through  a  two-stage  downconversion 
technique.  As  low-noise  frequency  conversion  is  only  possible  with  long-wavelength 
pumping,  it  is  impossible  to  achieve  low-noise  conversion  when  >  2cui  in  a  single 
step.  The  use  of  cascaded  frequency  conversion  eliminates  this  constraint,  enabling 
possible  uses  of  this  technique  even  for  extremely  large  frequency  translations,  for  ex¬ 
ample  between  UV  photons  emitted  from  trapped-ion  qubits  and  the  1.5-/xm  telecom 
band. 

The  second  device  concept  was  multichannel  frequency  conversion.  In  this  config¬ 
uration,  multiple  pumps  can  be  used  to  either  probabilistically  or  deterministically 
convert  an  input  signal  to  one  of  several  output  frequencies.  As  discussed  in  Sec.  7.2, 
we  demonstrated  the  concept  with  a  two-channel  upconverter  for  converting  1.3-ywm 
single-photon  signals  to  two  wavelengths  near  710  mu.  Combining  the  two  channel 
converter  with  alternating  pumps  and  a  spectral  filter  at  the  output,  we  demonstrated 
a  multi- wavelength  sampling  experiment  whereby  the  timing-jitter-induced  maximum 
clock  rate  of  a  single-photon  communication  system  could  be  bypassed  by  the  effec¬ 
tive  time-to-space  conversion  provided  by  the  multichannel  device.  Additionally,  we 
used  a  volume  Bragg  grating  to  filter  the  upconverted  emission.  The  narrow  band¬ 
width  (0.04  nm)  of  the  VBG  compared  with  the  acceptance  bandwidth  of  the  device 
(0.3  nm)  enabled  a  strong  suppression  of  noise  counts  due  to  anti-Stokes  Raman 
scattering:  we  measured  a  noise  count  rate  of  approximately  500  counts/s  [144],  a 
factor  of  30  lower  than  previously  observed  with  these  wavelengths  [30] .  We  also  pro¬ 
posed  a  technique  to  extend  the  discrete-time  time-to-space  converter  to  continuous 
time,  potentially  enabling  very  high-bandwidth  quantum  communications  through 
time-division  multiplexing  or  pulse-position  modulation. 

My  initial  projects  in  the  Fejer  group  were  focused  on  modeling  a  synchronously 
pumped  OPO  for  terahertz  generation.  A  ps-pulsed  1064-nm  laser  pumped  a  type-II 
PPLN  crystal  producing  signal  and  idler  beams  spaced  by  the  desired  THz  frequency. 
The  signal  and  idler  then  interacted  to  produce  the  THz  radiation  in  an  QPM  GaAs 
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crystal  placed  inside  the  OPO  cavity.  The  modeling  involved  both  x*'2'*  and 
nonlinear  effects,  dispersion,  cavity  phase,  and  multiple  wave-mixing.  The  measured 
properties  of  the  intracavity  pulses  were  found  to  agree  well  with  the  model,  and  en¬ 
abled  a  more  detailed  understanding  of  possible  improvements  to  the  THz  generation 
efficiency  [155].  I  additionally  assisted  in  the  modeling  and  experimental  demonstra¬ 
tion  of  supercontinuum  generation  in  QPM  waveguides  [90,  91],  further  details  of 
which  can  be  found  in  the  dissertation  of  Chris  Phillips  [156]. 

8.2  Outlook  and  future  directions 

The  results  obtained  during  the  course  of  this  research  suggest  a  number  of  potential 
future  directions.  As  quantum  communication  techniques  become  more  sophisticaed 
as  the  field  develops,  one  imagines  that  the  same  sorts  of  functionality  that  x1'2'* 
mixers  bring  to  classical  communications  systems  will  find  analogous  applications  in 
the  quantum  realm.  These  classical  techniques  have  been  reviewed  in  [9],  and  here 
we  outline  a  few  exciting  future  opportunities. 

As  the  fabrication  model  for  APE  and,  subsequently,  RPE  waveguides  developed, 
it  becamse  possible  to  build  complex  integrated  optical  circuits  with  wavelength  con¬ 
version  functionality.  Linear-optics  quantum  computing  [98],  which  has  to  date  been 
attempted  largely  using  bulk  optics  with  limited  approaches  to  large-scale  scalable 
information  processing,  could  make  use  of  such  integrated  functionality  in  which  not 
only  the  interferometer  but  the  photon  source  (based  on  SPDC)  could  be  integrated 
onto  the  same  chip  [157].  The  ability  to  do  high-efficiency  fiber  coupling  with  RPE 
waveguides  also  facilitates  outcoupling  to  the  needed  single-photon  detectors. 

The  possibilities  outlined  in  the  general  idea  of  quantum  waveform  conversion  of¬ 
fer  intriguing  possibilities  when  combined  with  the  quantum  dot  single  photon  source 
[158].  Quantum  dot  photons  are  long  (~  ns)  on  communications  timescales:  mixing 
the  single-photon  with  a  chirped  pump  to  increase  its  bandwidth,  followed  by  com¬ 
pression  using  dispersion,  can  create  shorter  single-photon  wavepackets  which  could 
then  be  used  to  increase  the  transmission  bandwidth  of  quantum  communications 
systems. 
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Lastly,  all  the  work  described  in  this  dissertation  has  been  done  in  the  LiNbC>3  ma¬ 
terial  system.  QPM  devices  based  on  orientation-patterned  (OP)  zincblcnde  semicon¬ 
ductors  have  made  tremendous  progress  in  recent  years  [159]  .  In  addition  to  the  high 
nonlinearity  of  orientation  patterned  semiconductors,  a  critical  point  is  that  they  lack 
birefriengence.  As  such,  birefringent  phasematching  is  impossible  in  semiconductors, 
but  the  lack  of  birefringence  opens  up  intriguing  possibilities  for  polarization-diverse 
entangled-photon  generation  and  quantum  frequency  conversion  if  QPM  is  possible. 
As  a  classical  example,  optical  parametric  oscillation  with  a  depolarized  pump  has 
been  demonstrated  [160].  It  would  surely  be  interesting  to  study  the  polarization 
properties  of  spontaneous  parametric  downconversion  from  OP  semiconductor  de¬ 
vices.  With  the  recent  development  of  OP-GaP  [161],  such  applications  now  become 
feasible  at  customary  wavelengths  for  quantum  optics  in  the  visible  and  near-infrared 
where  good  single-photon  detectors  exist. 

Additionally,  as  ferroelectrics  can  only  achieve  QPM  for  one  polarization  at  a  time, 
polarization  diversity  has  typically  been  achieved  by  using  two  crystals  back-to-back 
[36].  With  semiconductor-based  QFC  devices,  both  polarizations  can  be  converted 
within  the  same  device,  which  should  enable  preservation  of  polarization  entanglement 
without  a  need  for  two  devices.  As  material  and  device  quality  continues  to  improve, 
orientation-patterned  semiconductors  will  be  an  interesting  platform  for  quantum 
frequency  conversion. 


Appendix  A 

Efficiency  pedestal  in  QPM  devices 
with  random  duty-cycle  errors 


In  this  Appendix,  we  calculate  the  ensemble-averaged  statistical  properties  of  the 
spatial  frequency  pedestal  that  arises  due  to  RDC  errors.  We  consider  a  three-wave 
difference-frequency  generation  (DFG)  process  involving  continuous,  plane-waves  with 
negligible  signal  gain  or  pump  depletion.  This  formulation,  combined  with  an  input 
noise  spectrum,  can  provide  a  semiclassical  description  of  spontaneous  parametric 
fluorescence  as  discussed  in  Chapter  4.  Our  results  are  also  directly  applicable  to 
second  harmonic  generation  (SHG)  when  the  first  harmonic  is  undepleted. 

Under  the  above  assumptions,  the  idler  at  the  end  of  the  QPM  grating  is  related 
to  the  signal  input  by  a  transfer  function  which  is  related  to  the  spatial  Fourier 
transform  of  the  grating  [162],  The  propagation  equation  for  the  generated  idler 
envelope  is  given  by  [12] 


dAj 

dz 


A*  ApC~iAk'z, 


rite 


(A.l) 


where  subscripts  ( i,s,p )  correspond  to  the  idler,  signal,  and  pump  respectively,  ujj  is 
the  carrier  frequency  of  wave  j,  and  rij  is  the  refractive  index  evaluated  at  uir  The 
envelopes  Aj  are  defined  in  relation  to  the  electric  field  as  E  —  \  J2jAjexp[i(ujjt  — 
kjZ )]  +  c.c.  where  summation  is  performed  over  j  G  {i,  s,  p },  and  kj  =  LOjiij/c.  The 
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phase  mismatch  is  defined  as  A k'  =  kp  —  ks  —  ki,  and  d(z )  =  ±d0  is  the  value  of  the 
relevant  second-order  nonlinear  coefficient  based  on  the  x ^  tensor  and  the  polariza¬ 
tion  directions  of  the  interacting  fields,  and  is  inverted  periodically  or  aperiodically 
in  the  QPM  grating. 

The  idler  held  Ai  can  be  written  in  a  form  normalized  to  its  value  at  the  output  of 
an  ideal  (phasematched,  with  no  RDC  errors)  periodic  QPM  grating  (for  which  the 
domain  size  would  be  given  by  l  =  7r / Ak'): 


4 


.  2  ujid0 

—i - 

7T  71  jC 


A*ApLai, 


(A.2) 


where  a*  is  the  normalized  idler  held  and  L  is  the  length  of  the  QPM  grating.  With 
this  definition,  Eq.  (A.l)  can  be  written  as 


dai  _  7 r  d(z)  iAk,z  _  ,.  iAkiz 

dz~2Ld0  9[Z) 


(A.3) 


where  the  normalized  grating  prohle  g(z)  is  defined  implicitly  in  Eq.  (A.3),  with  g(z)  = 
0  for  z  <  0  and  z  >  L.  By  integrating  Eq.  (A.3),  it  can  be  seen  that  the  normalized 
idler  output  a,(L)  is  simply  the  Fourier  transform  of  g(z),  g(k),  evaluated  at  spatial 
frequency  k  =  Ak.  The  square  magnitude  of  g  gives  the  idler  intensity  relative  to 
the  “maximum”  intensity  achievable  in  a  QPM  grating  of  the  same  length,  assuming 
that  the  pump  and  signal  (or  hrst  harmonic  in  the  case  of  SHG)  are  undepleted.  By 
writing  g(z)  as  a  sum  of  individual  domains  as  in  Eq.  (4.4),  g{k )  can  be  written  as 
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ikz[n] 


(A.4) 


where  z  is  a  vector  which  specihes  the  domain  boundary  positions.  The  subscript  z 
is  included  to  indicate  the  domain  boundary  vector  being  considered.  We  will  denote 
gz(k)  as  the  normalized  transfer  function. 

Next,  we  approximate  the  edge  terms  in  Eq.  (A.4),  i.e.  those  related  to  z [0] 
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and  z [N],  in  order  to  write  gz(k)  in  a  simpler  form  as  a  summation  of  complex 
exponentials.  We  note  that  z  [0]  and  z [N]  are  defined  by  the  start-  and  end-points 


of  the  device  and  are  therefore  not  random  variables.  The  resulting  (approximate) 
normalized  transfer  function  is  given  by 


N 


(A.5) 


71=1 


In  the  presence  of  random  perturbations  to  each  z[n],  this  summation  can  be  viewed  as 
a  random  walk  in  the  complex  plane.  The  ideal  vector  of  domain  boundary  positions 
(in  the  absence  of  any  RDC  errors)  is  denoted  z0;  for  ideal  periodic  poling,  z0[n]  =  nl 
for  domain  boundary  length  l.  For  a  phasematched  interaction  in  such  an  ideal 
grating,  k  =  A k  =  7r/Z,  in  which  case  the  summand  is  simply  unity  for  each  term, 
and  so  gZo(Tr/l )  =  1,  consistent  with  the  expected  result. 

Following  Ref.  [57],  we  assume  independent  and  identically  distributed  (HD)  er¬ 
rors  in  the  domain  lengths;  we  assume  a  normal  distribution  with  variance  erf  and 
mean  0.  A  schematic  of  RDC  errors  is  illustrated  in  Fig.  4.2.  In  most  cases,  RDC 
errors  are  introduced  during  electric-field  poling  following  lithographic  fabrication  of 
the  QPM  electrode  pattern,  such  that  their  distribution  is  fixed  for  a  given  device. 
This  assumption  implies  that  in  any  particular  device,  a  single  vector  z  applies.  How¬ 
ever,  in  order  to  understand  the  effects  of  RDC  errors  it  is  useful  to  perform  ensemble 
averaging  of  g(k),  treating  z  as  a  random  vector. 

The  ensemble-averaged  mean  of  the  normalized  transfer  function  is  given  by 


(A.6) 


where  the  ensemble-averaging  is  performed  over  random  vectors  z  and  the  relation 


—  g -k2crf/2eik(z) 


(A.7) 


which  is  a  Gaussian- integral  identity  for  scalar  z,  has  been  used.  From  Eq.  (A.6), 
the  expected  value  of  the  transfer  function  is  reduced  compared  to  that  of  the  ideal 
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transfer  function.  For  example,  the  phasematched  DFG  conversion  efficiency  averaged 
over  many  different  QPM  gratings  with  RDC  errors  would  be  reduced  by  a  factor 
exp  (-k2a2). 

In  addition  to  reducing  the  peak  efficiency,  RDC  errors  give  rise  to  additional 
spatial  frequency  components  far  from  nominal  phasematching.  To  see  this  behavior 
we  look  at  the  ensemble-averaged  value  of  \gz(k)\2,  which  is  given  by 


JV  N 

(ISzMI2)  =  (tt)  (££(-!) 


kL 
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e-ik(z0[n]-zQ[m]) 

n— 1  m=  1 


(A.8) 


where  the  second  relation  can  be  found  by  separating  the  first  relation  into  summation 
terms  with  m  —  n  and  those  with  m  —  n  are  independent  of  z[n]  and  hence 

do  not  acquire  an  exp[— (fccq)2]  factor.  We  can  write  Eq.  (A.8)  in  a  simpler  form  in 
terms  of  the  ideal  transfer  function  in  the  absence  of  RDC  errors  as 


(\9z(k)\2)  =  e  k2a'\gZ0(k)\2  +  N  (^r) 


(A.9) 


where  the  second  term  is  a  QPM  noise  pedestal  originating  from  the  RDC  errors. 
The  factor  N  associated  with  this  noise  floor  is  accurate  provided  that  there  are 
enough  domains  N  that  terms  0(1),  which  arise  due  to  the  neglected  edge  terms,  can 
be  neglected.  For  very  thin  crystals  with  small  N,  these  edge  terms  would  not  be 
negligible. 

For  periodic  QPM  gratings  with  a  50%  duty  cycle,  there  is  a  particularly  simple 
analytical  form  for  the  ideal  transfer  function  when  the  edge  terms  are  included. 
Assuming  that  z [n]  =  nl  and  defining  A k  =  k  —  7r/Z,  it  can  be  shown  that 


\~9*0(k)\2 
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The  first  factor  square  brackets  represents  the  reduced  Fourier  amplitudes  in  higher 
orders  of  the  QPM  grating  relative  to  to  the  first  order.  The  second  factor  is  a  “comb 
of  sines” ,  which  gives  rise  to  tuning  curves  at  each  odd  order  of  quasi-phasematching. 

Eq.  (A. 9),  when  combined  with  the  QPM  tuning  of  an  ideal  device  described  by 
Eq.  (A.  10),  gives  a  complete  description  of  the  effects  of  RDC  errors  for  all  orders  of 
quasi-phasematching,  and  is  useful  for  descriptions  of  highly  nonlinear  interactions 
involving  multiple  coupled  nonlinear  processes  [156].  In  this  dissertation,  we  have 
primarily  been  concerned  with  the  behavior  of  the  device  in  the  vicinity  of  first-order 
QPM.  With  these  assumptions,  Eq.  (A. 9)  can  be  substantially  simplified.  Near  first- 
order  QPM, 

\gZ0(k)\2  «  sine2  (A.ll) 

and  k  ~  n/l.  Therefore,  we  can  rewrite  Eq.  (A. 9)  as 

(&(Ai)|2}  =  e-'2".W  sine2  (PD)  +  b  (l  -  .  (A.12) 


This  equation,  presented  in  the  main  text  as  Eq.  (4.6)  gives  the  value  of  the  QPM 
pedestal  far  from  the  QPM  peak,  the  main  goal  of  this  Appendix. 

As  was  seen  in  the  simulation  results  of  Fig.  4.3,  an  individual  grating  exhibits 
large  fluctuations  in  the  generation  efficiency  far  from  phasematching,  as  might  be 
expected  for  a  random-walk-type  process.  Here,  we  calculate  the  variance  of  the 
efficiency  pedestal.  This  quantity  is  obtained  by  evaluating  both  the  second  and 
fourth-order  moments  (\gz(k)\2)  and  (\gz(k)\4)  using  the  same  approach  used  to  derive 
Eq.  (A. 8).  Specifically,  the  fourth-order  moment  is  given  by 


N- 1  N-l  N- 1  N- 1 
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Following  lengthy  algebraic  simplification,  the  fourth-moment  result  is  given  by 


i4>  =  ©4 


—  [  1  -  e~k2(r‘ 
N2  ' 


+  ^  -  e“4fc2°*  +  3e-3fc2ff2  -  6e~2k2(T‘  + 


(A.  14) 


When  N  is  large,  the  l/N2  term  will  dominate  over  the  1/N3  term.  The  variance  a 


of  the  normalized  transfer  function  is  calculated  as  a2  =  (\gz(k)\4)  —  (\gz(k)\2)2.  The 
standard  deviation  a„  can  then  be  approximated  as 


a',aN(uYh-e~*f 


(A. 15) 


This  approximation  applies  for  large  N  and  values  of  k  for  which  the  ideal  transfer 
function  |gZo(A;)|2  is  negligible.  This  relation  is  the  same  as  the  mean  (|^|2),  as  would 
be  expected  for  a  random  walk  process  for  the  generated  idler  field.  The  magnitude 
of  the  spatial  Fourier  spectrum  of  a  particular  grating  with  RDC  errors  is  thus  with 
high  probability  within  a  few  times  the  ensemble-averaged  mean. 


Appendix  B 

Effects  of  Cerenkov  modes  on 
parametric  interactions  in 
waveguides  with  random 
duty-cycle  errors 


We  consider  parametric  fluorescence  in  waveguides  with  QPM  random  duty-cycle 
errors  including  effects  of  radiation  modes.  In  planar  waveguides,  parametric  am¬ 
plification  involving  an  idler  mode  which  is  radiated  into  the  substrate  can  have 
substantially  larger  bandwidth  than  an  interaction  involving  only  bound  modes  [163] . 
Modeling  of  large-bandwidth  parametric  interactions  in  this  Cerenkov-idler  configu¬ 
ration  has  been  published  in  work  by  Rastogi  and  coworkers  [164],  For  simplicity, 
in  this  treatment  we  consider  planar  waveguides  as  in  previous  modeling.  This  as¬ 
sumption  is  adequate  to  describe  wave-mixing  processes  in  RPE  waveguides,  which 
typically  have  widths  suitable  for  lateral  confinement  but  modes  which  may  leak  into 
the  substrate.  The  analysis  is  in  any  case  straightforwardly  generalized  to  channel 
waveguides,  as  given  in  Eq.  (B.14).  The  growth  of  a  particular  idler  radiation  mode 
with  wavevector  /%  by  the  Cerenkov  process  is  given  as 

=  7(/y  av(z)a,(z)e-'^>  (B.l) 
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where  the  coupling  coefficient  7 (/3j)  and  nonlinear  polarization  V,  are  dehned  as 

eo  des 


7(A)  = 


up(x)u*s(x)u*(x)  dx 


ViU*(x )  dx 


(B.2a) 

(B.2b) 


where  we  see  that  the  coupling  coefficient  7(A)  can  be  viewed  as  the  projection  of 
the  nonlinear  polarization  onto  the  idler  radiation  mode  with  propagation  constant 
fa.  The  nonlinear  polarization  Ti(x)  has  been  defined  implicitly  in  Eq.  (B.2b). 

In  the  low-gain  regime  corresponding  to  parametric  fluorescence,  we  can  integrate 
Eq.  (B.l)  straightforwardly  over  the  length  L  of  the  crystal  since  the  signal  and  pump 
holds  are  unchanged  from  their  initial  values: 

-iA/3(/3i)L  _  1 

a* (A, L)  =  7(AK(°K(0)  _jA/3(ff.)  •  (B-3) 


The  photon  flux  summed  across  the  idler  modes  is  then  given  by 


Ni{L) 


ai{Pi,L) | 2  d ^ 


(B.4) 


where  /3min  and  /3max  are  the  minimum  an  maximum  values  of  the  range  of  radiation 
mode  propagation  constants.  Expanding  Eq.  (B.4)  gives 


Ni(L)  =  |op(0)as(0)|  — 


L 2 


7 (A)  sine 


AA(A  )L 


dA 


(B.5) 


We  are  interested  in  evaluating  the  integral  in  Eq.  (B.5)  in  two  different  regimes, 
corresponding  to  the  case  in  which  a  phasematching  feature  for  a  Cerenkov  idler 
mode  exists,  or  whether  such  a  phasematching  feature  does  not  exist,  but  there  is 
rather  a  weak  phasematching  pedestal  independent  of  Pi.  It  is  the  latter  case  that 
most  concerns  us  here,  but  the  calculation  of  the  former  is  illustrative. 

We  first  consider  the  case  in  which  a  phasematched  point  exists  in  the  range  of 
propagation  constants  between  /3min  and  Anax-  Assuming  that  the  crystal  is  long 
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enough  that  the  coupling  coefficient  7 (fa)  varies  slowly  over  the  width  of  the  phase¬ 
matching  feature  at  fa  =  (3pm,  we  can  pull  it  out  of  the  integral  over  fa  in  Eq.  (B.5), 
and  are  left  with 

Ni(L)  =  M<>H(0)|2  T  |7(ft»,)|2  (B.6) 

where  we  have  made  use  of  the  definite  integral 


(B.7) 


As  the  total  number  of  photons  generated  into  the  all  the  idler  modes  is  equal  to  that 
added  to  the  signal  mode,  we  can  write  the  signal  gain  as 
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Here,  as  in  most  treatments  of  Cerenkov-idler  parametric  amplihcation  [164],  the  gain 
is  dominated  by  only  the  idler  radiation  mode  that  is  phasematched. 

We  now  consider  the  opposite  limit,  in  which  there  are  no  phasematching  features 
over  the  range  of  idler  propagation  constants  (3min  to  /3max.  In  this  limit,  Eq.  (B.5) 
becomes 


Ni(L)  =  |ap(0)a£ 


29lo~r  /  Wi)\2dfa, 


(B.9) 


where  we  have  inserted  the  QPM  pedestal  height  gfa  of  Eq.  (4.9).  We  can  evaluate 
the  integral  of  Eq.  (B.9)  by  inserting  the  definition  of  the  coupling  constant  from 
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Eq.  (B.l).  Expanding,  we  have 


\i{fa)\2dfa  =  /  dfa 


Vi{x)u*{x)  da: 
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/oo  poo  p 

d x  Vi(x)  /  dx'V*(x')  /  d/^w-^x^/^x').  (B.lOd) 
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To  evaluate  the  integral  over  fa  in  Eq.  (B.lOd)  we  can  use  the  completeness  relation 
for  the  waveguide  eigenfunction  expansion: 

8{x—x')  =  y^uiin(x)u*n(x')  +  f  Ui(fa,x)u*(fa,x')  dfa,  (B.ll) 

n 

where  the  sum  is  over  the  bound  modes  of  the  waveguide  and  the  integral  is  over  the 
radiation  modes.  Applying  Eq.  (B.ll),  we  have: 


|t(A)|2  d/3j  =  /  d xVi{x)  /  dx'V*(x') 
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where  the  sum  extends  over  the  bound  modes  of  the  idler  wave.  Eq.  (B.12c)  follows 
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from  Eq.  (B.12b)  via  replacing  Vi(x)  with  its  definition  in  Eq.  (B.2b).  The  simplih- 
cation  in  Eq.  (B.12d)  is  valid  when  no  bound  modes  at  the  idler  wavelength  exist,  so 
the  sum  can  be  dropped. 

Therefore,  in  the  limit  when  any  phasematching  feature  can  be  neglected,  we  cal¬ 
culate  the  gain  (and  correspondingly  the  signal-photon  generation  rate)  via  inserting 
Eq.  (B.12d)  into  Eq.  (B.9): 


(B.13) 


Eq.  (B.13)  represents  all  photons  created  into  all  the  idler  modes,  and  has  now  been 


expressed  without  reference  to  the  modal  held  distributions  of  the  radiating  idler 


modes,  making  it  much  simpler  to  calculate.  The  integrand,  when  generalized  to  a 
2D  structure,  is  the  Cerenkov-idler  mode-overlap  integral: 


(B.14) 


where  we  define  Ai  as  the  interaction  area.  This  expression  is  used  to  calculate  the 
expected  SPDC  rate  and  the  factor  r  for  the  discussion  of  random  duty-cycle  error 
tolerances  from  Sec.  4.1. 
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